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Abstract 


The  most  widely  accepted  method  for  specifying  flying 
qualities  of  highly-augmented' aircraft  is  the  low-order  equivalent 
systems  technique.  This  technique  matches  the  frequency  response 
of  a  high-order  system  (HOS)  with  a  lower-order  equivalent  system 
(LOES).  The  LOES  is  generally  in  the  form  of  a  'classical^ 
(unaugmented)  aircraft  so  that  comparisons  with  the  classical  data 
base  can  be  made.  Most  LOES  research  has  been  accomplished  using 
frequency  response  matching,  however  there  are  also  methods 
available  for  matching  the  time  responses  of  the  HOS  and  LOES. 

Time  response  matching  is  an  attractive  option  for  flight  test 
applications  since  all  that  is  required  is  output  versus  input 
data . 

This  thesis  develops  and  tests  a  Least  Squares  LOES  time 
response  matching  program.  The  program  was  tested  analytically  at 
the  Air  Force  Institute  of  Technology,  Wright-Patterson  AFB,  OH. 
During  the  analytical  testing,  a  HOS  was  modeled  using  an  A-4D 
aircraft  with  servo-actuator  and  feel  system  dynamics.  Lower 
order  equivalent  parameters  were  calculated  using  the  Least 
Squares  program  and  the  results  were  compared  to  those  obtained 
using  a  frequency  response  matching  program  (LONFIT).  Various 
flight  conditions  and  feel  system  dynamics  were  used  to  change  the 
HOS  dynamics  and  the  effects  of  these  changes  on  the  two  matching 
techniques  were  investigated. ‘  ~"'ns 

The  Least  Squares  LOES  program  Vas  also  tested  during  a 
flight  test  program  at  the  USAF  Test 
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Pilot  School,  Edwards  AFB, 

'  kt)  r— 


CA.  During  this  test  a  variable  stability  NT-33A  aircraft 
operated  by  the  CALSPAN  Corporation  was  used  to  model  the  higher- 
order  aircraft  dynamics.  Thirteen  configurations  were  tested 
using  system  identification  and  offset  landing  tasks.  The 
configurations  were  chosen  to  span  Level  1,  2,  and  3  flying 
qualities  by  varying  wSp,  ?3p»  and  pre-filter  dynamics.  The  LOES 
Least  Squares  program  was  used  to  predict  the  flying  qualities  of 
the  configurations  based  on  MIL-STD-1797  guidance  for  equivalent 
£Sp,  w3p.  and  tq.  Parameter  identification  was  also  accomplished 
using  the  actual  flight  test  output  versus  input  data.  The 
results  of  the  flying  qualities  predictions  were  compared  to 
predictions  given  by  Hoh’s  Bandwidth  method  for  specifying  flying 


qualities. 


DEVELOPMENT  OF  A  LEAST  SQUARES  TIME  RESPONSE 


LOWER-ORDER  EQUIVALENT  SYSTEMS  TECHNIQUE 

I.  Introduction 

Background 

The  use  of  complex  augmentation  systems  for  modern  airplanes 
has  necessitated  the  development  of  various  methods  for  predicting 
and  analyzing  their  flying  qualities.  The  Military  Specification 
for  Flying  Qualities  of  Piloted  Airplanes,  MIL-F-8785B  (1)  was 
developed  in  the  1960’s  for  unaugmented  airplanes  and  did  not 
account  for  high-order  augmentation  systems.  As  a  result,  some 
modern  airplanes  were  designed  without  direct  benefit  of  the 
guidelines  of  MIL-F-8785B  since  these  guidelines  were  not 
considered  to  be  applicable.  MIL-F-8785C  (2),  published  in  1980, 
did  acknowledge  the  existence  of  high-order  systems  and  the  current 
specification,  MIL-STD-1797  (3)  includes  detailed  recommendations 
of  methods  for  specifying  the  flying  qualities  of  these  systems. 
MIL-STD-1797  states  that  the  preferred  method  for  specifying  flying 
qualities  of  highly  augmented  aircraft  is  the  low-order  equivalent 
systems  approach.  This  approach  makes  use  of  the  substantial  data 
base  that  exists  for  unaugmented  ’’classical”  aircraft  and  retains 
the  requirements  of  MIL-F-8785. 

Hodgkinson  (4)  demonstrates  the  need  for  equivalent  systems  by 
first  reviewing  conventional  methods  for  specifying  flying 
qualities.  In  the  past,  the  short  term  pitch  rate  dynamics  were 
typically  represented  by  a  linear  second  order  response  to  stick 
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force.  This  representation,  known  as  the  short  period 
approximation,  is  shown  in  equation  1. 

6  Kq ( s  +  1/Te2)  (1/Te2) 

-  =  -  =  -  (1) 

Fs  ( 21CgpWgpS  +  Wsp^ ^  ^sp’wsp^ 

Texts  such  as  McRuer  et  al.  (5)  define  all  of  the  terms  in 
this  function  using  aerodynamic  characteristics,  aircraft  speed, 
inertia,  etc.  This  expression  ignores  high-order  terms  such  as 
structural  modes,  stick  dynamics,  and  the  interaction  of  other 
aerodynamic  modes  of  motion  such  as  the  phugoid.  The  effects  of 
these  terms  were  considered  small,  as  were  the  effects  of 
linearization;  or  if  they  were  significant,  they  could  be 
considered  separately.  Thus,  this  transfer  function  defined  a 
well-accepted  "simulation  space”.  That  is,  the  parameters  Kq, 
1/t02»  ?ap<  and  w3p  provided  all  the  information  necessary  to 
simulate  the  short  period  response  of  the  aircraft.  Additionally, 
these  parameters  appear  (explicitly  or  implicitly)  in  MIL-F-8785 
and  in  MIL-STD-1797  meaning  that  they  also  define  a  well-accepted 
"specification  space”  for  short  term  pitch  dynamics.  For  past 
aircraft,  the  dimension  of  the  simulation  space  and  specification 
space  was  the  same. 

Transfer  functions  required  to  simulate  modern  aircraft 
dynamics  have  grown  in  complexity.  Figure  1,  for  example,  shows 
the  pitch  rate  transfer  function  which  emerged  during  the 
development  of  the  F-18  digital  flight  control  system.  Even  with 
the  same  simplifying  assumptions  used  in  equation  1,  the  simulation 
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space  required  to  describe  short-term  pitch  rate  dynamics  for  the 
F-1S  has  S9  dimensions!  Using  the  same  approach  as  in  the  past, 
this  means  that  a  89-dimensional  specification  space  would  be 
required  to  specify  the  longitudinal  short  period  dynamics! 
Clearly,  another  approach  is  needed. 


d  [-0. 5; 1273  [ -0. 3; 833  [-0.07;61]  [1;0.9]  £1;1]  Cl; 1. 23 


Fa  (0.5)  (0.8)  (0.9)  [  1 ;  1 . 1  ]  (2)  [1;2.4]  [0.7;4.1]  (3.4) 

[ 1 ; 2. 3]  C 1 ; 2. 8]  (4)  [0.1;53]  (13)  (14)  [0.4;34]  [0.34;46] 

x  — - - - - 

(3.7)  (6.7)  [0.2; 531  [1;131  C0.4;33]  [0.7;201  [0.3;421 

[0. 55; 39]  [1;22]  [0.7;37]  [0.9;34]  [0.9;34]  (35) 

x  - - - 

[0.9;  16]  [0.6;  27  ]  [0.3;66]  [0.4;85]  [0.4-.105]  [0.9;47] 

1 

x  - 

[0.3; 140]  [0.4; 180 ]  [0.7;102]  (72)  [0.5;228]  (112)  [0.6;279] 

1 

x  - 

[0.8; 329 ]  (259)  [0.9;378]  [1;410] 

x  TERMS  DUE  TO  DIGITAL  SYSTEM 

Figure  1.  Modern  High-Order  Response 
(Reference  4) 

Since  it  is  impractical  to  change  the  MIL  Standard 
requirements  to  correspond  to  the  complexity  of  these  high-order 
systems,  analytical  tools  have  emerged  so  that  high-order  systems 
can  be  specified  in  fewer  dimensions.  Several  methods  are 
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presented  in  MIL-STD-1797  including  the  low-order  equivalent 
systems  approach,  which  approximates  high-order  responses  with  low- 
order  responses;  the  Neal-Smith  method,  which  uses  a  mathematical 
model  of  the  pilot;  and  the  bandwidth  method,  which  takes  simple 
frequency  response  measurements.  If  the  gain  parameter  is 
specified  separately,  the  equivalent  system  of  equation  1  has  a 
three-dimensional  specification  space.  The  Neal-Smith  and 
bandwidth  methods  are  two-dimensional.  Human  pilots  process  these 
systems  to  produce  a  one-dimensional  specification  of  flying 
qualities  -  the  Cooper-Harper  rating. 

The  Low-Order  Equivalent  System  Method 

The  low-order  equivalent  sysf“”*  method  provides  an  order 
reduction  by  matching  (in  a  least  squares  sense)  the  frequency 
response  of  the  high-order  system  (HOS)  to  that  of  a  low-order 
equivalent  system  (LOES).  The  LOES  pitch  rate  and  normal  load 
factor  transfer  function  forms  to  be  used  for  the  matching  are 
given  in  MIL-STD-1797  as  follows: 

0  Kes(s  +  1/T01)(s  +  l/Te2)e_TS 

-  =  -  (2) 

Fa  (s2  2rpwps  +  wp2)(s2  +  2r3pwsps  +  wsp2) 

n2  Kn(s  +  l/Thl)a-TS 

-  =  -  (3) 

Fa  <s2  +  2Tpwps  +  wp2)(s2  +  2?spwaps  +  wsp2 ) 

These  expressions  are  linearized,  reduced-order  models  of  the 
actual  aircraft  response.  The  MIL-Standard  states  that  in  most 


4 


cases  the  phugoid  and  short  period  modes  are  sufficiently  separated 
that  further  order  reduction  is  possible  (i.e.  for  the  pitch  rate 
LOES)  as  follows: 


§  Kp( s  +  1/Tei) 

Phugoid:  -  =  - 

Fg  (s2  +  2rpwps  +  Up2) 


(4) 


e  Ksp(s  +  l/T02)e-TS 

Short  Period:  -  =  - 

Fg  (S2  +  2tgpUgpS  +  Wgp2 ) 


(5) 


Equations  4  and  5  are  universally  recognized  as  pitch  rate 
models  of  phugoid  and  short  period  dynamics  and  may  normally  be 
used  in  place  of  the  three-degree-of-f reedom  equations  over  their 
appropriate  respective  frequency  ranges.  The  normal  load  factor 
LOES  can  be  reduced  in  the  same  manner.  General  guidelines 
regarding  when  to  use  the  two-degree-of-f reedom  versus  the  three- 
degree-of-f  reedom  equations  are  discussed  in  more  detail  in  chapter 
II. 

Computer  programs  accomplish  the  frequency  response  matching 
by  varying  the  LOES  parameters  so  as  to  minimize  an  error  function 
between  the  HOS  and  LOES  over  the  frequency  range  of  interest. 

This  frequency  range  is  defined  by  MIL-STD-1797  as  follows: 


u^  -  Normally  0.1  rad/sec  but  >  Up 

ug  -  Normally  10.0  rad/sec  but  >  the 

resulting  equivalent  usp  and  1/T02 

and,  the  general  form  of  the  error  function  is: 
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M  =  S  (AG)2  +  K  S  (A*)2 

*  2  <SHOS  "  gL0ES*2  +  K  2  <*HOS  ~  *L0ES>2 

where:  M  is  the  sum-of-squares  "mismatch”  function 

G  is  the  amplitude  in  decibels 
#  is  the  phase  angle  in  radians 
K  is  a  weighting  factor 

AG  and  A$  are  calculated  at  discrete  frequencies  between 
and  ui2  evenly  spaced  on  a  logarithmic  scale.  Note  that  the 
mismatch  function  is  not  well  defined  unless  the  number  of  discrete 
frequencies  is  specified.  For  this  reason,  20  frequencies  are 
generally  used,  or  else  M  is  normalized  by  20/n  where  n  is  an 
arbitrary  number  of  frequencies  selected  for  the  match  (6).  This 
will  be  seen  in  chapter  III  where  the  error  function  for  the 
computer  programs  LONFIT  and  MACFIT  is  defined.  The  weighting 
factor  K  is  generally  specified  at  approximately  0.017.  This  gives 
approximately  the  same  importance  to  one  decibel  of  gain  mismatch 
as  to  eight  degrees  of  phase  mismatch  (7). 

Mismatch  Criteria 

Because  low-order  equivalent  systems  are  only  an  approximation 
to  the  actual  high-order  system  there  is  a  need  to  define 
acceptable  versus  unacceptable  levels  of  mismatch  between  the  two 
(i.e.  how  "good"  does  the  match  have  to  be  for  the  equivalent 
system  to  provide  an  acceptable  approximation?)  The  need  to  define 
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acceptable  mismatch  criteria  is  often  cited  as  a  fundamental 
drawback  of  the  equivalent  system  technique.  Other  criteria  such 
as  the  Neal-Smith  and  bandwidth  methods  simply  use  the  response  of 
the  high-order  system  and  do  not  involve  a  matching  process. 
Therefore  it  is  said  that  the  question  of  mismatch  does  not  arise, 
which  is  an  advantage.  Hodgkinson  (4)  argues  that  any  method  for 
specifying  high-order  systems  involves  a  reduction  in  dimension 
which  means  that  some  high-order  effects  are  ignored.  These 
effects  are  quantified  as  "mismatch”  for  low-order  equivalent 
systems,  but  can  be  similarly  quantified  for  the  Neal-Smith  and 
bandwidth  methods  using  the  concept  of  a  "minimum  order  system. ” 

In  fact,  the  more  a  specification  method  reduces  the  dimension 
(i.e.  the  lower  the  order  of  its  minimum  order  system)  the  larger 
the  mismatch  can  be. 

Regarding  the  specification  of  acceptable  mismatch  criteria, 

MIL-STD-1797  states  the  following: 

There  is  currently  insufficient  data  to  place  definitive 
requirements  on  mismatch  between  the  HOS  and  LOES.  It  should 
be  noted,  however,  that  the  question  of  mismatch  is  inherent 
in  any  n-dimensional  specification  of  an  m-dimensional 
response,  when  n  <  m.  (3) 

Hodgkinson  and  Johnston  (8)  addressed  the  question  of  mismatch 
with  the  USAF/CALSPAN  NT-33  aircraft  in  an  investigation  of 
longitudinal  and  lateral  dynamics  during  landing.  The  flying 
qualities  of  various  high-order  systems  and  their  low-order 
equivalents  were  compared.  Pilot  ratings  were  used  to  determine 
how  closely  the  LOES  must  approximate  the  HOS  to  be  a  valid  flying 
qualities  prediction  tool.  It  was  found  that  rather  large 
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analytical  mismatch  values  proved  unnoticeable  to  the  pilots.  The 
report  concludes  with  a  tentative  guideline  of  a  mismatch  value  of 
200  or  less  as  a  LOES  criterion  for  representative  augmented 
dynamics.  This  guideline  is  also  quoted  in  MIL-STD-1797.  This 
guideline  is  considerably  more  lenient  than  the  previous  "limit”  of 
10  which  was  based  solely  on  the  visual  appearance  of  the  frequency 
response  mismatch. 

Another  study  used  by  the  MIL-Standard  to  define  acceptable 
mismatch  criteria  was  conducted  by  Hodgkinson  and  Wood  (9).  These 
investigations  compared  pilot  rating  differences  between  pairs  of 
configurations  in  previous  NT-33  experiments.  Each  pair  of 
configurations  consisted  of  an  unaugmented,  low-order  aircraft 
response  and  a  high-order  system  formed  by  adding  terms  to  the  low  - 
order  response.  From  these  comparisons,  mismatch  envelopes  of 
"maximum  unnoticeable  added  dynamics”  were  derived  (figure  2). 
MIL-STD-1797  states  that  mismatches  between  the  HOS  and  the  LOES  in 
excess  of  the  values  shown  in  figure  2  would  be  cause  to  suspect 
that  the  equivalent  parameters  may  not  accurately  predict  pilot 
opinion. 

Equivalent  Time  Delay 

Equivalent  system  studies  by  DiFranco,  Neal-Smith,  and  MCAIR 
(10,11,12)  determined  that  the  specification  space  of  classical 
systems  is  insufficient  to  adequately  model  most  high-order 
systems.  This  is  true  primarily  because  the  phase  lags  of  the 
high-frequency  modes  are  not  accounted  for.  Since  the  major  effect 
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of  these  high-frequency  modes  (i.e.  frequencies  well  beyond  the 
crossover  frequency)  is  an  equivalent  time  delay,  the  solution  to 
the  problem  is  to  add  an  equivalent  time  delay  term  to  the  low- 
order  system.  (This  is  the  e-TS  term  shown  in  equations  2,3,  and 
5). 


Frequency  (rod/tec) 


Figure  2.  Envelopes  of  Maximum  Unnoticeable  Added  Dynamics 
(Reference  3) 


Equivalent  time  delay  can  be  used  to  model  the  effects  of  many 
different  flight  control  components  such  as  stick  prefilters,  feel 
systems,  actuators,  and  structural  filters.  The  DiFranco  and  Neal- 
Saith  investigations  (10,11)  found  that  these  time-delays  - 
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consistently  degraded  pilot  ratings  down.  Time  delay  therefore  has 
a  significant  effect  on  flying  qualities  and  must  be  included  in 
any  viable  equivalent  system  model.  This  adds  one  extra  dimension 
to  the  specification  space  in  the  MIL-Standards . 

Problems  Encountered  in  LOES  Methods 

The  encouraging  results  from  equivalent  systems  techniques 
hold  much  promise  for  their  ability  to  specify  flying  qualities  of 
highly  sophisticated  aircraft.  There  are,  however,  a  few 
unresolved  questions  to  be  addressed  concerning  the  LOES  approach. 
The  most  significant  question  is  on  the  definition  of  ’’equivalence” 
itself.  In  some  cases  it  has  been  shown  that  the  LOES  parameters 
are  not  necessarily  equivalent  to  their  classical  counterparts. 
Before  discussing  this  point  further,  it  will  be  helpful  to  review 
the  MIL-Standard  requirements  for  longitudinal  equivalent  system 
parameters. 

The  requirements  on  equivalent  pitch  time  delay  (T0>  and  on 
equivalent  short  period  damping  (Eap)  are  straightforward  and  are 
specified  as  follows  (3): 

Equivalent  Time  Delay  and  Short  Period  Damping  Ratio  Criteria: 

Level  Allowable  Delay  Allowable  Limits 

1  0.10  (sec)  0.35-1.30 

2  0.20  (sec)  0.25-0.35  or  1.30-2.00 

3  0.25  (sec)  <  0.25  or  >  2.00 

The  requirements  on  equivalent  short  period  natural  frequency 
(uSp)  are  given  in  terms  of  the  control  anticipation  parameter 
(CAP)  defined  as  CAP  *  wap2/(n/a)  where  n/a.  =  ( V/g ) ( 1/Tq2) •  CAP  is 
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approximately  equal  to  the  ratio  of  initial  pitching  acceleration 
to  steady  state  normal  acceleration.  Because  of  the  time  lapse 
before  reaching  steady  state,  the  pilot  needs  an  earlier  indication 
of  the  response  to  control  inputs.  Both  the  initial  and  final 
responses  must  be  neither  too  sensitive  nor  too  insensitive  to  the 
commanded  flight-path  changes.  The  requirements  on  CAP  are  listed 
below  (3): 

Equivalent  Short  Period  Natural  Frequency  Criteria: 

Level  Allowable  CAP  Limits  (g~^sec-1) 

1  0.28-3.60 

2  0.16-0.28  or  3.60-10.00 

3  <  0. 16  or  >  10. 00 

The  ta p  and  CAP  requirements  shown  above  are  for  category  A 
(rapid  maneuvering,  precision  tracking,  or  precise  flight-path 
control)  flight  phases  only.  For  a  more  complete  discussion  of 
these  requirements,  the  interested  reader  is  referred  to  MIL-F- 
8785C  (2)  or  to  MIL-STD-1797  (3). 

The  available  data  base  of  H0S  flying  qualities  piloted 
evaluations  is  currently  quite  limited.  The  Neal-Smith  (11)  and 
Landing  Approach  High-Order  System  ”LAH0S”  (13)  studies  are  the 
most  widely  recognized.  Fortunately,  this  small  quantity  of  data 

covers  quite  a  large  range  of  flight  conditions  and  tasks.  The 
Neal-Smith  investigation  provides  good  examples  of  both  the 
encouraging  results  from  the  L0ES  method  as  well  as  problems  that 
have  been  encountered.  Of  the  51  configurations  flown  in  the  Neal- 
Smith  program,  41  were  either  unaugmented  or  included  only  first  or 
second  order  lags.  The  flying  qualities  of  28-  of  these  were. 
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correctly  predicted  by  the  LOES  method.  The  remaining  ten 
configurations  had  (possibly  unrealistic)  lead/lag  combinations, 
and  only  five  of  these  were  predicted  accurately.  For  the  five 
that  failed,  the  equivalent  dynamics  (£sp,  wsp,  1/Tq2*  and  tq  ) 
predicted  level  1  flying  qualities  but  were  rated  level  2  by  the 
pilots. 

There  are  two  potential  methods  for  dealing  with  lead/lag 
systems  like  those  in  the  Neal-Smith  study.  Unfortunately,  neither 
is  physically  very  appealing  and,  in  each,  there  is  an  underlying 
question  as  to  the  universality  of  the  equivalent  systems  approach. 

Redefine  Limits  on 

If  the  minimum  allowable  £sp  for  level  1  were  increased  from 
0.35  to  0.50,  four  of  the  five  configurations  would  meet  the 
requirements  (i.e.  predict  level  2).  But  applying  this  restriction 
to  unaugmented  aircraft  as  well  is  unappealing  since  the  lower  Tsp 
boundary  is  well  supported  by  flight  test  data  for  classical 
aircraft.  The  alternative  is  to  specify  two  sets  of  requirements  - 
one  for  unaugmented  aircraft  and  one  for  augmented  aircraft.  This 
is  especially  unattractive  since  it  is  an  admission  that  ”?e”  is 
not  equivalent  to  ?sp  and  "equivalent  systems”  is  a  misnomer. 
Additionally,  specifying  two  sets  of  requirements  presents  the 
problem  of  defining  the  specific  level  of  augmentation  at  which  the 
requirements  would  change  over.  For  example,  should  the  addition 
of  a  simple  high-frequency  stic1'  filter  (whose  only  major  effect  is 
to  increase  t)  suddenly  mean  that  the  airplane  must  meet  a  more 
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stringent  damping  requirement? 


Redefine  t 

Three  of  the  four  low  £gp  violators  discussed  above  also  had 
t  -  0  and  it  was  found  that  a  better  LOES  fit  could  have  been 
obtained  for  these  three  cases  if  t  were  allowed  to  be  less  than 
zero.  This  illustrates  another  fundamental  problem  with  the  LOES 
approach.  The  problem  is  in  interpreting  the  significance  of 
negative  time  delay  or  time  lead.  Physically,  it  might  be 
considered  to  represent  a  high-order  system  which  is  too  abrupt 
(i.e.  the  system  responds  t  seconds  before  the  input  is  made  or  has 
finite  magnitude  at  zero  time).  Unfortunately,  the  concept  of  a 
system  responding  prior  to  an  input  (i.e.  a  non-causal  system)  is 
not  very  appealing  physically  and  it  raises  some  question  as  to 
whether  or  not  the  equivalent  systems  approach  is  universally 
applicable. 

Another  important  question  regarding  the  LOES  approach  is 
whether  or  not  all  LOES  parameters  should  be  freed  in  the  matching 
process.  Specifically,  it  is  argued  that  the  pitch  numerator  term 
1/Tq2  (which  is  approximately  equal  to  La)  should  be  held  at  the 
basic  aircraft  value.  The  rationale  behind  keeping  it  fixed  is 
that  La  is  determined  by  wing  size  and  planform  which  are 
dimensions  that  can  not  be  varied  for  a  particular  aircraft.  In 
LOES  studies  where  La  has  been  freed,  however,  the  results  have 
produced  large  changes  in  La  from  the  basic  aircraft  values  (3). 
Again,  this  raises  questions  as  to  the  ’’equivalence”  of  the  LOES 
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parameters.  For  example,  how  is  the  equivalent  La  to  be 
interpreted  if  it  is  significantly  higher  (or  lower)  than  the  basic 
aircraft  value? 

It  should  be  mentioned  that  despite  the  problems  discussed  in 
this  section,  the  LOES  approach  to  the  analysis  of  highly  augmented 
aircraft  continues  to  show  promise.  The  success  obtained  with  HOS 
flight  tests  (11,13)  has  shown  that  the  pitch  short  period 
equivalent  dynamics  are  relatable  to  their  classical  counterparts 
for  most  cases.  An  additional  important  finding  resulting  from 
LOES  studies  is  the  significant  role  that  equivalent  time  delay  can 
play  in  the  degradation  of  both  longitudinal  and  lateral  flying 
qualities  (8,10,11.12). 

The  LOES  ?rs)roach,  as  defined  in  this  chapter,  is  accomplished 
by  minimizing  the  frequency  response  error  (equation  6)  between  the 
HOS  and  the  LOES.  Most  research  in  this  area  has  been  accomplished 
using  this  type  of  frequency  response  matching,  however  there  are 
also  methods  available  for  matching  the  time  response  of  the  HOS 
and  the  LOES.  In  certain  applications,  matching  the  time  response 
may  be  an  attractive  option  since  the  HOS  does  not  need  to  be 
modeled  in  the  frequency  domain.  Instead,  all  that  is  required  is 
the  time  history  of  HOS  input  and  response.  This  approach, 
therefore,  may  be  especially  well  suited  to  flight  test  system 
identification  requirements.  Shafer  (14),  for  example,  has 
conducted  investigations  using  maximum  likelihood  estimation  and 
the  input/output  time  histories  of  flight  test  data.  Her  studies 
used  data  from  the  highly  augmented  F-8  digital  fly-by-wire  (DFBW) 
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aircraft  to  create  a  LOES  model  of  the  aircraft’s  longitudinal 
system  response.  Shafer  reported  encouraging  results  in  that  the 
dynamics  of  the  HOS  were  excellently  matched  by  the  maximum 
likelihood  technique.  Further  investigations  using  this  technique 
are  in  progress. 

Another  possible  advantage  to  time  response  matching  is  that 
the  time  response  is  a  measure  of  the  real  world  behavior  (i.e.  it 
is  a  measure  of  how  the  aircraft  actually  ’’flies”).  Thus,  matching 
in  the  time  domain  would  apparently  ensure  a  LOES  that  ’’flies”  more 
like  the  HOS.  A  fundamental  question  in  this  regard  however,  is 
whether  or  not  a  good  match  in  the  time  domain  corresponds  to  a 
good  match  in  the  frequency  domain.  The  answer  to  this  question  is 
important  if  the  final  step  of  the  matching  technique  is  the 
extraction  of  frequency  response  data  (£sp,  wsp,  I/T92,  and  tq  )  in 
order  to  take  advantage  of  the  classical  data  base. 

The  purpose  of  this  thesis  is  to  develop  and  test  a  time 
response  matching  technique  adaptable  to  flight  test  parameter 
identification  requirements.  The  technique  will  be  compared  to  a 
frequency  response  matching  program,  Hoh’s  Bandwidth  method,  and 
actual  pilot  ratings.  Comparisons  between  the  techniques  will  be 
based  on  mathematical  analysis  such  as  fidelity  of  matching  and 
differences  between  LOES  parameters,  as  well  as  subjective  analysis 
such  as  problems  encountered  during  matching  and  adaptability  to 
flight  test  data  analysis  techniques. 
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II.  High-Order  Control  Systems 


It  has  been  known  for  some  time  that  control  system  dynamics, 
as  well  as  open-loop  airplane  dynamics,  effect  the  handling 
qualities  of  the  closed-loop  pilot  and  airplane  combination. 

Prior  to  the  late  1960’s,  however,  most  handling  qualities 
investigations  were  concerned  primarily  with  the  effects  of 
variations  in  certain  open-loop  airplane  parameters.  The  control 
system  characteristics  were,  in  general,  held  constant  and  their 
effect  on  airplane  response  to  pilot  inputs  was  not  investigated. 
By  making  the  control  system  sufficiently  ’’fast”  compared  to  the 
open-loop  airplane,  it  was  assumed  that  the  effects  of  the  control 
system  dynamics  could  be  neglected.  With  this  simplifying 
assumption,  the  handling  qualities  were  described  in  terms  of  the 
following  open-loop  airplane  parameters: 

(1)  The  short  period  frequency  (u)gp)  and  damping  ratio 


<  £sp) 

(2)  The  pitch  attitude  numerator  term  (La)  and  the 
response  ratio  (n/a) 

(3)  The  maneuvering  stick  force  gradient  (Fs/n) 

The  above  airplane  parameters  can  be  used  to  adequately 

describe  handling  qualities  if  the  effects  of  the  control  system 
are  not  too  significant.  In  recent  years,  however,  complex  flight 
control  systems  have  become  increasingly  common.  These  systems, 
with  various  types  of  feedback  and  feedforward  loops,  may 
introduce  significant  time  delay  and/or  additional  dynamic  modes 
which  have  natural  frequencies  close  to  the  airplane's  short 
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period  frequency.  When  this  is  the  case,  the  airplane’s  flying 
qualities  can  become  completely  unacceptable  even  if  the  short 
period  mode  itself  is  well  behaved. 

The  increased  complexity  of  flight  control  systems  has  been  a 
direct  result  of  increased  aircraft  performance.  Prior  to  the 
1940’s,  the  flight  and  loading  envelopes  of  most  airplanes  were 
limited  enough  that  stick  force  per  g  (Fs/n)  could  be  kept  within 
reasonable  limits.  The  control  systems  were  generally  low-inertia 
mechanical  systems,  with  restoring  and  damping  forces  derived  from 
aerodynamic  hinge  moments.  The  dynamics  of  these  control  systems 
were  quite  fast  compared  to  the  aircraft  dynamics,  and  there  was 
little  coupling  between  the  aircraft  and  control  system  modes. 
Friction  and  freeplay  could  cause  significant  control  system  lag 
for  small  control  inputs,  however,  resulting  in  control  problems 
during  precision  maneuvers. 

As  airplane  performance  increased,  it  became  difficult  to 
keep  maneuvering  forces  (Fs/n)  within  reasonable  limits.  The  high 
stick  forces  due  to  increased  elevator  hinge  moments  were  overcome 
by  the  introduction  of  aerodynamic  devices  such  as  spring  tabs, 
horn  balances,  and  set  back  control  surfaces.  In  order  to 
minimize  variations  in  Fs/n,  normal  acceleration  bobweights  were 
often  employed. 

As  performance  was  increased  even  further,  it  became 
necessary  to  resort  to  hydraulic  actuator  controls  to  overcome  the 
increased  hinge  moments.  The  use  of  hydraulic  powered  control 
introduced  another  control  system  problem:  there  was  no  feedback 
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of  control  surface  hinge  moments  to  the  stick  for  feel.  To  make 


the  aircraft  ’’feel  right”  to  the  pilots,  feel  springs  were  added. 
Pilots  normally  fly  by  physical  association  between  applied  force 
and  maneuvering  response.  At  low  dynamic  pressure,  the 
predominant  cue  is  pitch  rate  and,  at  high  dynamic  pressure  the 
predominant  cue  is  acceleration.  Without  feel  devices,  the 
maneuvering  response  cues  were  still  present,  but  the  stick 
’’flopped"  arounr  the  cockpit,  and  the  need  for  force  producers  at 
the  controls  became  critical.  Feel  springs  were  able  to  solve 
part  of  the  problem,  especially  at  low  dynamic  pressure,  and  the 
bobweight  was  able  to  provide  additional  force  at  the  control 
stick  during  high  g  maneuvering  flight. 

The  problem  with  feel  springs,  hydraulic  actuators, 
bobweights,  and  other  control  augmentation  systems  is  that  they 
may  introduce  unwanted  dynamics  into  the  flight  control  system. 
This  in  turn  may  lead  to  control  difficulties  during  maneuvering 
flight  and  in  some  cases  lead  to  severe  oscillations  of  the 
aircraft.  From  an  equivalent  systems  standpoint,  the  control 
system  dynamics  may  cause  difficulties  in  the  matching  process. 

The  aircraft  dynamics  may  be  changed  by  the  presence  of  the 
control  system  to  the  point  that  it  is  no  longer  possible  to  match 
the  system  to  conventional  low-order  forms.  These  matching 
difficulties  seem  to  be  compounded  when  the  control  system 
dynamics  are  ’’close”  to  the  aircraft  dynamics  and  there  is 
coupling  between  the  control  system  and  aircraft  roots. 

The  next  sections  of  this  chapter  develop  high-order  control 
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and  aircraft  systems  that  are  typical  of  systems  currently  in  use 
on  modern  high-performance  fighters.  This  will  be  accomplished 
using  feel  system  and  servo-actuator  dynamics  combined  with 
conventional  aircraft  dynamics  (A-4D).  Since  the  primary  purpose 
for  modeling  these  high-order  systems  is  to  compare  equivalent 
systems  methods,  the  intent  here  is  not  to  create  overly  complex 
systems.  There  are  certainly  other  causes  of  high-order  effects 
besides  feel  systems  and  servo-actuators  (i.e.  stability 
augmentation,  model  following,  structural  dynamics,  bobweights, 
etc.)  that  will  not  be  considered  here.  The  available  data  base 
for  HOS  flying  qualities  piloted  evaluations  is  currently  quite 
limited,  but  some  of  the  most  widely  recognized  programs  have  used 
a  feel  system  and  actuator  to  model  the  high-order  effects.  The 
development  presented  here  is  adapted  from  these  sources  (10,11). 

Control  System  Block  Diagram 

A  block  diagram  of  the  open-loop  airplane  control  through  a 
feel  system  and  actuator  is  shown  in  figure  3.  As  illustrated  in 
the  diagram,  control  is  initiated  by  the  pilot  applying  a  control 
force  (F3)  to  the  stick.  Through  the  dynamics  of  the  elevator 
feel  system,  this  force  results  in  a  stick  displacement  (8es) 
which  is  an  input  to  the  elevator  actuator.  The  elevator 
deflection  (8e)  is  a  function  of  the  actuator  dynamics,  and 
commands  the  response  of  the  airplane  as  determined  by  the 
airplane  dynamics. 
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Figure  3.  Block  Diagram  of  Open-Loop  Airplane  Control 

The  overall  aircraft  response  to  pilot-applied  stick-force 
inputs  can  be  represented  as  a  transfer  function  which  is  the 
product  of  the  transfer  function  of  the  feel  system,  actuator,  and 
airplane.  Thus,  the  pitch  rate  (0)  and  normal  acceleration  (nz) 
transfer  functions  for  stick  force  inputs  (Fa)  can  be  expressed  as 
fol lows: 


0(s)  0<s)  Se(s)  Bes(s) 

-  _  -  x  -  x  - 

Fa(s)  6e(s)  Se3(s)  F3(s) 

nz(s)  nz(s)  8e(s)  8ea(s) 

-  =  -  x  -  x  - 

F3(s)  Se( s )  Se3(s)  F3( s ) 


(7) 


(8) 


The  parameters  of  each  of  the  elements  (feel  system,  actuator,  and 
airplane)  in  the  above  transfer  functions  can  be  varied 
independently  during  ground  or  in-flight  simulations  to  assess  the 
effects  on  handling  qualities  of  various  combinations  of  these 
elements.  A  detailed  description  of  each  of  these  elements  is 
given  in  the  following  sections. 
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Feel  System  Dynamics 


The  feel  system  dynamics  can  be  approximated  reasonably  well 
with  a  ’’spring-mass-damper  system”  as  shown  by  Biezad  (16).  This 
model  results  in  the  following  second-order  transfer  function: 


_ KfsU)fs2 _ 

s2  +  2rfau)fas  +  wfa2 


(9) 


The  natural  frequency  and  damping  ratio  (w£a  and  Tfa) 
determine  the  dynamics  of  the  feel  system  and  can  be  varied 
independently  of  the  aircraft  dynamics  to  evaluate  their  effects 
on  aircraft  response.  The  natural  frequency  uj£S  can  be  varied 
from  approximately  6  to  31  rad/sec  which  effectively  selects 
’’slow”  to  ’’fast”  feel  system  response  times.  During  the  matching 
process,  uifa  was  varied  from  6  rad/sec  (’’close”  to  the  aircraft 
roots)  to  31  rad/sec  ("far”  from  the  aircraft  roots).  This  was 
done  to  investigate  the  effects  of  control  system/aircraft 
coupling  on  equivalent  systems  techniques. 

The  damping  ratio  ZfB  was  varied  from  0.45  to  1.00  to 
determine  values  which  would  provide  the  desired  HOS  time  response 
characteristics.  With  the  ’’slow”  feel  system,  it  was  found  that  a 
value  of  tfB  -  1.00  resulted  in  the  greatest  coupling  between  the 
control  system  and  aircraft  roots.  The  ’’moderate”  feel  system, 
and  t£B  -  1.00,  resulted  in  a  relatively  uncoupled  (pure  time 
delay)  HOS  response.  Since  the  intent  was  not  to  evaluate  changes 
in  the  damping  ratio  from  one  configuration  to  another,  and 
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because  tfa  =  1.00  was  found  adequate  to  provide  the  desired 
coupled  and  uncoupled  response  characteristics,  tfa  was  held 
constant  at  1.00  for  all  HOS  configurations. 

A  complete  specification  of  the  feel  system  requires 
specification  of  the  steady-state  spring  gradient  (Fa/Sea)aa. 
Applying  the  final  value  theorem  for  Laplace  transforms  to  eq.  9 
gives  the  steady-state  spring  gradient  for  a  step  input  as 
(Fa/Ses)aa  =  1/Kf 3.  During  the  matching  process,  the  feel  system 
i jin  was  held  constant  at  K^g  =  0.046  which  corresponds  to  a 
steady-state  spring  rate  of  21.7  lbs/in.  This  is  a  typical  value 
of  the  spring  gradient  for  fighter  "up  and  away”  evaluations  (11). 

The  variation  in  feel  system  response  time  is  illustrated  in 
figure  4  which  compares  "moderate”  (Wf3  =  18.5  rad/sec)  and  "slow” 
(w£a  =  6  rad/sec)  feel  systems.  In  each  case,  the  response  is 
shown  for  a  step  input  with  Zf3  ~  1.00  and  Kfa  =  0.046.  To 
demonstrate  the  spring  gradient  calibration,  a  step  input  of  Fa  = 
21.7  lbs  is  used  which  results  in  the  correct  steady-state  output 
of  (8ea)aa  =  1.00  in. 


Servo-Actuator  Dynamics 

The  servo-actuator  dynamics  can  be  approximated  by  a  first- 
order  "lag”  response  as  shown  by  Biezad  (15)  and  Roskam  (16).  The 
transfer  function  used  to  represent  the  servo-actuator  is  shown 
below: 


K 


gr 


Tas  +  1 


(10) 
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DELTA-ES  (IN)  DELTA-ES  (IN) 


TIME  (SEC) 

(b)  Slow  Feel  System  Response  (uj£s  =  6.0  rad/s) 


Figure  4.  Feel  System  Response  Characteristics 

(step  input  Fa  =  21.7  lbs,  £fa  =  1.00, 
and  Kfa  =  0.046  in/lb) 
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The  time  constant  Ta  determines  the  speed  of  the  actuator 
response.  Although  a  value  of  1/Ta  =  20  is  often  used,  this  may 
be  an  optimistic  view  for  the  dynamic  performance  of  most 
actuators.  In  a  study  by  Gentry  and  Pujara  (17),  for  example,  an 
elevator  actuator  from  the  YF-16  was  modeled  using  an  equivalent 
system  method  over  a  frequency  range  of  0.1  to  10.0  rad/sec  (as 
suggested  for  use  with  MIL-F-8785C ) .  This  equivalent  system 
approach  resulted  in  a  more  conservative  and  realistic  first-order 
model  of  13/(s  +  13).  For  the  matching  techniques  used  in  this 
thesis,  therefore,  the  servo-actuator  time  constant  was  fixed  at 
1/Ta  =  13.  Variations  in  the  control  system  roots  were 
accomplished  only  by  changing  the  feel  system’s  natural  frequency 
(ufs) . 

Elevator-to-Stick-Force  Gearing 

The  gain  Kgr  shown  in  equation  10  represents  the  elevator-to- 
stick-force  gearing  selection  which  determines  the  stick  force  per 
g  (F3/n)  for  a  particular  configuration.  In  flying  qualities 
evaluations,  the  selection  of  elevator-to-stick-force  gearing  can 
be  difficult  for  the  pilot  because  it  requires  a  compromise 
between  the  conflicting  demands  of  satisfactory  initial  forces  for 
good  precision  tracking  capability,  and  satisfactory  steady  forces 
for  gross  fighter  maneuverability.  Consider,  for  example, 
configurations  with  initial  response  characteristics  that  are 
described  as  sluggish  or  slow,  either  because  of  low  wap  or 
control  system  lags.  In  these  configurations,  the  pilot  will  use 
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large  initial  inputs  to  "overdrive  the  airplane”  in  order  to 
achieve  the  desired  response.  This  calls  for  a  high  elevator— to— 
stick-force  gearing  to  keep  the  initial  forces  reasonable.  These 
same  configurations,  however,  may  have  an  unpredictable  final 
response  because  the  airplane  tends  to  "dig  in”.  These 
characteristics  demand  low  elevator  gearing  to  provide  adequate  g 
protection  during  the  gross  fighter  maneuvers;  and  thus,  a 
compromise  gearing  selection  is  required.  Conf igurations  with 
initial  response  characteristics  described  as  abrupt  or  too 
sensitive,  on  the  other  hand,  require  the  opposite  gearing 
compromise.  Because  the  initial  forces  for  these  configurations 
are  considered  "too  light”  by  the  pilot,  a  low  elevator  gearing  is 
required  to  prevent  inadvertent  inputs.  Such  a  gearing  selection, 
however,  tends  to  produce  steady  forces  that  are  too  high  for  the 
gross  fighter  maneuvers;  and  thus,  the  opposite  compromise  is 
required. 

In  the  Neal-Smith  evaluations  (11),  the  pilots,  when  faced 
with  a  compromise  in  the  elevator  gearing  selection,  were  not 
willing  to  vary  F3/n  over  a  very  large  range.  In  general,  the 
pilots  would  not  compromise  their  ability  to  pull  large  load 
factors  (i.e.  achieve  satisfactory  gross  maneuverability)  even  if 
the  resulting  Fs/n  was  not  compatible  with  the  initial  forces 
required  for  precision  tracking.  The  values  of  Fa/n  selected  by 
the  pilots  in  the  Neal-Smith  study  ranged  from  approximately  4.5 
to  7  lbs/g  although  a  few  excursions  (as  low  as  3  and  as  high  as 
14  lbs/g)  were  noted  in  some  parts  of  the  program. 
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Based  on  the  findings  of  the  Neal-Smith  report,  the  gain  Kgr 
was  selected  to  give  a  stick  force  per  g  of  6  Ibs/g  during  the 
equivalent  systems  matching.  To  find  the  required  value  of  Kgr, 
the  stick  force  per  g  was  calculated  using  the  short  period 
approximation  for  nz/Fs  given  by  eq.  11. 

nz  Kfswfs^  Kgr  Knz(s2  +  2£nzwnzs  +  “nz^^ 

-  =  -  X  -  x  -  (11) 

Fa  s2  +  2rfawfss  +  wfa2  Tas  +1  s2  +  2tapWaps  +  u>ap2 

The  steady  state  value  of  nz/Fa  for  a  step  input  was  used  to 
determine  the  required  gain.  Using  eq.  11,  the  final  value 
theorem  for  Laplace  transforms,  and  solving  for  the  gain  gives  the 
fol lowing: 


“sp2( nz/FS^33 

Kfs^nzwnz2 


(12) 


Eq.  12  was  used  to  calculate  the  gain  needed  for  a  stick 
force  per  g  of  6  lbs/g.  This  resulted  in  a  value  of  Kgr  =  0.2710 
rad/in  for  flight  condition  1  and  a  value  of  Kgr  =  0.0498  rad/in 
for  flight  condition  2.  (See 

a  description  of  flight  conditions  1  and  2).  To  check  that  these 
values  of  gain  were  correct,  they  were  used  in  eq.  11  along  with 
the  other  appropriate  parameters  for  flight  conditions  1  and  2. 

The  time  response  was  then  calculated  for  a  step  input  of  Fa  =  6.0 
lbs.  Figure  5  shows  that  the  step  input  results  in  the  desired 
steady  state  value  of  1.0  g  in  each  case. 
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Longitudinal  Airplane  Dynamics 

The  mathematical  models  used  to  simulate  longitudinal  aircraft 
response  to  elevator  inputs  are  developed  by  McRuer  et  al.  (5). 
The  derivation  begins  with  the  equations  of  motion  describing  how 
the  aircraft  deviates  from  trimmed,  steady-state  cruising  flight 
when  affected  by  small  aerodynamic  perturbations  or  limited 
control  inputs  and  results  in  the  following  fourth-order  transfer 
functions  for  pitch  rate  and  normal  acceleration: 


A$s(s  +  1/T@^)(s  +  I/T02) 

(s2  +  2?pWps  +  u>p2)(s2  +  2?gpWaps  +  w3p2) 


(13) 


Anzs(s  +  ^^TnZ2)(s2  +  2^nzwnzs  +  wnz^^ 
(s2  +  2?pWpS  +  Wp2)(s2  +  2rspu)SpS  +  wSp2) 


(14) 


Equations  13  and  14  are  the  universally  accepted  3  degree  of 
freedom  transfer  functions  used  for  describing  longitudinal 
aircraft  response  to  elevator  control  inputs.  For  most  flight 
conditions,  the  phugoid  and  short  period  natural  frequencies  (wp 
and  uiSp)  are  sufficiently  separated  that  their  dynamic  modes  are 
relatively  uncoupled  (table  I).  Because  of  this,  equations  13  and 
14  can  be  simplified  into  the  short  period  and  phugoid 
approximations  which  are  valid  over  their  respective  frequency 
ranges.  Generally,  the  phugoid  response  of  an  aircraft  is  very 
slow  compared  to  the  short  period  response  (i.e.  <<  wsp).  In 

an  approximate  sense,  the  phugoid  mode  describes  the  long-term 
translatory  motions  of  the  aircraft’s  center  of  mass,  whereas  the 
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short  period  mode  describes  rotations  about  the  center  of  mass. 


The  short  period  approximation  involves  setting  the  variation 

in  forward  velocity  u  equal  to  zero  with  the  assumption  that  it  is 

relatively  insignificant  in  the  short  period  mode  and  that  the  two 

degrees  of  freedom  a  and  6  are  dominant.  McRuer,  et  al.  state: 

"...for  frequencies  above  the  phugoid  the  two  degree  of 
freedom,  short  period  transfer  functions  are  very  good 
approximations  in  both  amplitude  and  phase.  Furthermore,  the 
two-  and  three  degree  of  freedom  time  responses  in  w  and  6 
are  in  excellent  agreement  for  times  shorter  than  about  10 
sec.”  (5) 


Table  I 


Values 

of  Up  and 

usp  at  Various 

Flight  Conditions* 

Mach  No 

Alt  (ft) 

Wt  (lbs) 

Up  ( rp3 ) 

Ugp  (rps 

0.20 

0 

22.058 

0. 1520 

1.558 

0.85 

0 

17,578 

0.0696 

7.350 

0.  40 

15,000 

17,578 

0.0980 

2.450 

0.60 

15,000 

17,578 

0.0635 

3.690 

0.50 

35,000 

17,578 

0.0861 

1.941 

0.70 

35,000 

17,578 

0.0752 

2.770 

*Data  taken  from  McRuer  (5)  for  the  A-4D  aircraft 
with  c.g.  location  at  0.25  X  c. 


Thus,  the  short  period  characteristics  are  dominant  during  the 
time  immediately  following  an  initial  perturbation  (1-10  seconds) 
after  which  the  phugoid  characteristics  begin  to  take  over.  For 
flying  qualities  specification,  the  primary  interest  is  in  the 
aircraft's  short  period  response  characteristics.  This  is 
evidenced  by  the  extensive  requirements  in  MIL-STD-1797  for  short¬ 
term  pitch  response  (para  4.2. 1.2  -  98  pages)  compared  to  the 
relatively  sparse  requirements  for  long-term  pitch  response  (para 
4. 2. 1.1  -  4  pages).  Also,  MIL-STD-1797  requires  that  equivalent 
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systems  be  matched  over  a  frequency  range  of  0. 1  to  10  rad/sec 

which,  for  typical  flight  conditions,  corresponds  to  the  frequency 

range  in  which  the  short  period  mode  is  dominant.  For  these 

reasons,  equivalent  systems  techniques  will  often  use  the 

appropriate  short  period  transfer  functions  as  the  low-order 

equivalent  systems.  The  low-order  equivalent  systems  used  by 

program  LONFIT  for  example  (chapter  III)  are  derived  from  the 

short  period  approximations  for  0/8e  and  nz/Se.  The  primary 

advantage  of  this  approach  is  that  it  provides  less  complicated 

low-order  systems  to  work  with  and,  because  wel 1 -accepted  short 

period  forms  are  used,  the  transfer  functions  give  meaningful 

short  period  parameters  when  matched  to  the  high-order  systems. 

MIL-STD-1797  suggests  using  the  3  degree  of  freedom  transfer 

functions  (eqs.  13  and  14)  in  equivalent  systems  work  but  states 

that  in  most  cases  the  phugoid  and  short  period  modes  are 

sufficiently  separated  that  "within  appropriate  respective 

frequency  ranges"  the  2  degree  of  freedom  approximations  can  be 

used.  The  requirements  state  that 

"While  no  specific  guidance  on  the  lower  frequency  bound  of 
the  matching  region  is  offered,  the  phugoid  and  short  period 
are  generally  separated  by  at  least  a  factor  of  10,  which 
should  be  adequate  to  consider  them  separately.  The 
assumption  of  widely  separated  phugoid  and  short  period  modes 
breaks  down  at  low  values  of  static  stability  (i.e.,  Mc  =  0) 
such  as  for  conventional  aircraft  with  extreme  aft  center  of 
gravity  locations  and  on  most  STOL  conf igurations. ”  (3) 

The  flight  conditions  and  associated  aircraft  parameters 

chosen  to  represent  the  basic  aircraft  dynamics  during  the 

equivalent  systems  matching  are  shown  in  table  II.  These  data  are 

taken  from  McRuer,  et  al.  (5)  for  the  A-4D  aircraft. 
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Table  II 


Flight 

Conditions  and 

Associated  Aircr 

aft  Parameters 

Fit  Cond 

Alt  (ft) 

Mach  No. 

vTo 

(fps)  q 

(psf  ) 

1 

35,000 

0.70 

681 

171 

2 

0 

0.85 

960 

945 

Flight 

Condition  1 

Flight 

Condition 

2 

^sp 

0.2250 

^sp 

=  0.4360 

”sp 

= 

2.7700 

"sp 

=  7.3500 

rp 

= 

0. 1177 

rp 

=  0.2260 

W— 

p 

= 

0.0752 

wp 

=  0.0696 

Ae 

= 

-11.33 

Ae 

=  -64.00 

1/T01 

= 

0.0091 

1/Tei 

=  0.0287 

1/T02 

2 

0.4280 

1/T82 

=  2.0800 

Anz 

2 

2. 1353 

Anz 

=  12.0690 

l'Tnal 

= 

0.0000 

1/Tnzl 

=  0.0000 

1/Tnz2 

= 

-0.0023 

1/Tnz2 

=  0.0268 

^nz 

= 

0.0260 

^nz 

=  0.0461 

wnz 

6.8800 

“nz 

=  18.0500 

These  flight  conditions  were  chosen  primarily  because  of  the 
values  of  uiSp  as  will  be  explained  in  the  following  section.  Note 
that  in  both  flight  conditions  the  natural  frequencies  wp  and  uSp 
are  sufficiently  separated  (i.e.  by  a  factor  of  10  or  more)  so 
that  the  short  period  and  phugoid  modes  can  be  considered 
separately.  Also,  since  wp  is  less  than  0.1,  the  short  period 
mode  will  be  dominant  in  the  frequency  range  of  interest  (0.1  to 
10  rad/sec).  The  short  period  approximation  should  therefore  be 
appropriate  to  use  during  the  equivalent  systems  matching  process. 
This  will  be  true  as  long  as  the  control  system  roots  do  not  cause 
unfavorable  coupling  effects  with  wSp  thereby  degrading  the 
equivalent  systems  matching  ability. 
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HOS  Configurations 


The  effect  of  control  system  dynamics  is  to  raise  the  order 
of  the  airplane  response  to  pilot  inputs.  The  airplane  responses 
(eqs.  13  and  14)  are  fourth  order  responses  for  elevator  control 
inputs.  The  additional  dynamics  of  the  elevator  response  to 
control  stick  displacements  (servo-actuator)  and  the  control  stick 
response  to  control  stick  forces  (feel  system)  will  increase  the 
order  of  the  overall  airplane  response  to  control  stick  forces. 

In  addition,  the  feel  system  and  actuator  roots  may  be  near  the 
airplane  roots.  Increased  order  and  closeness  of  roots  can  alter 
the  responses  of  the  airplane  to  stick  force  inputs,  affect  the 
pilot’s  closed  loop  control,  and  change  the  airplane’s  flying 
qualities. 

Studies  of  many  different  types  of  flight  control  systems 
have  shown  that  most  systems  have  certain  similar  effects  on  the 
airplane's  overall  maneuvering  characteristics.  In  general,  it 
can  be  said  that  a  flight  control  system  will  do  one  of  two  things 
to  the  overall  aircraft  response:  1)  introduce  a  ’’pure”  time 
delay  without  changing  the  basic  shape  of  the  time  response,  or 
2)  introduce  time  delay  and  change  the  shape  of  the  time  response. 
The  first  case  (pure  time  delay)  is  characterized  by  control 
system  roots  with  relatively  high  natural  frequencies  compared  to 
the  aircraft  roots.  This  type  of  control  system  will  have  a 
"fast”  or  at  least  ’’moderately  fast”  response  compared  to  the 
basic  aircraft  response.  The  second  case  is  characterized  by 
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insufficient  separation  of  the  flight  control  system  and  airplane 
roots  which  causes  a  change  in  the  basic  aircraft  response 
characteristics.  This  type  of  control  system  will  have  a  "slower" 
response  because  of  its  lower  natural  frequency.  In  either  case, 
the  time  delay,  if  it  becomes  excessive,  is  generally  cited  by 
pilots  as  the  most  objectionable  characteristic  of  flight  control 
systems. 

In  theory,  the  equivalent  system  matching  process  should  be 
relatively  straightforward  for  a  high-order  system  with  pure  time 
delay.  In  this  type  of  system,  the  natural  frequencies  of  the 
control  system  roots  are  high  compared  to  the  aircraft  roots  and 
their  effects  can  be  neglected  in  the  frequency  band  of  interest 
(0.1  to  10  rad/sec).  The  only  source  of  mismatch,  therefore,  is 
from  the  high-frequency  phase  contributions  of  the  control  system 
roots;  and  this  can  be  approximated  by  a  time  delay.  In  the 
DiFranco  (10)  high-order  system  evaluations,  for  example,  it  was 
found  that  the  frequency  and  damping  ratio  of  the  equivalent 
systems  were  very  close  to  the  airplane’s  short  period 
characteristics  when  the  control  system  roots  were  high  relative 
to  ui3p.  Thus,  the  high-order  systems  could  be  closely  matched  to 
equivalent  second-order  systems  (i.e.  the  short  period 
approximation)  using  a  time  delay  term. 

In  the  second  type  of  high-order  system,  the  time  response 
has  been  altered  because  the  control  system  roots  are  close  to  the 
aircraft  roots.  This  type  of  system  may  not  match  well  to  a 
lower-order  equivalent  system  because  of  the  combined  effects  of 
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the  aircraft  and  control  system  roots.  One  of  the  conclusions 
reported  in  the  DiFranco  evaluations  for  example  states  the 
fol lowing: 

"With  higher-order  control  system  dynamics,  the  airplane 
short  period  response  to  step  stick  force  inputs  can  be 
reasonably  well  represented  by  a  time  delay  and  an  equivalent 
second-order  response. .. This  simplified  representation  is 
poorest  when  the  lowest  frequency  of  the  control  system  is 
near  the  airplane  short  period  frequency.”  (10) 

This  problem  was  also  reported  in  the  Neal-Smith  evaluations  which 

tested  various  configurations  with  the  control  system  roots  close 

to  the  aircraft  roots.  Comparing  the  pure  time  delay  type  of 

high-order  system  with  the  systems  used  in  their  study,  the  report 

states  that:  “Since  the  control  system  dynamics  studied  in  the 

present  experiment  significantly  alter  the  shape  of  the  airplane 

response  to  pilot  inputs,  determining  an  equivalent  system  is  a 

more  difficult  proposition. ”  The  report  goes  on  to  state  that 

because  of  the  practical  difficulties  associated  with  the 

equivalent  systems  approach,  a  pi lot-in-the-loop  analysis  was  used 

instead  (11). 

Based  on  the  results  of  these  investigations,  it  was 
concluded  that  to  adequately  evaluate  equivalent  systems  matching 
techniques  it  would  be  necessary  to  use  high-order  systems 
representative  of  both  types  discussed  above.  The  aircraft  and 
feel  system  parameters  were  therefore  chosen  to  simulate  high- 
order  systems  with  "fast”  response  times  (resulting  in  pure  time 
delay)  and  high-order  systems  with  "slow”  response  times 
(resulting  in  time  delay  plus  change  in  time  response).  The 
flight  conditions  and  parameters  selected  for  these  systems  are 
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summarized  in  table  III. 


Table  III 


High-Order 

System  Conf igurations 

HOS  Des¬ 
ignation 

Flight 

Condition 

HOS 

Type 

Kgr 

( rad/in ) 

wfs 

( rad/s ) 

sp 

( rad/s ) 

HOSi 

(§/Fa) 

1 

Uncoupled 

0.2710 

18.  5 

2.77 

hos2 

( nz/F a ) 

1 

Uncoupled 

0. 2710 

18.5 

2.77 

HOSg 

(§/Fa) 

2 

Coupled 

0.0498 

6.0 

7.35 

hos4 

( nz/Fa ) 

2 

Coupled 

0.0498 

6.0 

7.35 

(Kfa  =  0.046  in/lb  and  1/Ta  =  13  rad/s  for  all  HOS  configurations) 


Note  that  for  clarity  throughout  the  equivalent  systems 
matching  process,  the  high  order  systems  with  pure  time  delay  will 
be  designated  as  HOS^  (0/ Fa )  and  HOS2  (nz/Fa)  and  the  higher-order 
systems  with  coupling  between  the  control  system  and  aircraft 
dynamics  will  be  designated  as  HOSg  (&/Fa)  and  H0S4  (nz/Fa)  as 
listed  in  table  III.  Although  ui£a  will  be  varied  throughout  the 
range  of  6  to  31  rad/sec  resulting  in  a  variety  of  other  HOS 
configurations  (as  explained  in  chapter  III)  the  HOS^  -  H0S4 
configurations  will  be  the  primary  (baseline)  systems  used 
throughout  the  frequency  domain  and  time  domain  comparisons. 

The  high-order  configurations  were  chosen  because  of  the 
relative  values  of  1/Ta,  ui£a,  and  wap.  The  dynamics  of  HOS^  and 
HOSg*  for  example,  are  the  result  of  the  relatively  low  short 
period  frequency  of  the  aircraft  (wap  =  2.77)  together  with  the 
relatively  high  natural  frequency  of  the  feel  system  ( u>£a  =  18.5) 
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and  corner  frequency  of  the  servo  (1/Ta  *  13).  This  combination 
provides  good  separation  between  the  control  system  and  aircraft 
roots  and  results  in  a  relatively  uncoupled,  pure  time  delay, 
high-order  system  response.  The  magnitude  plots  for  HOS^  and  H0S2 
(figures  6a  and  7a)  show  that  the  dynamics  of  the  feel  system  and 
servo  are  essentially  removed  from  the  frequency  range  of  interest 
(0.1  to  10  rad/sec).  The  only  difference  between  the  basic 
aircraft  and  the  HOS,  in  this  frequency  range,  is  an  increase  in 
gain  for  the  HOS  equal  to  ( Kf3Kgru)£a2 ) /Ta.  The  HOSj^  and  H0S2 
phase  plots  (figures  6b  and  7b)  show  the  phase  lags  of  the  control 
system  dynamics  which  correspond  to  the  time  delay.  Figure  8 
compares  the  time  responses  of  HOS^  and  H0S2  with  the 
corresponding  basic  aircraft  responses.  The  time  responses  are 
shown  for  a  step  input  to  the  HOS  of  Fa  =  5  lbs  and  an  equivalent 
step  input  to  the  aircraft  of  8e  =  3.57  deg.  Note  that  the  time 
responses  of  the  high-order  systems  show  a  time  delay  but  have 
essentially  the  same  shape  as  the  basic  aircraft  responses. 

The  HOSg  and  HOS4  dynamics  are  influenced  by  the  relatively 
high  short  period  frequency  of  the  aircraft  (usp  =  7.35)  combined 
with  the  relatively  low  natural  frequency  of  the  feel  system  (u)£S 
*  6.0).  This  combination  gives  very  little  separation  between  the 
aircraft  and  control  system  roots  and  results  in  a  coupled  high- 
order  system  response.  Compared  with  the  HOS^  and  H0S2  Bode 
plots,  the  HOS3  and  HOS4  plots  (figures  9  and  10)  show  the 
presence  of  W£a  within  the  0. 1  to  10  rad/sec  frequency  band  and 
indicate  increased  phase  lags  corresponding  to  a  greater  time 
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delay.  The  greater  time  delay  is  expected  because  of  the 
relatively  slow  fee  1  system  response  at  uifa  =  6.0  rad/sec.  The 
presence  of  u»£S  within  the  0. 1  to  10  rad/sec  frequency  band 
results  in  coupling  effects  between  the  control  system  and 
aircraft  dynamics.  Figure  11  shows  the  time  responses  of  HOSg  and 
HOS4  compared  to  the  corresponding  basic  aircraft  responses.  The 
time  responses  are  for  a  step  input  to  the  HOS  of  Fa  =  5  lbs  and 
an  equivalent  step  input  to  the  aircraft  of  8e  =  0.66  deg.  Note 
that  the  HOSg  and  HOS4  time  responses  show  an  increase  in  time 
delay  (compared  to  HOS^  and  HOSg)  as  well  as  a  change  in  the  basic 
shape  of  the  time  response. 
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PHASE  (DEG)  MAG  (DB) 


.01  .1  1  10  100 

FREQ  (RPS) 


(a)  Magnitude 
HOSi  - 


FREQ  (RPS) 


(b)  Phase 


Figure  6.  Frequency  Response  Comparison,  HOS^  and  Basic  Aircraft 
(Flight  Condition  1  and  W£S  -  IS. 5  rad/sec) 
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(b)  Phase 


Frequeny  Response  Comparison,  H0S2  and  Basic  Aircraft 
(Flight  Condition  1  and  wfa  *  18.5  rad/sec) 


NZ  (G) 


PHASE  (DEG)  MAG  (DB) 


HOSg 


FREQ  (RPS) 


(a)  Magnitude 


.01  .1  1  10  100 

FREQ  (RPS) 

( b )  Phase 


Figure  9.  Frequeny  Response  Comparison,  H0S3  and  Basic  Aircraft 
(Flight  Condition  2  and  u»£s  =  6.0  rad/sec) 
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PHASE  (DEG)  MAG  (DB) 


hos4 


FREO  (RPS) 

(a)  Magnitude 


H0S4 


FREQ  (RPS) 

(b)  Phase 


Figure  10.  Frequeny  Response  Comparison,  H0S4  and  Basic  Aircraf 
(Flight  Condition  2  and  u i£s  =  6.0  rad/sec) 
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(f))  7N 


TIME  (SEC) 

(a)  HOS3  versus  Basic  Aircraft 
(Flight  Condition  2  and  iofs  =  6.0  rad/s) 


TIME  (SEC) 

(b)  HOS^  versus  Basic  Aircraft 
(Flight  Condition  2  and  u>£s  =  6.0  rad/s) 


Time  Response  Comparisons  (step  input  F3  *  5.0  lbs  to 
HOS  and  equivalent  step  input  8P  =  0.66  deg  to  A/C) 


Figure  11. 


III.  Frequency  Response  Computer  Matching 


Computer  Programs 

The  frequency  domain  matching  was  accomplished  using  the 
computer  programs  LONFIT  and  MACFIT  (18,19).  These  programs, 
developed  by  McDonnell  Douglas,  are  applicable  to  MIL-F-8785C 
para  3.2.2. 1.1  (short  period  frequency  and  acceleration 
sensitivity)  and  para  3. 2. 2. 1.2  (short  period  damping) 
requirements  (2). 

LONFIT  matches  high-order  pitch  rate  and  normal  acceleration 
transfer  functions  with  the  following  equivalent  low-order 
systems: 


0  Ke(s  +  La)e  ST 

Fs  +  2£SpW3pS  +  wsp2 


K  e“ST 
^nz® 


s2  +  2tapu)sps  +  w3p2 


(15) 


( 16) 


where  the  parameters  Kq,  Knz,  La,  £sp,  u3p,  and  t  are  determined 
by  simultaneously  matching  the  frequency  response  of  the  above 
transfer  functions  to  those  of  the  airframe  plus  flight  control 
system  (high-order  system). 

LONFIT  is  configured  for  longitudinal  systems  only  and  the 
equivalent  low-order  forms  are  restricted  to  the  specific  transfer 
functions  of  equations  15  and  16.  MACFIT,  on  the  other  hand,  is  a 
general  purpose  matching  program  which  can  accept  arbitrary  HOS 
and  LOSS  forms.  Although  LONFIT  (due  to  the  pre-chosen  LOES 
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forms)  is  somewhat  easier  to  use,  MACFIT  has  the  advantage  of 


versatility  in  cases  where  the  structure  of  LONFIT  is  inadequate. 

LONFIT  and  MACFIT  both  use  a  Rosenbrock  multivariable  search 
algorithm  to  match  Bode  plots  of  the  low-order  and  high-order 
systems  by  minimizing  the  weighted  sum  of  the  squares  of  the 
differences  in  magnitude  and  phase  angle  between  the  systems  at  n 
discrete  frequencies.  The  following  cost  function  is  minimized: 


COStf 


u> 


n 


20/ 


„  2 


u>. 


2 

WT1( Again) 


2 

WT2( Aphase ) 


I 


(17) 


where:  cost_£  is  a  frequency  domain  mismatch  function 

WT1  =  1.0 
WT2  =  0.01745 

Again  =  gainHQS  “  gainL0S  (in  decibels) 

Aphase  =  phase^os  “  PhaseL0S  < in  degrees) 
and  n  is  the  number  of  frequencies 


To  begin  the  matching  process,  the  coefficients  (or  roots)  of 
the  HOS  numerator  and  denominator  are  required  as  input.  The  user 
then  has  the  flexibility  to  specify  initial  values  for  the  LOS 
parameters,  which  parameters  the  program  may  vary,  the  number  of 
frequencies  (n),  the  frequency  range,  the  number  of  iterations, 
whether  or  not  to  include  time  delay,  and  whether  or  not  to  permit 
the  resulting  transfer  functions  to  be  unstable. 

The  subscript  f  in  equation  17  is  used  to  denote  that  "costf" 
is  a  measure  of  the  mismatch  in  the  frequency  domain.  Because  it 
will  be  necessary  to  compare  time  domain  and  frequency  domain 
matching  in  this  study,  it  will  be  helpful  to  define  a  second  cost 
function  ”costt”  which  will  be  used  to  measure  the  mismatch  of  the 
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HOS  and  LOES  time  responses.  Cost^.  will  be  defined  as  follows: 


costt 


*k 

1/k  2 


2 

( Athetadot ) 


(18) 


where:  cost^.  is  a  time  domain  mismatch  function 

Athetadot  =  8  ^OS  -  0  j^gg  (in  degrees) 
and  k  is  the  number  of  discrete  data  points 

The  time  response  (9)  in  equation  18  is  calculated  using  a 
step  input  of  5  lbs.  Although  cost^.  has  been  defined  using  the 
mismatch  in  0  a  similar  definition  could  be  used  for  nz  mismatch. 
In  this  thesis,  however,  comparisons  between  time  domain  and 
frequency  domain  LOES  matching  techniques  will  only  be  made  for 
the  6  matching  results.  The  time  response  mismatch  has  been 
defined  such  that  the  magnitude  of  costt  will  be  much  less  (i.e. 
on  the  order  of  100  less)  than  the  magnitude  of  cost^  for  a  given 
H0S/L0ES  mismatch.  This  has  been  done  to  avoid  the  temptation  of 
making  direct  comparisons  between  the  magnitudes  of  costf  and 
costj..  For  example,  a  time  domain  mismatch  of  cost^.  <  200  can  not 
be  used  as  a  guideline  for  a  "good”  match  as  is  generally  the  case 
for  cost£. 


Matching  Strategy 

The  HOS^  and  HOSg  transfer  functions  derived  in  chapter  II 
were  matched  to  the  corresponding  LOS^  and  LOSg  forms  shown  in 
eqs.  15  and  16  using  the  program  LONFIT.  The  following  matching 
strategy  was  taken  primarily  from  the  LONFIT  users  guide  (18) 
although  it  has  been  modified  somewhat  for  this  particular 
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example. 


(1)  A  frequency  range  of  0.1  to  10  rad/sec  was  used  in 
accordance  with  MIL-F-8785C  equivalent  systems  criteria  and  is 
intended  to  cover  the  frequency  range  of  major  interest  to  the 
pilot.  The  number  of  discrete  frequencies  (n)  was  set  to  20  as 
suggested  in  the  LONFIT  users  guide. 

(2)  La,  £sp,  and  wsp  were  initially  set  equal  to  their 
high-order  system  values  (taken  from  the  short-period  roots  of  the 
basic  aircraft  transfer  function).  Initial  values  of  gain  and 
time  delay  were  arbitrarily  set  to  1.0  and  0.1  respectively. 

(3)  To  investigate  the  effects  of  La  and  time  delay  che 
H0S  frequency  response  was  matched  by  varying  gain,  £sp,  and  tuSp 
with  the  following  combinations  of  La  and  time  delay: 


H0St  Matching  H0Sg  Matching 

La  Time  Delay  Time  Delay  LOS  Form 


LOSia 

Fixed 

No 

L0S2A 

No 

LONFIT 

LOSib 

Fixed 

Yes 

J-0S2B 

Yes 

LONFIT 

LOSic 

Free 

No 

LOS2C 

No 

MACFIT 

LOS  id 

Free 

Yes 

L0S2D 

Yes 

MACFIT 

(4)  The  LOS2C  ®nd  LOS2D  matches  were  carried  out  using 
program  MACFIT  and  a  second  order  numerator,  second  order 
denominator  form  for  the  equivalent  nz/Fg  low-order  system.  This 
is  because  the  LONFIT  form  of  eq.  16,  which  is  an  abbreviated  form 
of  the  short-period  approximation,  did  not  match  well  with  the 
high-order  nz/Fs  transfer  function.  The  complete  short-period 
approximation  was  therefore  used  and  produced  good  results.  (This 
is  explained  further  in  the  section  on  HOS2  matching). 
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HOS-)  Matching 


The  8/Fs  high-orJer  system  (HOS^)  was  matched  to  the  LOSj^  - 
LOSid  equivalent  systems  using  the  first  three  steps  outlined 
above.  The  resulting  LOS  parameters  and  cost£  functions  are  shown 
Table  IV.  The  cost  function  trends  seen  in  Table  IV  demonstrate 
that  increasing  the  the  number  of  free  variables  (i.e.  increasing 
the  dimension  space  of  the  low-order  system)  generally  results  in 
a  better  match  to  the  high-order  system.  Based  on  these  results, 
the  logical  approach  to  achieving  the  best  match  would  apparently 
be  to  free  all  the  available  parameters.  The  problem  with  this 
approach,  however,  is  that  questions  can  arise  regarding  the 
meaning  of  some  free  variables.  Disagreement  continues  in  the 
flying  qualities  community,  for  example,  as  to  whether  the 
parameter  La  should  be  fixed  at  the  basic  aircraft  value  or  be 
determined  by  the  match  (3).  The  rationale  behind  keeping  it 
fixed  is  that  La  is  determined  by  wing  size  and  planform  which  are 
dimensions  that  can  not  be  varied  significantly  for  a  particular 
aircraft.  It  is  argued,  therefore,  that  La  should  be  kept 
constant  at  the  basic  aircraft  value.  Furthermore,  in  studies 
where  it  has  been  allowed  to  vary,  the  results  have  produced 
unreasonable  fluctuations  in  La  from  one  configuration  to  another. 
The  point  is  that  even  though  the  lowest  cost  function  is  achieved 
by  freeing  all  the  parameters,  this  may  not  be  the  best  overall 
approach  if  ambiguity  arises  regarding  the  meaning  of  some  terms. 
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Table  IV 


§/Fs  Matching  Results 

Flight  Condition  1  and  u)£S  =  18.6  rad/s  (HOS^) 


System 

La 

Ke(A0) 

^sp 

usp 

La 

T 

COSt£ 

A-4D* 

- 

- 

(-11.33) 

0.  225 

2.770 

0.428 

- 

— 

los1a 

Fixed 

No 

-0. 117 

0.  180 

2.435 

0.428 

0 

441 

LOSib 

Free 

No 

-0.099 

0.  128 

2.617 

0.786 

0 

374 

LOSic 

Fixed 

Yes 

-0. 133 

0.238 

2.601 

0.428 

0.  164 

82 

losiD 

Free 

Yes 

-0. 120 

0.  193 

2.686 

0.595 

0.  156 

60 

*  Basic  A-4D  aircraft  parameters  (without  fee 1  system  and  servo 
dynamics)  taken  from  McRuer  (5). 

For  some  equivalent  system  configurations  there  may  be  no 
choice  but  to  vary  all  the  parameters.  To  achieve  an  acceptable 
match  for  the  longitudinal  high-order  systems  in  the  Neal-Smith 
evaluations,  for  example,  the  report  stated  that  all  parameters, 
including  La,  had  to  be  varied  (11).  In  the  present  example, 
however,  an  excellent  match  was  obtained  using  the  LOS-^c 
configuration  with  La  fixed. 

Magnitude,  phase  angle,  and  time  response  comparisons  between 
HOS^  and  the  LOS^q  and  LOS^q  configurations  are  shown  in  Figures 
12-14.  To  obtain  the  phase  angle  and  time  response  plots  it  was 
necessary  to  represent  the  time  delay  term,  e-ST,  by  poles  and 
zeros  directly.  This  was  done  using  a  first  order  Fade' 
approximation  as  shown  below: 


1  -  st/2 

e-ST  =  -  (19) 

1  +  st/2 
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The  effect  of  time  delay  can  be  seen  in  Figure  15  which 
compares  the  phase  mismatch  between  HOS^  and  LOS^c  both  with  and 
without  the  Pade'  approximation.  Figure  15b  which  does  have  the 
Pade'  approximation  for  e~ST  illustrates  that  equivalent  time 
delay  approximates  the  phase  lags  of  the  high-frequency  terms 
quite  well. 

Figure  16  shows  the  corresponding  time  domain  comparisons 
between  HOS^  and  LOS^q  for  a  step  input  of  Fs  =  5  lbs.  Figures 
15a  and  16a  do  not  have  the  time  delay  term  and  illustrate  that  a 
poor  match  in  the  high-frequency  region  corresponds  to  a  poor 
match  in  the  transient  time  response  region.  On  the  other  hand. 
Figures  15b  and  16b,  which  both  include  time  delay,  show  that  a 
good  match  in  the  high-frequency  region  corresponds  to  a  good 
match  in  the  transient  time  response  region.  These  results 
demonstrate  the  initial  value  theorem  for  Laplace  transforms. 

This  theorem  states  that  the  behavior  of  f(t)  in  the  neighborhood 
of  t  =  0  is  related  to  the  behavior  of  sF(s)  in  the  neighborhood 
of  a  **  •».  This  behavior  is  clearly  evident  in  Figures  15  and  16. 

HOSp  Matching 

The  nz/Fs  high-order  system  for  flight  condition  1  and  u»£S  = 
18.5  <H0S2>  was  matched  to  LOS2A  anc*  LOS2B  using  program  LONFIT 
and  the  first  three  steps  outlined  in  the  matching  strategy.  The 
resulting  LOS  parameters  and  cost  functions  are  shown  in  Table  V. 
The  time  delay  was  again  applied  using  a  first  order  Pade' 
approximation  to  obtain  magnitude  and  phase  comparisons  for  these 
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Figure  12. 

Frequency  Response  Comparison,  HOSa  vs  LOS 
(Costf  =  82) 
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Figure  13.  Frequency  Response  Comparison,  HOS-i  vs  LOS 
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(a)  LOS^c  without  Pade  Approximation 
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(b)  L0S1C  with  Pade  Approximation 


Figure  15.  LOS^c  Phase  Angle  Comparisons  -  Effects  of  Time  Delay 
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THETADOT  (RPS) 


configurations.  These  plots  (Figures  17  and  18)  and  the 
corresponding  cost  functions  in  Table  V  show  a  relatively  poor 
match  for  the  LOS2A  and  LOS2B  systems. 

_ Table  V _ 

nz/Fs  Matching  Results 
Flight  Condition  1  and  ui£S  =  18.5  (HOS2) 


System 

LOS 

Form 

Knz 
<  Anz ) 

^nz 

wnz 

^sp 

wsp 

T 

cost 

A-4D* 

- 

(2. 135) 

0.026 

6.  880 

0.225 

2.  770 

- 

- 

L0S2A 

LONFIT 

0.749 

- 

- 

0.  254 

2.  166 

0 

679 

l0s2B 

LONFIT 

0.713 

- 

- 

0.229 

2.  110 

-0.074 

604 

los2C 

MACFIT 

0.  176 

0.104 

7.790 

0.  193 

2.386 

0 

394 

los2D 

MACFIT 

0.  174 

0.022 

6.999 

0.238 

2.  601 

0.  161 

87 

*  Basic  A— 4D  aircraft  parameters  (without  feel  system  and  servo 
dynamics)  taken  from  McRuer  (S). 


The  results  of  the  LOS2A  and  LOS2B  matching  show  that  a 
better  match  is  obtained  when  time  delay  is  used  but  only  if  it  is 
allowed  to  be  negative.  The  question  of  whether  or  not  to  allow 
negative  time  delay  has  come  up  in  piloted  evaluations. 

Equivalent  systems  data  from  the  Neal-Smith  program  (11)  show  that 
three  of  the  five  lead/ lag  systems  that  incorrectly  predicted 
flying  qualities  (i.e.  predicted  level  1  but  were  rated  level  2) 
could  have  produced  a  better  match  if  t  was  allowed  to  be  less 
than  zero.  The  Neal-Smith  program,  however,  did  not  allow 
negative  time  delay  and  t  was  set  to  zero  for  these 
configurations.  In  an  analysis  of  these  data,  Mitchell  and  Hoh 
(20)  argue  that  if  t  was  allowed  to  be  negative  the  correlation 
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between  predicted  flying  qualities  and  pilot  ratings  could  have 
been  improved  for  these  configurations. 

Although  negative  time  delay  did  improve  the  match  of  the 
LOS2B  configuration,  the  improvement  still  did  not  result  in  a 
"good”  match  to  the  high-order  system.  The  mismatch  seen  in 
Figures  17  and  18  suggests  that  a  more  fundamental  system  change 
(other  than  time  delay)  is  needed.  The  problem  is  in  the  basic 
form  of  the  low-order  equivalent  system  used  by  LONFIT.  This 
form,  with  a  simple  gain  in  the  numerator,  is  an  approximation  to 
the  short  period  approximation  for  nz/Be.  Although  this 
approximation  is  commonly  used  in  equivalent  systems  studies, 
McRuer,  et  al.  (5)  show  that  this  approximation  is  valid  only  at 
frequencies  below  the  numerator  break  frequencies  as  shown  below: 


l/g(MSeZa  -  MaZSe)  1  1 

-  ,  s  <  -  ,  -  ,  or  wnz  (20) 

S2  +  STgpWgpS  +  Ulgp2  T^  T^ 


Thus,  the  LOS2A  and  LOS2B  systems  do  not  match  well  to  the  high- 
order  system  because  of  the  quadratic  numerator  term  at  wnz  -  6.88 
rad/sec.  The  frequency  range  of  interest  (0.1  to  10  rad/sec)  lies 
above  this  numerator  break  frequency  and  eq.  20  is  therefore  not  a 
good  approximation  in  this  case.  For  frequencies  above  the  break 
frequencies,  McRuer,  et  al.  (5)  state  that  the  full  short  period 
approximation  (eq.  21)  should  be  used. 


Knz(s2  +  2r„^u inws  +  u> 


nz“n  z- 


2) 
'nz  ' 


8< 


i2  +  2£anw 


(21  ) 


spwsp 


s  +  w 


sp 
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Equation  21  was  used  as  the  low-order  equivalent  form  for  the 


LOS2C  and  LOS2D  systems  and  matched  to  HOS2  using  program  MACFIT. 
On  MACFIT,  this  simply  involved  specifying  a  second  order 
numerator  and  second  order  denominator  as  the  “guets”  for  the  low- 
order  system  form.  The  excellent  results  (Table  V  and  Figures  19- 
21)  show  that  the  full  short  period  approximation  is  indeed  an 
appropriate  low-order  system  to  use  for  this  situation. 

The  difficulties  encountered  during  the  Nz/Fs  matching 
demonstrate  the  ’’guess  work”  involved  in  certain  parts  of  the 
matching  process.  The  fundamental  problem,  for  example,  became 
one  of  finding  the  best  low-order  system  to  use.  The  form 
typically  used  in  equivalent  system  work  (eq.  16)  was  found  to  be 
inadequate  and  a  higher-order  system  ( eq .  21)  was  required.  Since 

the  higher-order  system  improved  the  match,  a  logical  question 
which  follows  is  whether  or  not  an  even  better  LOS  form  exists. 
Certainly,  a  better  match  is  possible  as  the  order  of  the 
equivalent  system  is  increased,  but  carrying  this  process  too  far 
may  defeat  the  purpose  of  the  equivalent  system  approach.  That 
is,  the  existing  data  base  for  flying  qualities  parameters  is  for 
classical  (unaugmented)  airplane  systems  and  a  certain  amount  of 
mismatch  must  be  tolerated  in  order  to  get  the  high-order  system 
into  one  of  these  low-order  forms.  Among  the  classical  forms, 
however,  there  may  be  several  options  to  choose  from  (as  we  have 
seen  in  this  example)  and  so  the  approach  is  not  always 
straightforward.  The  individual  performing  the  match  must 
consider  what  LOES  parameters  are  desired  and  decide  what  level  of 
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Figure  19.  Frequency  Response  Comparison,  H0S2  vs  L0S2C 
(Costf  =  394) 
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Figure  20.  Frequency  Response  Comparison,  HOS2  vs  LOS20 


mismatch  is  acceptable.  With  these  constraints  in  mind  the  LOS 
form  can  be  chosen  from  the  available  options. 

HOS<a  and  H0S4  Hatching 

The  0/Fs  and  nz/Fs  high-order  systems  for  flight  condition  2 
and  ui£S  =  6.0  (  HOSq  and  H0S4 )  were  matched,  using  the  first  three 
steps  outlined  in  the  matching  strategy,  to  the  following  low- 
order  systems: 


HOSq  Matching  H0S4  Matching 

La  Time  Delay  Time  Delay  Los  Form 


LOSqa 

Fixed 

No 

los4A 

No 

LONFIT 

LOSqb 

Free 

No 

los4B 

Yes 

LONFIT 

LOSqc 

Fixed 

Yes 

LOSqd 

Free 

Yes 

The  resulting  LOS  parameters  and  costf  functions  for  the  9/Fs 
(HOS3)  matching  are  shown  in  Table  VI.  Note  that  the  relatively 
slow  feel  system  dynamics  (ui£S  =  6.0)  of  these  HOS  configurations 

_ Table  VI _ 

™  0/Fs  Matching  Results 
Flight  Condition  2  and  u)£a  =  6.0  (HOSq) 


System 

La 

Time 

Delay 

Ke 

<A0) 

^sp 

wsp 

La 

T 

cost 

A— 4D* 

- 

- 

(-64.00) 

0.436 

7.350 

2.080 

- 

- 

LOSqa 

Fixed 

No 

-0.030 

0.449 

3.  194 

2.080 

- 

511 

LOSqb 

Free 

No 

-0.016 

0.251 

4.837 

9.616 

- 

274 

L0SqC 

Fixed 

Yes 

-0.059 

0.720 

4.524 

2.080 

0.220 

58 

L0SqD 

Free 

Yes 

-0.027 

0.  359 

5.887 

8.  197 

0.  168 

33 

*  Basic  A-4D  aircraft  parameters  (without  feel  system  and  servo 
dynamics)  taken  from  McRuer  (5). 
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resulted  in  a  significant  increase  in  time  delay.  Low-order 


equivalent  systems  without  a  time  delay  approximation  therefore 
did  not  match  well  to  these  high-order  systems. 

The  equivalent  system  parameters  and  cost^  functions  for  the 
nz/Fs  (H0S4)  matching  are  shown  in  Table  VII.  Since  wnz  =  18.05 
rad/sec  for  these  configurations  (which  is  well  above  the  0.1  to 
10  rad/sec  frequency  band),  eq.  20  shows  that  the  simple  LOES  form 
used  by  LONFIT  should  be  adequate  to  approximate  the  short  period 
dynamics.  The  excellent  results  obtained  using  LONFIT  show  that 
this  approximation  is  indeed  valid.  Thus,  there  is  no  need  to  use 
the  full  short  period  approximation  (and  program  MACFIT)  as  was 
the  case  for  flight  condition  1  where  u>nz  =  6.88  rad/sec. 
Magnitude,  phase,  and  time  response  plots  for  L0S3C  and  L0S4B  are 
shown  in  Figures  22-24. 


Table  VII 


n_/Fa  Matching  Results 

Flight 

Condition 

2  and  u>£S 

=  6.0 

(H0S4) 

System 

LOS 

^nz 

^  sp 

U)— — 

sp 

T 

COSt£ 

Form 

(  Anz ) 

A-4D* 

- 

( 12.059) 

0.436 

7.350 

- 

- 

los4A 

LONFIT 

1.679 

0.460 

3.066 

0 

457 

los4B 

LONFIT 

2.888 

0.694 

4.076 

0.  196 

80 

*  Basic  A-4D  aircraft  parameters  (without  feel  system  and  servo 
dynamics)  taken  from  McRuer  (5). 


A  cursory  comparison  between  the  "pure" 
systems  of  HOS^  and  HOSg  and  the  time  delay 


time  delay  high-order 
plus  change  in  time 
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(0)  ZN 


response  systems  of  HOSg  and  HOS^  does  not  show  any  significant 
differences  between  the  coupled  and  uncoupled  systems  as  was 
expected.  In  fact,  considering  only  the  costf  functions,  the 
coupled  systems  (HOS3  and  H0S4 )  actually  matched  better  than  the 
uncoupled  systems  (HOS^  and  HOS2)  since  the  former  had  slightly 
lower  cost£  functions  -  indicating  a  closer  frequency  domain 
match.  The  cost£  functions,  however,  were  all  below  100  for  the 
best  matches  which,  from  a  flying  qualities  perspective,  indicates 
an  excellent  match.  The  HOS^  -  HOS4  systems,  therefore,  have  been 
useful  in  providing  a  detailed  look  at  the  matching  process  but 
have  provided  only  a  limited  look  at  the  effects  of  coupling 
between  the  control  system  and  aircraft  dynamics.  For  this 
reason,  a  more  complete  range  of  control  system  dynamics  was 
investigated.  The  results  of  this  analysis  are  given  in  the 
following  section. 

Effects  of  Control  System  Dynamics 

To  investigate  the  effects  of  coupling  between  the  control 
system  and  aircraft  dynamics  on  equivalent  systems  matching 
techniques,  the  feel  system  natural  frequency  u)£S  was  varied 
throughout  the  range  of  6  to  31  rad/sec  for  each  flight  condition. 
This  was  done  using  the  pitch  rate  8/Fs  high-order  system  with  La 
both  fixed  and  free  which  resulted  in  a  total  of  12  HOS  and  24 
LOES  configurations.  The  tabulated  results  of  this  matching  are 
given  in  Appendix  A,  and  selected  plots  of  the  data  are  presented 
here  for  consideration.  The  results  of  this  matching  demonstrate 
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trends  associated  with  the  movement  of  the  control  system  roots 
and  illustrate  some  of  the  problems  in  equivalent  systems  matching 
when  there  is  coupling  between  the  control  system  and  aircraft 
dynamics. 

In  general,  the  results  of  this  matching  show  that  as  the 
feel  system’s  natural  frequency  is  moved  away  from  the  aircraft’s 
natural  frequency  a  better  equivalent  system  match  is  possible. 
This  fact  is  perhaps  best  demonstrated  by  Figures  25  and  26  which 
show  the  cost  functions  versus  ui£S  for  both  flight  conditions. 

The  fact  that  a  better  match  is  possible  with  a  ’’faster”  control 
system  is  not  surprising.  With  sufficiently  fast  control  system 
dynamics,  the  HOS  configuration  responds  more  like  the  basic 
aircraft  (without  control  system);  and,  it  is  precisely  this  basic 
aircraft  form  (i.e.  the  short  period  approximation  in  this  case) 
that  is  being  matched  to  the  high-order  response. 

Figures  27-29  show  that  the  high-order  systems  do  indeed 
behave  more  like  the  basic  aircraft  as  the  control  system  dynamics 
are  made  faster.  Figures  27  and  2£  show  that  the  equivalent  short 
period  parameters,  u)ap  and  ?sp,  approach  their  basic  aircraft 
values  as  W£a  is  increased;  and,  Figure  29  shows  the  expected 
result  that  the  equivalent  time  delay  is  reduced  as  u>£s  is 
increased. 

In  addition  to  increased  mismatch,  another  limitation 
associated  with  coupling  between  the  control  system  and  aircraft 
roots  was  in  the  ’’behavior"  of  the  matching  programs  themselves. 
The  problem,  encountered  during  the  computer  matching,  was  that 
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the  search  vector  (i.e.  the  vector  made  up  of  the  matching 
variables)  stopped  at  ’’local’*  minimums.  Changing  the  initial 
condition  vector  slightly  would  generally  solve  the  problem  by 
allowing  the  search  vector  to  bypass  the  local  minimum  and  arrive 
at  a  smaller  minimum  which  gave  the  desired  best  match.  In  most 
cases,  it  was  apparent  when  a  local  minimum  was  encountered 
because  the  equivalent  system  parameters  would  take  on  extreme 
values,  costf  would  be  unreasonably  high,  and  (in  some  cases)  the 
low-order  equivalent  system  would  have  an  unstable  time  response. 
In  some  cases,  however,  it  was  difficult  to  determine  if  a  local 
minimum  had  been  encountered,  or  if  there  was  some  other  problem 
precluding  a  good  match.  For  configurations  with  coupling  between 
the  control  system  and  aircraft  roots  the  programs  had  a  greater 
tendency  to  arrive  at  a  local  minimum.  In  these  conf igurations, 
the  programs  were  more  sensitive  to  variations  in  the  initial 
condition  vector,  and  the  input  initial  conditions  had  to  be 
closer  to  the  correct  final  values  in  order  to  arrive  at  the  best 
match. 

A  comparison  of  the  La-free  versus  La-fixed  matching  data 
raises  questions  regarding  the  meaning  of  equivalent  La.  Figure 
30  shows  that  the  equivalent  La  approaches  the  aircraft  value  for 
configurations  with  fast  control  system  dynamics,  and  assumes 
considerably  higher  values  for  configurations  with  slow  control 
system  dynamics.  These  results  not  only  raise  questions  regarding 
the  meaning  of  equivalent  La  but  also  raise  questions  regarding 
the  meaning  of  the  other  equivalent  parameters  in  the  same  LOS 
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configuration.  Consider,  for  example,  the  HOS3  configuration  with 
flight  condition  2  and  Wfs  =6.0  (Table  VI  and  Figure  30b).  For 
this  configuration,  the  equivalent  La  of  S.197  is  unreasonably 
high  given  that  the  basic  aircraft  value  is  2. OS.  The  value  of  LQ 
=  8.197  is  a  value  picked  by  the  matching  routine  that  produced 
the  best  mathematical  match  but  is  not  necessarily  representative 
of  the  actual  stability  derivative  of  the  high-order  system.  For 
the  same  HOS  conf iguration,  £Sp  is  nearly  doubled  and  there  are 
significant  changes  in  wap  and  r  between  the  La-fixed  and  LQ-free 
conf igurations.  The  question  is,  which  set  of  parameters  (La, 

£3p,  aiap,  and  r)  actually  characterizes  the  HOS  short  period 
dynamics;  or,  are  the  actual  HOS  characteristics  somewhere  in 
between  the  two  extremes?  The  fact  that  both  sets  of  equivalent 
parameters  give  an  excellent  match  to  the  HOS  indicates  that  they 
both  characterize  the  HOS  dynamics.  If  that  is  the  case,  however, 
which  set  of  parameters  should  be  related  to  the  requirements  of 
MIL-F-8785C  in  order  to  specify  and  predict  flying  qualities? 
Figures  27-30  show  that  this  "dilemma"  is  greatest  (i.e.  there  is 
the  greatest  change  in  equivalent  parameters  between  the  La-fixed 
and  La-free  configurations)  when  the  control  system  dynamics  are 
the  slowest. 

Considerable  insight  could  be  gained  through  a  simulator  or 
flight  test  program  set  up  to  evaluate  a  large  number  of  La-fixed 
and  La-free  configurations.  The  configurations  should  be  similar 
to  the  example  discussed  above  in  which  there  is  a  significant 
change  in  equivalent  parameters  between  the  La-fixed  and  LQ-free 
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Figure  30.  Equivalent  La  versus  u»fa 


configurations.  If  the  results  of  such  evaluations  show  that  over 


a  large  range  of  flight  conditions  and  control  system  dynamics, 
the  HOS  flying  qualities  are  best  characterized  by  the  La-fixed 
parameters,  this  would  lend  credence  to  the  idea  that  La  should  be 
fixed  during  equivalent  system  matching.  On  the  other  hand,  such 
evaluations  could  show  that  La  should  be  free  during  equivalent 
system  matching,  or  that  there  are  certain  configurations  for 
which  La  should  be  free,  etc. 

The  NT-33  piloted  evaluations  conducted  by  Neal-Smith  (11) 
evaluated  a  total  of  51  HOS  configurations,  many  with  significant 
coupling  between  the  control  system  and  aircraft  dynamics.  The 
report  states  that  reasonable  correlation  between  pilot  rating 
data  and  equivalent  systems  data  was  achieved  by  plotting  the 
product  WgpTgg  against  the  time  delay  t.  Due  to  numerous  problems 
with  the  equivalent  systems  approach  however,  the  authors 
developed  a  pi lot-in-the-loop  analysis  that  was  used  instead.  In 
summary  of  the  equivalent  systems  matching  results,  the  report 
states  the  following: 

Although  the  correlation  of  the  PR  data  with  this 
equivalent  system  approach  is  reasonable,  several 
limitations  should  be  noted... In  the  region  of  low  wSpTQ2 
values  (say  wgpTg2  <  2.5),  large  trade-offs  between  Tq2  and 
usp,  or  some  cases  between  ?Sp  and  u)sp,  can  be  made  with 
little  discernible  difference  in  the  accuracy  of  the  analog 
match  achieved.  In  the  present  experiment,  this  problem 
occurred  for  those  configurations  having  0  responses  with 
little  or  no  overshoot  (configurations  with  low  wSp  or 
large  control  system  lags)... The  important  point  to  note  is 
that  this  lack  of  precision  occurs  in  an  area  of  primary 
practical  importance,  since  the  lower  wgpTg2  boundary  would 
represent  a  design  limit  on  aft  center-of-gravity  travel 
for  many  airplanes.  (11) 

Some  of  the  problems  quoted  .in  the  Neal-Smith  report  are 
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similar  to  the  problems  found  in  this  study.  The  La-£ixed  and  La- 
free  configurations  also  had  ” large  trade-offs”  between  the 
parameters  with  "little  discernible  difference  in  the  accuracy  of 
the  analog  match. "  As  in  the  Neal-Smith  study,  these  "trade-offs” 
were  especially  apparent  in  configurations  with  large  control 
system  lags. 

The  Neal-Smith  report  also  states  the  following  with  regard 
to  equivalent  time  delay  determination}  "Accurate  determination 
of  the  time  delay  t  is  often  difficult?  and  for  time  delays 
greater  than  0.1  sec,  small  variations  in  t  can  mean  significant 
changes  in  PR"  (11).  Similarly,  in  the  HOSg  matching  discussed 
above,  t  varied  from  0.168  to  0.220  between  the  LQ-fixed  and  La~ 
free  configurations  (Table  VI)  and  there  is  some  question  as  to 
which  value  is  the  "correct"  one  to  use  for  specifying  flying 
qualities. 

The  discussion  of  these  problems  is  not  meant  to  imply  that 
equivalent  systems  techniques  are  not  suitable  for  specifying 
flying  qualities.  Despite  the  limitations  of  the  equivalent 
systems  method  it  will  continue  to  be  used  in  flying  qualities 
work  and  the  requirements  in  MIL-F-8785C  regarding  equivalent 
systems  are  well  documented.  The  purpose  here  is  to  point  out 
some  of  the  limitations  of  the  method  and  discuss  the  areas  that 
need  more  work.  By  discussing  these  problems  for  both  the 
frequency  domain  and  time  domain  techniques  it  will  be  possible  to 
make  subjective  as  well  as  objective  comparisons  between  the  two 
methods.  The  objective  comparisons  will  be  based  on  mathematical 


differences  (i.e.  mismatch  comparisons),  and  the  subjective 
comparisons  will  be  based  on  the  consideration  of  problems  such  as 
those  discussed  in  this  section. 


81 


IV.  Time  Response  Computer  Matching 


The  time  response  matching  theory  discussed  in  this  chapter 
is  a  least  squares  technique  taken  from  Franklin  and  Powell  (22) 
and  adapted  for  use  in  flight  test  system  identification  work. 

The  final  product  of  the  time  response  matching  is  a  lower  order 
equivalent  system  (LOES)  which  is  a  second  order  approximation  of 
a  higher  order  system  (HOS). 

To  begin  the  time  domain  matching  the  higher  order  system  is 
first  expressed  in  state  space  form  using  vector  notation  as 
f ol lows: 

x  =  Ax  +  Bu  ( 22 ) 

y  =  Cx 

where  A  e  Rn,  B  e  Rm,  and  C  is  chosen  to  give  the  desired  output 
y.  Note  that  n  equals  the  number  of  states  and  m  equals  the 
number  of  inputs  for  the  HOS. 

The  first  step  in  the  time  domain  matching  requires  discrete 
HOS  input  and  output  data.  To  develop  a  discrete  aircraft  model 
from  the  continuous  model,  equation  22  is  integrated  over  a 
sampling  time  T  to  yield  (21): 

x [ ( k+1 )T]  =  F(T)x(kT)  +  G(T)u(kT)  (23) 

where  the  matrices: 

F(T)  =  eAT 

T 

G(T)  =  $eATdTB 
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and  u  is  assumed  constant  over  the  sampling  interval.  Equation  23 
forms  a  discrete  aircraft  model  which  gives  the  aircraft  states  x 
at  time  in  terms  of  the  state  x  and  input  u  at  the  previously 

sampled  time  T^.  Equation  23  will  be  referred  to  in  shortened 
form  as  follows: 

x(  k+1 )  =  Fx(k)  +  Gu(k)  (24) 

y(k+l )  =  Hx( k+1 ) 

where  the  H  matrix  can  be  chosen  to  give  the  desired  output  y. 

One  appealing  aspect  of  the  least  squares  time  response 
matching  technique  is  that  the  discrete  data  can  be  experimentally 
derived.  This  means  that  input  versus  output  time  histories  of 
the  HOS  are  all  that  is  required  to  extract  a  low-order  equivalent 
model.  This  technique,  therefore  would  seem  to  be  well  suited  for 
flight  test  data  collection  and  parameter  identification 
requirements.  For  the  present  development,  however,  the  HOS  has 
been  analytically  modeled  so  that  equation  24  can  be  used  to 
generate  ’’synthetic"  flight  test  data. 

ARMA  Canonical  Form 

At  this  point  it  will  be  necessary  to  introduce  the  Auto 
Regressive  Moving  Average  (ARMA)  canonical  state  space  form.  This 
canonical  representation  has  one  feature  especially  appropriate 
for  identification,  and  will  be  used  to  develop  the  least  squares 
matching  technique  in  the  next  section.  For  a  complete  discussion 
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of  the  ARMA  and  other  canonical  forms,  the  interested  reader  is 
referred  to  Franklin  and  Powell  (22). 

To  introduce  the  ARMA  canonical  form,  consider  the  discrete 
transfer  function  matrix  //(z)  which  can  be  derived  from  the  state 
space  matrices  using  the  following  relation: 

Hi  z)  =  Y(z)/U(z)  =  H(zl  -  F)_1G  (25) 

Equation  25  implies  that  for  a  given  set  of  state  matrices  F, 
G,  and  H  there  is  a  unique  transfer  function  representation  Hiz) . 
Note,  however,  that  for  a  given  transfer  function  matrix  there  is 
not  a  unique  set  of  state  matrices.  In  general,  there  may  be  many 
combinations  of  F,  G,  and  H  that  have  the  same  transfer  function. 
Thus,  the  definition  of  the  state  is  not  unique  and  can  be 
selected  in  such  a  way  which  makes  the  task  of  system 
identification  as  easy  as  possible.  The  standard  way  to  do  this 
is  to  define  the  state  so  that  F,  G,  and  H  are  in  accordance  with 
one  of  the  canonical  forms  for  transfer  functions.  For  the  least 
squares  system  identification  presented  in  the  next  section,  the 
ARMA  canonical  form  will  be  used. 

To  look  at  a  specific  example  of  the  ARMA  representation, 
suppose  the  state  equations  represent  a  single  input/single  output 
(SISO)  system.  The  transfer  function  matrix  would  then  reduce  to 
a  scalar  and,  for  a  third-order  system  may  be  represented  by  the 
following  equation: 
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b-iz  1  +  boz  2  +  boz  3 

tf(z)  =  -  (26) 

1  -  a^z-1  -  a2 z~2  ~  a3z_3 


The  ARMA  canonical  form  for  equation  26  has  six  states  and  is 
given  by  the  matrices: 


~al 

a2 

a3 

bl 

b2 

b3 
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0 
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1 
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0 

0 
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0 
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1 

0 

0 

0 
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0 
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0 
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0 

0 

H 

= 

a2 

a3  b 

1  b2 

b3] 

The  system  described  by  equation  27  is  seen  to  have  six 
states  to  describe  a  third-order  transfer  function  and  thus  to  be 
"nonminiial . ”  In  fact,  the  F  matrix  can  be  shown  to  have  three 
poles  at  z  =  0  that  are  not  observable  for  any  values  of  or  b^ 
However,  the  system  does  have  one  remarkable  property:  the  state 
is  given  by 


xk  =  [yk-i  yk-2  yk-3  uk-i  uk-2  uk-3]T  (28) 

In  other  words,  the  state  is  exactly  given  by  the  past  inputs  and 
outputs.  Thus,  if  the  set  of  uk  and  yk  are  known,  the  state  is 
known  also  since  it  is  merely  a  listing  of  the  six  members  of  the 
set.  All  the  "action”  takes  place  in  the  output  equation  which  is 

yk  «  Hxk 

=  ®iyk-l  +  a2yk-2  +  a3yk-3  +  bluk-l 

+  b2uk-2  +  b3,vik-3  (29) 
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There  is  actually  no  need  to  carry  any  other  equation  along  since 
the  state  equations  are  trivially  related  to  this  output  equation. 


Least  Squares  and  Equation  Error 

Before  discussing  the  actual  least  squares  matching 
technique,  a  definition  of  ’’goodness”  of  fit  will  be  made  using 
the  concept  of  equation  error.  To  begin  with,  a  parameter  vector 
3  will  be  defined  as 

3  =  £®l  a2  a3  ^1  ^2  (30) 

where  the  a.^  and  b.^  are  the  coefficients  of  the  discrete  transfer 
function  (equation  26).  Next,  imagine  that  a  set  of  input  sample 
values  u(k)  and  a  set  of  corresponding  output  sample  values  y(k) 
are  observed,  and  that  these  come  from  a  plant  which  may  be 
described  by  the  transfer  function  (equation  26)  for  some  ’’true” 
value  of  the  parameters,  3°«  The  identification  problem  is  to 
compute  from  these  u(k)  and  y(k)  an  estimate  £  which  is  a  ’’good” 
approximation  of  3°-  To  do  this,  it  is  necessary  to  define  some 
idea  of  ’’goodness  of  fit”  of  a  proposed  value  of  3  to  the  true  3°. 
Because,  by  the  very  nature  of  the  problem,  3°  is  unknown,  it  is 
unrealistic  to  define  a  direct  parameter  error  between  3  and  3°. 
The  error  must  be  defined  in  a  way  which  can  be  computed  from  the 
u(k)  and  y(k).  Three  criteria  which  have  been  proposed  and 
studied  extensively  are  known  as  equation  error,  output  error,  and 
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output-prediction  error. 

Equation  error  (the  method  to  be  used  in  the  present 
development)  requires  the  complete  equations  of  motion  in  a  state 
variable  description.  In  general,  consider  the  continuous  time 
state  description  with  parameter  vector  £  which  can  be  written  as 
fol lows: 


x  =  .f(x,u?0)  (31) 

First,  it  is  assumed  that  the  form  of  the  vector  functions  f  are 
known  but  not  the  parameters  0°,  which  describe  the  plant.  Next, 
it  is  assumed  that  it  is  possible  to  measure  not  only  the  controls 
u  but  also  the  states  x  and  the  state  derivatives  x.  Thus, 
everything  is  known  about  the  equations  except  for  the  particular 
parameter  values.  An  error  can  therefore  be  defined  by  the  extent 
to  which  these  equations  of  motion  fail  to  be  true  for  a  specific 
value  of  0  when  used  with  the  specific  actual  data  xa,  xa,  and  ua. 
The  error  is  written  as 

xa  -  f(xa,ua;e)  =  e(t;0)  (32) 

and  e(t,0°)  =  0  where  0°  are  the  true  plant  parameters.  The 
vector  e(t,0)  are  the  equation  errors,  and  are  used  to  form  a 
nonnegative  cost  function: 

T 

J(0)  =  $eT(t;0)e(t;0)dt  (33) 
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With  the  cost  function  defined,  a  search  over  3  is  made  until 
£  is  found  such  that  J(£)  =  0,  at  which  time  the  parameter  set  (3 
is  equivalent  to  0°.  If  a  unique  parameterization  has  been 
selected,  then  only  one  parameter  set  will  make  e(t,3°)  =  0  and 
thus  @  =  0°. 

The  assumption  that  enough  instruments  are  available  to 
measure  the  total  state  and  all  state  derivatives  is  a  strong 
assumption  and  often  not  realistic  in  continuous  model 
identification.  However,  in  discrete  linear  models  there  is  one 
case  where  it  is  immediate,  and  that  is  the  case  of  an  ARMA  model. 
The  reason  for  this  is  not  hard  to  find:  in  an  ARMA  model  the 
state  is  no  more  than  recent  values  of  input  and  output!  To  be 
more  explicit  about  it,  the  linear  discrete-model  equation  error 
can  be  written: 

xa(k+l)  -  Fxa( k )  -  Gua(k)  =  e(k;3)  (34) 

where  F  and  G  are  functions  of  the  parameters  3.  Substituting  the 
values  from  equation  27,  the  ARMA  model,  gives  the  following: 

X1  al  a2  a3  bl  b2  b3  X1  0 

x2  100000  x2  0 

Xg  0  1  0  0  0  0  Xg  0 

x4  -  0  0  0  0  0  0  x4  -  1  u(k) 

x6  000100  x5  0 

xe  000010  x6  0 

J  k+1  L  J  k 

When  the  substitution  from  equation  28  is  made  it  is  found  that 

for  any  3  (4. e.  for  any  values  of  a^  and  b^ )  the  elements  of 
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equation  error  are  all  zero  except  e^ .  This  element  of  error  is 
given  by: 

x^k+l)  -  a1x1(k)  -  a2x2(k)  -  agxg(k)  -  b1x4Ck) 

-  b2x5(k)  -  bQxe(k)  =  e^k;^)  (36) 


or,  from  equation  28: 

ya(  k )  -  aj^ygtk-l)  -  a2ya(k-2)  -  a3ya(k-3)  -  b1ua(k-l) 

-  b2ua(k-2)  -  b3ua( k-3 )  =  ejJkje)  (37) 

The  discrete  cost  function  becomes: 

N 

y  2 

J(6)  =  ,  ^  ei  (k;6)  (38) 

k=0  x 

Again,  the  subscript  a  has  been  placed  on  the  observed  data 
to  emphasize  the  fact  that  these  are  actual  data  which  were 
produced  via  the  plant  with  parameter  values  8°  and  in  equations 
36-38  an  error  results  because  8  differs  from  6°. 

Least  Squares  Matching  Technique: 

To  begin  the  discussion  of  the  least  squares  matching 
technique,  consider  the  ARMA  model  and  equation  error  which  leads 
to  equation  37,  repeated  below  for  the  nth-order  case.  (The 
subscript  a  is  understood  but  omitted  here). 
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y(k)  -  a^yOc-l)  -  a2y(k-2)  -  •  •  •  -  any(k-n) 

-  b1u(k-l)  -  •  *  *  -  bnu(k-n)  =  e^(k;9)  (39) 

Next,  assume  that  the  following  set  of  inputs  and  outputs  have 
been  observed 

[u(0),u(l),u(2) . u(N),y(0),y(l) . y(N)]  (40) 

from  which  it  is  desired  to  compute  values  for 

9  =  [ax  • • •  an  bx  bn]T  (4l) 

which  will  best  fit  the  observed  data.  Since  y(k)  depends  on  past 
data  back  to  n  periods  earlier,  the  first  error  that  can  be  formed 
is  e(n;9).  Suppose  a  vector  of  errors  is  defined  by  writing 
equation  39  over  and  over  for  k  -  n,  n+1,  N.  The  results 

would  be 

y ( n )  =  xT(n)9  +  e(n;9> 
y(n+l)  -  x^(n+l)9  +  e(n+l;9) 

y(N)  =  xT(N)9  +  e( N { 9 )  (42) 

where  the  definition  of  the  state  of  the  ARMA  model  is 

x(k)  =  [y(k-l)  y( k-2) . . . u( k-1 ) . . . u( k-n  )  1T  (43) 
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To  make  the  error  even  more  compact,  another  level  of  matrix 
notation  is  used  by  defining  the  following: 

Y (N )  =  Cy ( n )  y(n+l). . .y(N)]T 
X(N )  =  Cx( n )  x(  n+1 ) . . . x( N ) 1T 
e(NO)  =  Ce(n)  e( n+1 ) . . . e( N ) ) T  (44) 

Note  that  X(N)  is  a  matrix  with  2 n  columns  and  N-n+1  rows.  In 
terms  of  these  definitions,  the  equation  error  can  be  written  as 

Y  =  X0  +  e(N;£)  (45) 


Least  squares  is  a  prescription  that  one  should  take  that 
value  of  6  which  makes  the  sum  of  the  squares  of  the  e(k)  as  small 
as  possible.  In  terms  of  equation  42  the  cost  function  to  be 
minimized  is  defined  as 


JO) 


2 

e  (kO) 


(46) 


and  in  terms  of  equation  45  this  is 


JO)  =  eT( N O  )e(  N ;  0  )  (47) 


It  is  desired  to  find  @ls»  )eas^  squares  estimate  of  6°, 

which  is  that  0  having  the  property 
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(48) 


J($LS)  $  J(8) 

But  J(8)  is  a  quadratic  function  of  the  2n  parameters  in  8,  and 
from  calculus,  a  necessary  condition  on  that  the  partial 

derivatives  of  J  with  respect  to  8  at  8  a  @^g  should  be  zero. 

This  is  done  as  follows: 

J(8)  *  eTe 

=  Y  -  X8t( Y-X8 ) 

-  yTy  -  0TxTy  _  yTXp  +  8TXTX8  (48) 

and  applying  the  rules  developed  for  derivatives  of  scalars  with 
respect  to  vectors,  the  following  is  obtained: 

jp  =  aj/ae  =  -  2ytx  -  28TxTx  (so) 

Talcing  the  transpose  of  equation  50  and  letting  8  =  @LS’  the 
result  must  be  zero.  This  gives 

XtX@ls  ■  XtY  (51) 

These  equations  are  called  the  "normal  equations”  of  the  problem, 
and  their  solution  will  provide  the  desired  least  squares  @Lg. 

Whether  or  not  the  equations  have  a  unique  solution  depends 
mainly  on  how  8  was  selected  and  what  input  signals  u(k)  were 
used.  Consider  first  the  selection  of  8  and  the  "nonuniqueness” 
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of  state  space  representations.  If  the  state  space  model  of 
equation  24  represented  a  third  order  system,  for  example,  then 
equation  24  would  have  15  parameters  (nine  elements  in  F  and  three 
elements  each  in  G  and  H).  Recall,  however,  that  only  six 
elements  are  needed  to  completely  describe  the  input-output 
dependency  (i.e.  equation  26).  This  means  that  the  state  variable 
description  has  nine  parameters  that  are  in  some  sense  redundant 
and  may  be  chosen  rather  arbitrarily.  If  the  15  element  6  were 
used,  the  resulting  normal  equations  could  not  have  a  unique 
solution.  To  obtain  a  unique  parameter  set  a  canonical  form 
having  a  minimal  number  of  parameters,  such  as  the  observer  or 
ARMA  forms  must  be  selected.  By  way  of  definition,  a  parameter  6 
having  the  property  that  one  and  only  one  value  of  6  makes  J(6)  a 
minimum  is  said  to  be  '‘identifiable.  ”  Two  parameters  having  the 
property  that  J(@^)  =  Jf^g)  are  said  to  be  "equivalent.” 

As  to  the  selection  of  the  inputs  u(k),  consider  first  an 
extreme  case.  Suppose  u(k)  =  c,  a  constant,  for  all  k  (i.e.  a 
step  function  input).  Now  consider  equation  42  again  for  the 
third  order  case  to  be  specific.  The  errors  are 

y(3)  -  ajy(2)  -  a2y(l)  -  a3y(0)  -  b^  -  b2c  -  b3c  =  e(3) 

(52) 

y(N)  -  a^y(N-l)  -  a2y(N-2)  -  a3y(N-3)  -  b^c  -  b2c  -  b3c  =  e(N) 

It  is  obvious  that  in  equation  52  the  parameters  b^,  b2,  and  b3 
always  appear  as  the  sum  b^+b2+b3  and  that  separation  of  them  is 
not  possible  when  a  constant  input  is  used.  Somehow  the  constant 
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u  fails  to  ’‘excite”  all  the  dynamics  of  the  plant.  This  problem 
has  been  studied  extensively  by  Astrom  and  Boh  1 in  (23)  and  the 
property  of  ’’persistently  exciting”  has  been  defined  by  them  to 
describe  a  sequence  of  u(k)  which  fluctuates  enough  to  avoid  the 
possibility  that  only  linear  combinations  of  elements  of  0  will 
show  up  in  the  error  and  normal  equations.  Without  going  into 
more  detail  at  this  point,  an  input  is  said  to  be  persistently 
exciting  if  the  lower  right  (n  x  n)-matrix  component  of  X^X  (which 
depends  only  on  u(k))  is  nonsingular. 

Assuming  that  the  u(k)  are  persistently  exciting  and  that  the 
B  are  identifiable  (and  consequently  that  XTX  is  nonsingular), 
then  the  explicit  least  squares  solution  can  be  written  as 
f ol lows: 

Bls  =  (xTxi-ixTy  (53) 

Extracting  the  Continuous  Model: 

The  parameters  in  Bls  describe  a  low-order  discrete  model 
(i.e.  equation  26)  of  the  high-order  system.  The  next  step  in  the 
matching  process  is  the  extraction  of  a  low-order  continuous  model 
from  the  discrete  model.  The  discrete  transfer  function 
represented  by  and  equation  26  is  first  converted  into  a 

discrete  state  space  form  as  follows: 


al 

1 
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•The  tilda  denotes  the  Low-order  equivalent  state  matrices. 
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Equation  54  is  known  as  the  observer  canonical  form.  Note, 
however,  that  any  state  space  representation  of  the  transfer 
function  could  be  used  at  this  point.  This  is  because  the  final 
step  will  be  a  conversion  back  to  a  continuous  transfer  function 
form.  Thus,  the  state  space  form  (which  is  not  unique)  is  used 
only  as  an  intermediate  step  to  obtain  the  final  transfer  function 
form  (which  is  unique).  Since  the  final  transfer  function  is 
unique,  any  state  space  representation  should  give  the  same 
result. 


To  obtain  the  continuous  transfer  function  form,  the  discrete 
state  space  model  of  equation  54  must  first  be  converted  to  a 
continuous  state  space  model.  A  straightforward  and  accurate 
method  for  this  has  been  developed  by  Sinha  and  Lastman  (24).  The 


technique  begins  with  an  initial  estimate  A0  for  the  A  matrix 
based  on  the  derived  F  matrix: 


/v  'v  - 

AqT  =  0.5(F-F_1) 


(55) 


The  technique  then  begins  its  iterative  phase  by  computing  an 
estimated  F  matrix  Fj  as  follows: 


F  4  =  exp(A-iT) 


(56) 


where  j  is  the  iteration  counter.  Note  that  Fj  is  what  the  F 
matrix  should  be  if  A  =  Aj.  The  next  iterative  step  is  to  correct 
the  estimate  of  AjT  for  the  difference  between  F  and  Fj. 
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(57) 


/V  A/  «  A/  /W 

Ad+1T  =  AjT  +  F-l(F-Fj) 

The  technique  iterates  between  these  last  two  equations  until 
Aj+^T  is  sufficiently  close  to  AjT.  The  estimated  stability 
matrix  A  is  then  Aj+^.  The  estimated  control  matrix  B  is  found 
from: 


B  =  C(Aj+1)"1(exp(Aj+1T)_I]“1G  (58) 

The  Sinha  and  Lastman  technique  provides  extremely  accurate 
results  in  extracting  the  continuous  A  and  B  matrices  for  most 
applications.  For  the  technique  to  work  however,  the  inverted 
matrices  in  equations  57  and  58  must  be  full  rank.  For  this 
study,  a  trial  run  using  the  known  short  period  stability 
derivatives  for  the  A-4D  resulted  in  accuracies  of  10-14  for  the 
A  matrix  but  only  about  10~®  for  the  B  matrix.  It  was  found  that 
the  problem  was  in  the  estimated  stability  matrix  which  was  close 
to  being  singular  (determinant  close  to  zero)  and  equation  58 
required  that  this  matrix  (Aj+^)  be  inverted.  Because  of  this 
problem,  an  alternative  method  for  extracting  the  control  matrix 
was  found. 

In  cases  where  the  stability  matrix  is  singular,  an 
alternative  technique  is  to  extract  the  control  matrix  based  on 
the  power  series  expansion  of  the  matrix  exponential  function  eA^ 
which  is  given  by 
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eAT  =  I  +  AT  +  (A2T2)/2!  +  •••  (59) 

The  technique  is  based  on  the  definition  of  the  discrete  control 
matrix  6.  From  equation  23  we  have 

~  T  ~ 

G  =  $eATdTB  (GO) 

o 

which  can  be  written  as  a  power  series  using  equation  59  as 
fol lows: 

T 

G  =  ^(1  ♦  At  +  A2t2/2 !  +  *  •  •  )dT*B  (61) 

o 

After  integrating  and  inverting,  we  have  the  following: 

B  =  (T  +  AT2/2  +  A2T3/3*2  +  • • ♦ )-1G  (62) 

This  gives  a  method  for  estimating  the  control  matrix  B  which 
does  not  require  the  stability  matrix  A  to  be  full  rank.  Thus,  by 
substituting  the  estimated  A  matrix  obtained  from  the  Sinha  and 
Lastman  technique  into  equation  62  the  estimated  control  matrix  B 
was  obtained  for  the  A-4D  short  period  example  discussed  above. 

In  addition,  because  the  power  series  can  be  carried  out  as  far  as 
desired,  the  accuracy  of  this  method  exceeds  the  accuracy  of  the 
Sinha  and  Lastman  technique.  For  the  A-4D  short  period  example 
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the  B  matrix  was  estimated  with  an  accuracy  of  10“14,  a 
significant  improvement  over  the  previous  results.  For  this 
reason,  the  Sinha  and  Lastman  technique  was  used  to  extract  the  A 
matrices  and  the  alternative  technique,  based  on  the  power  series 
expansion  of  was  used  to  extract  the  B  matrices  throughout 

the  least  squares  matching. 

With  the  continuous  A  and  B  matrices  extracted  from  the 
discrete  model,  we  have  the  desired  LOES  in  state  space  form: 

x  =  Ax  +  Bu  (63) 

y  =  Cx 

~  /v 

where  A  e  Rn,  B  e  Rm,  and  C  is  chosen  to  give  the  desired  output 
y.  Note  that  n  equals  the  number  of  states  and  m  equals  the 
number  of  inputs  for  the  LOES. 

The  final  step  is  to  obtain  the  desired  LOES  transfer 
functions  (0/Fa,  nz/Fa,  etc)  from  the  transfer  function  matrix 
H(s): 


H(s)  *  Y(s)/U(s)  =  C(sl  -  A)-1B  (64) 

The  order  of  the  least  squares  LOES  is  determined  by  the 
number  of  columns  of  the  X  matrix  (see  equations  43  and  44)  which 
are  made  up  of  the  observed  inputs  u(k)  and  outputs  y(k).  The 
column  size  of  this  matrix  determines  the  size  of  the  resulting 
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parameter  vector  B^g  which  in  turn  determines  the  order  of  the 
final  LOES  transfer  function. 

Optimal  Input  for  System  Identification: 

The  studies  by  Astrom  and  Bohlin  (23)  have  shown  that  the 
choice  of  the  input  excitation  u  is  an  important  consideration  in 
time  response  matching  techniques.  In  general,  it  can  be  said 
that  if  the  input  fails  to  "excite"  all  the  dynamics  of  the  plant 
it  will  be  difficult  or  even  impossible  to  identify  the  plant 
correctly.  As  discussed  in  the  previous  section,  Astrom  and 
Bohlin  have  defined  an  input  to  be  "persistently  exciting”  if  the 
matrix  X^X  (which  depends  only  on  u(k))  is  nonsingular.  During 
the  least  squares  matching  of  this  study  the  inputs  used  to  excite 
the  plant  were  persistently  exciting  by  this  definition.  The 
question  arose,  however,  as  to  what  the  optimal  input  would  be. 
That  is,  of  the  types  of  inputs  that  qualify  as  "persistently 
exciting"  are  there  some  that  will  excite  the  plant  better  than 
others  and  thereby  achieve  a  better  identification  of  the  plant? 

To  determine  the  optimal  input  excitation  to  be  used  for  the 
least  squares  matching  a  variety  of  ramp  inputs  were  used. 
Variations  were  made  in  the  rise  time,  amplitude,  and  sinusoidal 
components  of  the  input  to  determine  the  effect  of  each  of  these 
on  the  quality  of  the  least  squares  match.  Figure  31  illustrates 
two  of  the  inputs  used  and  the  corresponding  plant  responses. 

Note  that  100  discrete  input  and  output  data  points  were  used 
throughout  the  LS  matching  (i.e.  N  *  100  in  equations  40-52). 
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The  range  of  rise  times  used  varied  from  N  =  0  (step  input)  to 
N  =  50.  Since  the  time  step  T  =  0. 1  sec,  this  corresponds  to  a 
variation  in  rise  time  of  0  to  5.0  seconds.  The  amplitude  of  the 
ramp  input  was  varied  from  0.1745  to  10.0  and  the  sinusoidal 
component  was  varied  from  0  to  0 . 1667sin( lOt ) .  (The  sinusoidal 
frequency  was  held  constant  at  10  rad/sec  while  the  sinusoidal 
amplitude  was  varied). 

As  the  rise  time,  amplitude,  and  sinusoidal  components  of  the 
input  were  varied,  corresponding  changes  in  the  frequency  domain 
and  time  domain  cost  functions  (costf  and  cost^. )  were  used  to 
determine  the  quality  of  the  match.  The  input  which  produced  the 
lowest  mismatch  as  defined  by  these  functions  was  considered  to  be 
the  ’’optimal”  input.  To  investigate  the  effects  of  rise  time, 
amplitude,  and  sinusoidal  components  one  of  the  three  elements  was 
varied  while  the  other  two  were  held  constant  (i.e.  rise  time  was 
varied  from  0  to  5.0  seconds  while  the  amplitude  and  sinusoidal 
components  were  held  constant  at  some  value,  etc.).  In  general, 
it  was  found  that  variations  in  amplitude  and  sinusoidal 
components  had  little  effect  on  the  quality  of  the  LS  match. 
Changes  in  the  rise  time,  however,  resulted  in  significant  changes 
in  the  cost  functions.  Figure  32  shows  the  variation  in  cost^  and 
costt  as  a  function  of  rise  time  (for  an  A-4D  H0S)  and  indicates 
that  there  is  indeed  an  optimal  rise  time  associated  with  this 
type  of  least  squares  matching. 
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The  fact  that  there  is  an  optimal  rise  tine  for  the  LS 
matching  technique  can  perhaps  be  explained  by  second-order 
dynamic  response  theory.  Craig  (25),  shows  that  a  second-order 
system  will  respond  quite  differently  to  ramp  inputs  depending  on 
the  rise  time.  Figure  33  shows  the  dynamic  response  of  a  second 
order  system  using  rise  times  of  tr  =  0.2Tn  and  tr  =  1.5Tn  where 
Tn  is  the  undamped  natural  period.  The  dynamic  effects  can  be 
ignored  (and  the  load  is  considered  "statically  applied")  if  the 
rise  time  is  longer  than  about  3Tn. 

It  makes  sense  that  if  an  input  is  "slowly  applied”  the 
dynamics  of  the  plant  will  not  be  excited  sufficiently  to  identify 
the  system.  On  the  other  hand,  it  has  been  shown  by  Franklin  and 
Powell  (22)  that  for  extremely  small  rise  times  (i.e.  step  input) 
the  dynamics  of  the  plant  are  also  not  excited  sufficiently  for 
system  identification  (see  equation  52  and  associated  discussion). 
With  these  considerations  in  mind,  it  makes  sense  that  there  would 
be  an  optimal  rise  time  between  the  two  extremes. 


Figure  33.  Response  of  Second  Order  System  to  Ramp  Inputs 
(Reference  25) 
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Calculation  of  Required  LQES  Parameters 

Once  the  LOES  transfer  function  (equation  64)  is  obtained, 
the  equivalent  parameters  La,  £sp,  uSp,  and  T3  are  extracted  and 
compared  with  the  corresponding  parameters  obtained  from  the 
frequency  response  matching.  The  extraction  of  the  parameters  La, 
X gp t  and  wSp  is  straightforward  since  these  are  just  the  roots  of 
the  LOES  numerator  and  denominator.  The  determination  of 
equivalent  time  delay  T3  however  is  somewhat  more  complicated. 
Since  there  is  not  a  straightforward  method  for  this  in  the  time 
domain;  Tq  was  determined  by  iteration  in  the  frequency  domain. 
This  was  accomplished  by  using  the  HOS  and  LOES  transfer  functions 
to  calculate  phase^s  and  phase^S*  The  iteration  was  done  by 
subtracting  57. 29578*w from  phaseLos  on  each  iteration  until 
the  cost^  function  was  minimized.  (Note  that  this  iteration  has 
no  effect  on  the  gain  component  of  the  cost  function).  The  change 
in  phase  before  and  after  this  iteration  in  time  delay  is  shown  in 
Figure  34  for  an  aircraft  configuration  used  in  the  flight  test 
program  (explained  in  chapter  V). 


Least  Squares  Watching  Results 

To  compare  the  time  response  matching  results  of  this  chapter 
with  the  frequency  response  matching  results  of  chapter  III,  the 
LS  technique  was  used  to  identify  low-order  equivalent  systems  for 
the  12  0/Fs  high-order  systems  used  in  chapter  III.  In  the 
frequency  response  matching,  the  12  0/Fa  HOS  configurations 


(derived  by  varying  ui£a  from  6.0  to  31.0  rad/sec  for  flight 
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Figure  34. 


Calculation  of  Equivalent  Time  Delay 
(Configuration  3-1) 


conditions  1  and  2)  were  all  matched  with  La  fixed  and  free 
resulting  in  a  total  of  24  LOES.  The  LS  matching  is  done  in  the 
time  domain,  however,  and  there  is  no  way  to  fix  La  at  a  specified 
value.  Thus,  only  the  La-free  LOES  configurations  will  be 
compared  between  the  two  matching  techniques. 

Based  on  the  analysis  of  "optimal”  LS  input  excitations 
(previous  section)  a  ramp  input  was  used  with  a  rise  time  of  1.1 
sec  for  the  flight  condition  1  matching  and  a  rise  time  of  0.5  sec 
for  the  flight  condition  2  matching.  The  amplitude  of  the  ramp 
input  was  5.0  and  there  was  no  sinusoidal  component.  The  results 
of  the  LS  matching  are  presented  in  tabular  format  in  Appendix  A 
and  various  graphical  results  are  presented  in  this  discussion. 
Analysis  of  the  cost  function  data  (Figures  35  and  36)  shows  that 
the  LS  matching  generally  had  less  mismatch  in  the  time  domain 
(lower  costf. )  but  greater  mismatch  in  the  frequency  domain  (higher 
costf)  when  compared  to  the  LONFIT  results.  This  makes  sense 
intuitively  since  the  LS  match  is  performed  using  HCS  time 
response  data  and  the  LONFIT  match  is  performed  using  HOS 
frequency  response  data.  It  makes  sense,  therefore,  that  the  LS 
match  is  closer  in  the  time  domain  and  the  LONFIT  match  is  closer 
in  the  frequency  domain.  A  visual  comparison  of  the  frequency 
domain  and  time  domain  mismatch  for  HOS^  and  HOS3  is  shown  in 
Figures  37-39  and  can  be  compared  with  the  corresponding  LONFIT 
results  in  chapter  III. 

Comparison  of  the  equivalent  parameters  (Figures  40-43)  show 
that  both  the  LS  and  LONFIT  parameters  approach  their 
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corresponding  basic  aircraft  values  as  u>£a  is  increased.  This  is 
expected  as  the  control  system  dynamics  are  made  faster.  As  ui£S 
is  increased,  the  effects  of  the  control  system  become  less 
significant  and  the  HOS  behaves  more  like  the  basic  aircraft  (see 
the  discussion  on  the  effects  of  control  system  dynamics  in 
chapter  III).  Figures  40-43  show  that  the  LS  parameters  are 
closer  to  the  basic  aircraft  values  for  La  and  tq  while  the  LONFIT 
parameters  are  closer  for  wap  and  £ap.  The  LS  equivalent  control 
anticipation  parameters  (CAP)  are  closer  to  the  basic  aircraft 
values  for  flight  condition  1  but  the  LONFIT  results  are  closer 
for  flight  condition  2  (Figure  44). 

These  comparisons  between  the  Least  Squares  time  response 
matching  program  and  the  LONFIT  frequency  response  matching 
program  have  shown  favorable  results  for  the  LS  program.  These 
comparisons  primarily  focused  on  the  equivalent  parameters 
generated  by  the  two  programs  for  a  given  high  order  system.  The 
equivalent  parameters  were  numerically  similar  between  the  two 
programs  and  approached  the  basic  aircraft  values  as  the  feel 
system  natural  frequency  was  increased.  This  is  consistent  with 
equivalent  systems  theory.  In  the  next  chapter  the  Least  Squares 
program  will  be  further  validated  and  compared  to  the  well- 
accepted  "Bandwidth'*  method  (described  in  Appendix  H)  by 
predicting  flying  qualities  of  higher  order  systems.  This 
comparison  will  be  made  based  on  the  ability  of  the  different 
techniques  to  accurately  predict  the  handling  quality  levels  of 
higher  order  systems  using  present  MIL-STD-1797  guidance. 
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Figure  35.  Cost  Functions  versus  ufa  for  HOS^  Matchi 
(Least  Squares/LONFIT  Comparison) 
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Figure  36.  Cost  Functions  versus  u>fa  for  HOS3  Matchi 
(Least  Squares/LONFIT  Comparison) 
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Figure  40.  Equivalent  u»ap  versus  ui£a 

(Least  Squares/LONFIT  Comparison) 


TAU  (SEC)  TAU  (SEC) 


(b)  Flight  Condition  2 

Figure  42,  Equivalent  Time  Delay  (tq)  versus  W£s 
(Least  Squares/LONFIT  Comparison) 
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Figure  44.  Equivalent  CAP  versus  W£S 

(Least  Squares/LONFIT  Comparison) 


V.  Flight:  Test  Results 


A  limited  flight  test  program  was  conducted  at  the  Air  Force 
Flight  Test  Center  (AFFTC),  Edwards  Air  Force  Base,  California,  to 
evaluate  the  Least  Squares  system  identification  program  developed 
in  this  thesis  and  to  evaluate  an  additional  AFIT  thesis  project 
involving  the  prediction  of  PIOs  and  handling  qualities  using  the 
Optimal  Control  Model  (OCM).  The  Commandant,  USAF  Test  Pilot 
School  (USAFTPS)  directed  the  test  program  in  partial  fulfillment 
of  the  USAFTPS  curriculum  requirements.  Three  pilots  and  two 
engineers  from  USAFTPS  Class  89A  comprised  the  test  team.  The 
USAFTPS  Special  Projects  Division  acted  as  the  System  Program 
Office.  Twenty-five  sorties  totaling  27.8  flight  hours  were  flown 
between  12  September  and  16  October  1989.  A  sortie  summary  is 
given  in  Appendix  C. 

The  conclusions  of  the  test  team  regarding  the  Least  Squares 
Lower  Order  Equivalent  System  matching  results  will  be  summarized 
in  this  chapter.  A  more  complete  description  of  the  flight  test 
program  can  be  found  in  the  NT-33  HAVE  CONTROL  final  flight  test 
report  (26). 

Objectives 

The  objectives  of  the  HAVE  CONTROL  flight  test  program 
pertaining  to  the  Least  Squares  LOES  technique  were: 

1.  Determine  if  the  LOES  time  response  matching  technique  is 
suitable  for  flight  test  applications.  This  evaluation  included 
ease  of  application  and  comparison  with  Hoh’s  Bandwidth  Method. 
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2.  Determine  the  optimal  ramp  input  rise  time  for  the  LOES 
time  response  matching  program. 

3.  Determine  the  effects  of  control  system  root  location  on 
system  identification  using  the  LOES  method. 

Test  Aircraft  Description 

The  NT-30A  variable  stability  test  aircraft,  S/N  51-4120,  was 
a  modified,  two  seat  Jet  trainer  operated  by  the  CALSPAN 
Corporation,  Buffalo,  New  York  and  owned  by  the  Flight  Dynamics 
Laboratory,  Wright-Patterson  AFB,  Ohio  (27,28).  The  aircraft  was 
capable  of  variable  dynamic  response  and  control  system 
characteristics  (29).  The  Variable  Stability  System  (VSS) 
modified  the  static  and  dynamic  responses  of  the  basic  NT-33A  by 
commanding  control  surface  positions  through  full  authority 
electro-hydraulic  servos.  A  programmable  analog  computer, 
associated  aircraft  response  sensors,  control  surface  servos,  and 
an  electro-hydraulic  force-feel  system  provided  the  total 
simulation  capability.  The  instructor/saf ety  pilot  varied  the 
computer  gains  through  controls  located  in  the  rear  cockpit, 
allowing  changes  in  airplane  dynamics  and  control  system 
characteristics  during  the  flight.  Appendix  D  contains  additional 
information  concerning  the  aircraft  systems,  capabilities,  and 
safety  provisions. 

The  front  cockpit  AVQ-7  Heads  Up  Display  (HUD)  displayed 
several  flight  parameters,  including  airspeed,  altitude,  angle  of 
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attack,  pitch  attitude,  heading,  and  the  velocity  vector.  The  HUD 
was  used  during  the  test  to  closely  simulate  a  representative 
fighter  aircraft. 

Instrumentation  and  Data  Reduction 

The  NT-33  test  instrumentation  system  contained  the  following 
items: 

1.  An  on-board  Ampex  AR  700  magnetic  tape  recording  system 
with  2.25  hours  recording  capability.  This  system  was  used  to 
record  aircraft  flight  conditions,  flight  control  positions,  pilot 
voice,  and  aircraft  states  from  the  aircraft  data  acquisition 
system  (DAS). 

2.  An  AN/ANH-2  voice  recorder  was  used  to  record  interphone 
and  UHF  radio  communications. 

3.  A  HUD  video  recorder  was  used  to  record  all  approaches 
and  landings. 

The  project  pilot  operated  the  HUD  and  on-board  voice 
recorder  system.  The  instructor /safety  pilot  operated  the 
magnetic  tape  system  and  the  HUD  camera.  The  AFFTC  photographic 
branch  provided  ground  videotape  coverage  of  each  landing  task. 

Following  each  NT-33  mission,  the  project  pilots  reviewed  the 
HUD  video  and  tape  recorder  audio  and  summarized  their  comments 
for  each  conf iguration  flown.  Project  pilot  comments  were  used  to 
qualitatively  describe  the  aircraft  PIO  tendencies  and  handling 
qualities  during  the  approach  and  landing  task  (this  task  is 
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defined  in  the  next  section).  In  addition,  pilot  comments  were 
used  to  ensure  project  pilots  used  similar  criteria  when  assigning 
pilot  handling  qualities  ratings  ( PHQR ) .  These  pilot  ratings  were 
given  in  accordance  with  the  Cooper-Harper  scale  as  described  in 
reference  30.  The  PHQRs  for  each  configuration  flown  are 
tabulated  in  the  results  and  analysis  section  of  this  chapter. 

In  accordance  with  the  Cooper-Harper  scale,  the  PHQRs  were 
used  to  form  three  levels  of  flying  qualities.  Level  1  consisted 
of  all  PHQR  values  from  one  to  three.  Level  2  consisted  of  PHQRs 
from  four  to  six.  Level  3  consisted  of  PHQRs  from  seven  to  ten. 
The  flying  qualities  levels  assigned  by  the  pilots  were  compared 
to  the  levels  predicted  prior  to  flight  by  the  Least  Squares  LOES 
and  Bandwidth  methods.  The  purpose  of  using  the  two  prediction 
methods  (Least  Squares  LOES  and  Bandwidth)  was  primarily  to 
compare  results  of  the  two  methods  and  also  to  verify  level 
boundaries  as  contained  in  MIL-STD-1797 . 

After  each  flight,  the  NT-33A  AR-700  magnetic  tape  data  was 
used  to  accomplish  LOES  matching  using  actual  flight  test  data. 

The  results  of  this  matching  were  compared  to  the  actual  pilot 
ratings  to  evaluate  the  suitability  of  the  LOES  program  to 
accomplish  equivalent  system  matching  using  flight  test  data. 

Teat  Methods  and  Conditions 

The  landing  longitudinal  PIO  tendencies,  flying  qualities, 
and  the  Least  Squares  LOES  identification  technique  were  evaluated 
at  three  pairs  of  short  period  natural  frequencies  and  damping 
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ratios.  All  short  period  dynasics  aet  MIL-STD-1797  Level  1 
requirements  for  the  landing  approach  (Category  C).  The  three 
pairs  of  baseline  dynamics  are  depicted  in  Figure  46  and  listed  in 
Table  VIII.  Also  shown  in  Table  VIII  are  the  first  and  second 
order  filters  that  were  used  to  vary  the  configuration  dynamics. 
The  dynamics  were  varied  to  achieve  predicted  handling  qualities 
to  span  the  spectrum  from  Level  1  to  Level  3.  Table  VIII  shows 
the  OCM  predicted  handling  qualities  for  each  of  the  13 
configurations  tested. 
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Figure  46.  NT-33A  Baseline  Dynamics 


The  configuration  dynamics  were  verified  by  the  HAVE  CONTROL 
test  team  and  by  the  CALSPAN  Corporation.  These  verifications  are 
shown  in  Appendix  E.  Only  the  short  period  and  filter  dynamics 
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were  varied  to  produce  different  aircraft  handling  qualities.  The 
phugoid,  lateral-directional,  servo-actuator,  and  feel  system 
dynamics  were  held  constant.  These  parameters  are  shown  in 
Appendix  F. 

After  takeoff,  the  project  pilot  took  control  of  the 
aircraft,  and  climbed  to  5000  feet  pressure  altitude  on  a  North 
turnout.  The  instructor /safety  pilot  set  the  short  period  and 
filter  dynamics  by  adjusting  the  appropriate  variable  stability 
gain  control  in  the  rear  cockpit.  The  project  pilot  established 
the  landing  configuration  and  accomplished  the  auto-step  and  auto¬ 
ramp  inputs  and  the  hand-step  and  hand-ramp  inputs  used  for  the 
system  identification  task.  Further  discussion  of  the  system 
identification  task  is  presented  in  the  next  section. 

After  accomplishing  the  open-loop  tasks,  the  project  pilot 
established  a  1000  feet  per  minute  descent  in  the  landing 
configuration.  Then  at  50  feet  above  a  4000  feet  mean  sea  level 
target  altitude,  the  pilot  simulated  a  landing  task  by  an 
aggressive  level  off.  As  a  safety  precaution,  if  the 
configuration  exhibited  a  divergent  FI0  or  other  Level  3 
characteristics,  that  configuration  was  not  tested  any  further. 

After  the  simulated  landing  was  accomplished,  the  project 
pilot  returned  to  the  pattern  and  flew  the  approach  and  landing 
without  an  offset.  After  touching  down,  the  instructor /safety 
pilot  disengaged  the  VSS  and  performed  the  takeoff.  The  project 
pilot  then  made  preliminary  comments  on  the  configuration  while 
the  instructor /safety  pilot  flew  the  aircraft  on  an  extended 
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downwind  (South  re-entry).  If,  during  the  straight  in  approach,  a 
divergent  PIO  occurred  or  adequate  performance  could  not  be 
achieved,  then  the  offset  landing  task  was  not  attempted.  Two 
visual  approaches  with  a  lateral  offset  were  then  flown,  with  one 
offset  to  each  side  of  the  runway.  After  the  first  offset 
approach,  the  project  pilot  added  to  his  preliminary  comments. 
After  the  second  offset  approach,  the  project  pilot  summarized  his 
overall  comments  and  assigned  a  PIO  rating  and  a  PHQR  for  that 
configuration.  If  the  evaluation  pilot  felt  confident  enough  to 
make  an  overall  evaluation  based  on  only  two  approaches  he  was 
allowed  to  eliminate  the  third  approach.  The  evaluation  pilot  was 
allowed  to  assign  separate  ratings  for  the  approach  and  flare  if 
he  deemed  it  necessary. 

The  landing  task  was  a  visual  approach  with  a  lateral  offset 
and  a  correction  to  centerline  prior  to  touchdown.  The  size  of 
the  lateral  offset  was  approximately  150  feet.  Since  the  width  of 
the  main  runway  at  Edwards  is  300  feet,  the  150  foot  offset  to  the 
left  was  made  by  aligning  with  the  left  edge  of  the  runway,  and 
the  150  foot  offset  to  the  right  was  made  by  aligning  with  the 
right  edge  of  the  runway.  The  aircraft  was  flown  on  glidepath 
using  the  instrument  landing  system  until  the  beginning  of  the 
overrun.  The  correction  to  centerline  was  initiated  at  100  feet 
above  ground  level.  The  safety  pilot  assisted  in  maintaining  a 
constant  offset  correction  between  the  three  project  pilots. 

The  touchdown  zone  was  1000  feet  long  starting  at  500  feet 
past  the  threshold.  The  touchdown  aimpoint  was  1000  feet  from  the 
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threshold  and  within  5  feet  of  centerline.  Each  landing  was 
treated  as  a  ’’must  land”  situation,  unless  the  instructor/saf ety 
pilot  and/or  project  pilot  determined  that  safety  of  flight  would 
be  compromised  in  an  attempt  to  land.  Table  IX  summarizes  the 
performance  criteria  used  to  assign  PHQRs  to  this  visual  landing 
task. 


Table  IX 


Offset  Landing  Task  Performance  Standards 


Desired 

No  PIOs 

Touchdown  within  5  feet  of 
centerline  (main  wheels 
on  centerline) 

Touchdown  aimpoint  +/-  250  ft 
Approach  airspeed  +/-  5  kts 


Adequate 


Touchdown  within  25  feet  of 
centerline  (tip  tank  on 
centerline ) 

Touchdown  aimpoint  +  /-  500  ft 
Approach  airspeed  +10/-5  kts 


System  Identification  Tasks 


The  dynamics  of  each  configuration  flown  in  the  approach  and 
landing  task  were  verified  using  the  Least  Squares  Lower  Order 
Equivalent  System  matching  technique.  This  required  a  ramp  input 
with  an  optimal  rise  time.  The  optimal  rise  time,  a  function  of 
aircraft  dynamics  and  flight  condition,  was  calculated  for  each 
configuration  prior  to  conducting  the  flight  tests.  To 
investigate  the  effect  of  rise  time  on  the  matching  technique, 
several  auto-ramp  and  hand-ramp  elevator  inputs  with  varying  rise 
times  were  used  for  each  change  in  aircraft  dynamics.  In 
addition,  auto-step  and  hand-step  inputs  were  used  for  each  set  of 
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aircraft  dynamics.  The  step  inputs  were  used  in  the  Least  Squares 


program  as  low  rise  time  ramp  inputs  because  in  actual  test  data  a 
step  input  has  a  finite  rise  time. 

The  ramp  and  step  inputs  were  accomplished  with  the  aircraft 
trimmed  for  straight  and  level  flight.  The  hand-step  input  was 
done  by  the  project  pilots  by  first  trimming  the  aircraft  for 
straight  and  level  flight,  then  running  the  trim  nose  down 
(approximately  5  lbs  Fs),  and  releasing  the  stick.  The  hand-ramp 
input  was  accomplished  like  the  hand  step  input  except  that  the 
stick  was  eased  slowly  to  the  trim  condition  using  a  predetermined 
(one  or  two  second)  rise  time.  The  pilots  were  instructed  that 
the  exact  rise  time  was  not  as  critical  as  accomplishing  a  smooth 
release  of  the  stick. 

The  instructor/safety  pilot  ran  the  data  acquisition  system 
(DAS)  during  these  inputs.  The  longitudinal  control  force,  and 
pitch  rate  response  (obtained  from  the  DAS)  were  used  in  the  Least 
Squares  LOES  program  to  extract  a  lower  order  equivalent  system 
model.  The  results  of  the  LOES  matching  using  the  actual  flight 
test  data  were  then  compared  to  results  obtained  prior  to  the 
flights  based  on  the  theoretical  configuration  dynamics.  These 


results  are  presented  in  the  remaining  sections  of  this  chapter. 


Results  and  Analysis 


The  Least  Squares  system  identification  program  was  used  to 
accomplish  lower  order  equivalent  system  (LOES)  matching  with  the 
13  aircraft  configurations  shown  in  Table  VIII.  The  matching  was 
conducted  prior  to  the  actual  flights  using  the  known  NT-33A 
phugoid,  short  period,  servo-actuator,  feel  system,  and  filter 
dynamics.  Only  the  short  period  and  filter  dynamics  were  varied 
from  one  configuration  to  another.  The  pre-flight  matching 
results  were  used  to  predict  handling  qualities  based  on  MIL-STD 
guidance  and  compared  with  the  actual  pilot  ratings  given  during 
the  offset  landing  task.  These  results  were  also  compared  to 
results  obtained  using  Hoh ’ s  Bandwidth  technique  (described  in 
Appendix  H)  and  to  the  LOES  results  using  actual  flight  test  data. 

The  Bandwidth  method  was  used  so  that  the  Least  Squares  LOES 
program  could  be  compared  with  a  well-accepted  handling  qualities 
estimation  technique.  The  comparison  between  the  two  techniques 
was  based  primarily  on  the  number  of  configurations  correctly 
predicted  (compared  with  the  pilot’s  actual  ratings).  The  results 
of  the  actual  flight  test  LOES  matching  were  important  to  verify 
the  analytical  LOES  results  and  demonstrate  the  ability  of  the 
program  to  work  using  discrete  flight  test  data. 

Pre-Flight  Matching  Results 

The  Least  Squares  LOES  matching  technique  (described  in 
detail  in  chapter  IV)  was  used  to  produce  a  least  squares  match  to 
the  higher  order  aircraft  dynamics.  The  higher  order  dynamics 
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were  represented  by  the  aircraft's  pitch  rate  time  response 
excited  by  a  ramp  input.  After  the  program  accomplished  the  time 
response  matching,  frequency  response  lower  order  transfer 
functions  were  extracted  in  the  form  shown  in  equation  15. 

Time  responses  and  frequency  responses  of  the  LOES  model  were 
compared  to  the  corresponding  responses  of  the  higher  order  system 
to  determine  the  closeness  of  the  match.  These  comparisons  were 
made  in  terms  of  cost  functions  which  mathematically  define  the 
quality  of  the  match.  The  frequency  domain  and  time  domain  cost 
functions  that  were  used  to  analyze  the  quality  of  the  LOES 
matching  are  shown  in  equations  17  and  18. 

For  a  given  aircraft  configuration,  the  rise  time  of  the  ramp 
input  and  the  length  of  the  matching  interval  were  varied  until 
the  best  identification  of  the  system  (defined  by  low  cost 
functions  and  realistic  values  of  La )  had  been  achieved.  In 
general,  higher  cost  functions  were  accepted  in  order  to  keep  La 
fixed  at  the  basic  aircraft  value  (approx.  0.70  sec”*).  This 
technique  produced  reasonable  results  for  all  the  LOES  parameters 
whereas  the  opposite  technique  (i.e.  freeing  La  to  reduce  the  cost 
functions)  resulted  in  unrealistic  values  for  some  LOES 
parameters. 

The  estimation  of  flying  qualities  levels  was  based  on  the 
criteria  for  equivalent  £sp,  u)3p  and  t@  given  in  MIL-STD-1797  (3). 
The  requirements  for  time  delay  (T0)  and  short  period  damping 
(Iflp)  are  specified  as  follows: 
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Level  Allowable  Delay  Allowable  gg-  Limits 

1  0.10  (sec)  0.35-1.30 

2  0.20  (sec)  0.25-0.35  or  1.30-2.00 

3  0.25  (sec)  <  0.25  or  >  2.00 


The  requirements  for  equivalent  short  period  natural 
frequency  (uSp)  are  given  In  terms  of  the  control  anticipation 
parameter  (CAP)  defined  as  CAP  =  wap2/(n/a)  where  n/a  = 

( V/g ) ( 1/Tq2 ) *  The  requirements  on  CAP  are  listed  below. 


Level 

1 

2 

3 


Allowable  CAP  Limits  (g  *sec  1 
0.16-3.60 

0.05-0.16  or  3.60-10.00 
<  0.05  or  >  10.00 


The  requirements  shown  above  are  for  Category  C  (approach, 
landing,  and  takeoff)  flight  phases.  The  results  of  the  L0ES 
matching  including  the  L0ES  parameters,  cost  functions,  and 
predicted  flying  qualities  based  on  the  above  criteria  are  shown 
in  Table  X.  Graphical  results  are  presented  in  Appendix  G. 


Table  x 


Least 

Square a 

LOSS 

Matching  Results  (Analytical  Data) 

C0NFXG 

LEVEL 

1/102  *9  *ep 

( l/aac ) (rpa/lb ) 

"ap 

(rps) 

Te 

(sac) 

CAP 

( 1/faec) 

coat  f 

co»tt 

1-1 

1 

0.708 

0.012 

0.599 

0.922 

0.056 

0.189 

45.37 

0.005 

1-3 

2 

0.706 

0.007 

0.396 

0.769 

0.185 

0.131 

156.25 

0.062 

1-10 

3 

0.707 

0.006 

0.335 

0.718 

0.340 

0.115 

459.62 

0.117 

2-1 

1 

0.690 

0.015 

0.574 

1.751 

0.044 

0.682 

74.33 

0.008 

2-D 

1 

0.709 

0.013 

0.543 

1.609 

0.070 

0.575 

108.24 

0.009 

2-2 

2 

0.715 

0.011 

0.471 

1.518 

0.107 

0.512 

186.35 

0.019 

2-5 

3 

0.719 

0.003 

0.383 

0.894 

0.205 

0.177 

250.02 

0.015 

2-7 

2 

0.700 

0.010 

0.444 

1.461 

0.143 

0.474 

343.43 

0.027 

3-1 

1 

0.710 

0.032 

0.420 

2.711 

0.038 

1.634 

107.92 

0.026 

3-3 

2 

0.696 

0.014 

0.350 

1.909 

0.140 

0.810 

387.44 

0.062 

3-5 

2 

0.709 

0.006 

0.445 

1.266 

0.186 

0.356 

391.75 

0.021 

3-6 

2 

0.7.04 

0.022 

0.344 

2.341 

0.105 

1.218 

381.30 

0.067 

3-8 

2 

0.706 

0.018 

0.293 

2.159 

0.168 

1.036 

548.45 

0.095 
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The  cost  functions  give  a  mathematical  measure  of  the 
"goodness”  of  the  match.  As  a  general  rule  of  thumb,  MIL-STD-1797 
states  that  a  frequency  domain  mismatch  function  of  200  or  less 
will  produce  LOES  dynamics  that  are  similar  enough  to  the  high 
order  system  that  pilots,  in  general,  will  not  be  able  to  tell 
them  apart  (3). 

The  matching  results  in  Table  X  show  that  cost^  was  less  than 
200  for  only  6  of  the  13  configurations.  A  costf  greater  than 
200,  however,  does  not  mean  that  the  LOES  does  not  approximate  the 
H0S  well.  The  MIL-STD  gives  200  as  a  guideline  above  which  pilots 
may  be  able  to  tell  some  difference  between  the  H0S  and  LOES 
dynamics.  Additionally,  the  fact  that  the  H0S  does  not  match  well 
to  the  LOES  is  an  indication  in  itself  that  the  flying  qualities 
will  not  be  "conventional”.  That  is,  the  LOES  is  in  the 
conventional  short  period  form,  and  a  H0S  that  can’t  match  this 
form  is  already  suspect  of  having  poor  flying  qualities.  Table  X, 
for  example,  shows  that  the  configurations  with  cost^  greater  than 
200  all  predicted  Level  2  or  3  flying  qualities.  The  flight  test 
results,  presented  in  Table  XI,  in  general,  agreed  with  these 
predictions. 

Note  that  the  "Actual  Level”  presented  in  Table  XI  was 
generally  determined  by  averaging  the  pilot  ratings  from  only  two 
evaluations.  Even  though  the  data  scatter  between  pilot  ratings 
was  small,  this  was  a  limited  sample  and  the  following  discussions 
are  presented  with  this  consideration  in  mind. 


Table  XI 


LOES 

Predicted  Handling 

Qua  1 ities 

and 

Flight  Test  Pilot 

Ratings 

Configuration 

Predicted 

Pilot 

Actual 

Level 

Ratings 

Level 

1-1 

1 

2/4 

1 

1-3 

2 

7/7 

3 

1-10 

3 

10 

3 

2-1 

1 

3/2 

1 

2-D 

1 

4.5/3 

2 

2-2 

2 

4/2 

1 

2-5 

3 

8/10 

3 

2-7 

2 

4.5/5 

2 

3-1 

1 

2/3 

1 

3-3 

2 

3/3 

1 

3-6 

2 

6/5 

2 

3-6 

2 

5/6 

2 

3-8 

2 

3/7/4/4 

2 

Table  XI  shows  that  the  LOES  matching  technique  correctly 
predicted  the  flying  qualities  for  9  of  the  13  configurations 
(69  percent  correct).  Of  the  four  that  were  incorrect,  three  were 
on  the  "borderline”  in  some  sense.  Configuration  2-2  for  example 
was  only  7  milliseconds  above  the  Level  1  boundary  for  time  delay 
(and  thus  rated  Level  2  by  the  criteria).  The  pilot  ratings  were 
also  on  the  borderline  between  Level  1  and  2  for  configuration  2-2 
(i.e.  pilot  ratings  of  2  and  4  could  actually  be  considered  as 
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Level  1  or  Level  2,  especially  when  considering  the  “noise”  level 
of  plus  or  minus  1  pilot  rating  often  cited  as  normal  data 
scatter).  With  these  considerations  in  mind,  the  LOES  technique 
apparently  correctly  identified  a  “borderline”  aircraft  response 
for  configuration  2-2. 

Similar  results  were  seen  with  configurations  1-3  and  2-D. 
Configuration  1-3  was  15  milliseconds  below  the  Level  3  boundary 
for  equivalent  time  delay  and  received  pilot  ratings  of  7/7. 
Configuration  2-D  was  30  milliseconds  below  the  Level  2  boundary 
and  received  pilot  ratings  of  4.5/3.  If  these  three 
configurations  are  considered  as  correctly  predicting  borderline 
cases  instead  of  incorrectly  predicting  the  flying  quality  level, 
12  of  the  13  configurations  were  correctly  predicted  by  the  LOES 
technique.  This  means  92  percent  were  correctly  predicted. 

These  results  do  not  suggest  that  the  LOES  criteria  need  to 
be  changed.  The  problem  of  predicting  flying  qualities  for 
“borderline"  cases  will  always  be  there.  Therefore  it  is 
important  to  test  configurations  that  are  borderline  in  nature  to 
better  define  what  the  boundaries  should  be.  The  recommendation 
of  the  HAVE  CONTROL  test  team  was  to  conduct  more  flight  testing 
using  configurations  that  are  borderline  in  equivalent  time  delay, 
equivalent  damping,  or  CAP  to  better  identify  LOES  boundaries. 

With  regard  to  borderline  configurations,  it  should  be 
recognized  that  the  flying  qualities  of  different  aircraft  can’t 
always  be  categorized  into  three  distinct  levels  but  are  better 
described  by  a  continuous  spectrum.  The  “bottom  line”  should  be 
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to  communicate  to  the  pilot  or  the  design  engineer  the  expected 
flying  qualities  of  the  aircraft  before  the  actual  flight. 
Therefore,  borderline  configurations  should  be  identified  as  such. 
This  is  especially  important  in  cases  such  as  configuration  1-3 
which  was  predicted  to  be  Level  2  (actually  on  the  borderline 
between  Level  2  and  3)  and  rated  Level  3  by  the  pilots. 

Bandwidth  Method  Correlation 

Hoh’s  bandwidth  method  was  used  to  compare  flying  qualities 
prediction  results  with  the  LOES  method.  The  bandwidth  technique 
predicts  levels  of  handling  qualities  based  on  phase  delay  and 
system  bandwidth.  Large  phase  delays  within  the  pilot’s  bandwidth 
of  control  (0.1  to  10  rad/sec)  indicate  the  presence  of  higher 
order  lag  dynamics.  In  general,  poor  handling  qualities  are 
associated  with  low  bandwidths  and  large  phase  delays. 

Figure  46  presents  Hoh’s  proposed  Level  1,  2,  and  3 
boundaries  based  on  phase  delay  and  bandwidth.  Also  shown  are  the 
test  configurations  along  with  their  associated  actual  pilot 
ratings.  A  detailed  discussion  of  the  bandwidth  method  and  the 
graphical  results  for  each  configuration  are  shown  in  Appendix  H. 
The  bandwidth  results  are  summarized  and  compared  with  the  actual 
pilot  ratings  in  Table  XII. 
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Bandwidth  (rad/ sec) 

Figure  46.  Handling  Qualities  Prediction  Based  on  Hoh's  Bandwidth  Criteria 


Table  XII 


Bandwidth  Predicted 

Handling  Qualities 

and 

Flight  Test 

Pilot  Ratings 

Conf iguration 

Predicted 

Pilot 

Actual 

Level 

Ratings 

Level 

1-1 

2 

2/4 

1 

1-3 

3 

7/7 

3 

1-10 

3 

10 

3 

2-1 

2 

3/2 

1 

2-D 

2 

4.5/3 

2 

2-2 

2 

4/2 

1 

2-6 

3 

8/10 

3 

2-7 

2 

4.  5/5 

2 

3-1 

1 

2/3 

1 

3-3 

2 

3/3 

1 

3-5 

3 

6/5 

2 

3-6 

2 

5/6 

2 

3-8 

3 

3/7/4/4 

2 

Hoh’s  theory  correctly  predicted  the  level  of  handling 
qualities  for  7  of  the  13  configurations  (54  percent).  Two  of  the 
configurations  which  were  not  correctly  predicted  were  close  to 
the  Level  2  boundary  (configurations  3-5  and  3-8).  If  the 
boundary  was  shifted  up  slightly,  these  configurations  would  have 
been  correctly  predicted.  Additionally,  shifting  the  Level  2 
boundary  up  brings  configuration  3-3  more  into  the  heart  of  the 
Level  2-  boundary  which  makes  sense.  Configuration  3-3  was  rated 
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Level  1  by  the  pilots  but  close  to  the  Level  2  boundary  (almost 
rated  Level  3)  by  the  bandwidth  criteria. 

Shifting  up  the  Level  2  boundary  is  further  justified  by  the 
fact  that  Tp  was  higher  than  T0  for  12  of  the  13  configurations 
tested  (a  comparison  between  Tp  and  Tg  is  made  in  the  following 
section).  Since  the  Level  2  limit  for  Tg  has  been  set  at  200 
msec,  and  Tp  tends  to  be  higher  than  Tg  for  a  given  conf iguration, 
it  is  reasonable  to  set  the  Level  2  limit  for  Tp  to  200  msec  as 
well.  Based  on  this  limited  flight  test  data.  Figure  47  presents 
the  HAVE  CONTROL  test  team’s  suggested  Level  2  boundary  change  for 
Hoh’s  bandwidth  criteria.  Using  the  new  boundary,  the  bandwidth 
method  would  have  correctly  predicted  9  of  the  13  configurations 
(69  percent). 

The  purpose  of  using  the  bandwidth  method  in  this  study  was 
to  compare  the  results  to  the  Least  Squares  LOES  technique.  With 
the  limited  flight  test  data  and  associated  considerations 
presented  in  this  discussion,  the  LOES  technique  correctly 
predicted  the  handling  qualities  for  92  percent  of  the 
configurations  tested  compared  with  only  69  percent  correctly 
predicted  by  the  bandwidth  criteria.  The  Least  Squares  LOES 
computer  program  performed  very  well  compared  to  the  wel 1 -accepted 
bandwidth  method.  Application  of  the  LOES  program  was  somewhat 
more  involved  than  using  the  simple  bandwidth  method,  but  the 
results  were  more  correct.  Further  comparisons  between  the  two 
methods  including  a  comparison  between  the  time  delay  parameters 
(Tp  and  tq)  will  be  made  in  the  following  section. 
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Proposed  Boundary  Changes  for  Hoh' s  Bandwidth  Criteria 


Flight  Teat  Matching  Results 

To  verify  the  analytical  LOES  matching,  the  actual  stick 
force  inputs  and  pitch  rate  response  outputs  were  used  to  produce 
lower  order  equivalent  matches  for  each  configuration.  The 
results  of  this  matching,  including  the  LOES  parameters,  the  cost 
function  (cost^),  and  the  predicted  flying  qualities  based  on  the 
MIL-STD-1797  criteria  are  shown  in  Table  XIII.  Graphical  results 
are  shown  in  Appendix  G. 

Table  XIII 


Laast  Squares  LOES  Matching  Results  (Flight  Test  Data) 


CONFIG 

LEVEL 

I/T02  Kg 

(l/aac)(rpa/lb) 

*«P 

“ap 

(rpa) 

TD 

(sac) 

CAP 

( 1/fsac ) 

coatt 

1-1 

2 

4.023 

0.084 

1.000 

0.863 

0.100 

0.166 

0.030 

1-3 

2 

7.846 

0.010 

0.876 

2.270 

0.180 

1.145 

0.024 

1-10* 

3 

- 

- 

- 

- 

- 

- 

- 

2-1 

2 

0.744 

0.031 

1.000 

1.267 

0.120 

0.357 

0.020 

2-D 

2 

0.694 

0.041 

1.000 

0.800 

0.100 

0.142 

0.030 

2-2 

2 

0.643 

0.032 

1.000 

0.888 

0.180 

0.175 

0.030 

2-5 

3 

-2.616 

-0.041 

1.000 

1.994 

0.200 

0.884 

0.004 

2-7 

2 

3.959 

0.010 

0.384 

2.675 

0.140 

1.590 

0.051 

3-1 

1 

0.493 

0.067 

0.959 

3.032 

0.040 

2.042 

0.056 

3-3 

2 

0.694 

0.028 

0.752 

2.264 

0.120 

1.140 

0.024 

3-5 

2 

0.518 

0.026 

1.000 

0.720 

0.160 

0.115 

0.018 

3-6 

2 

0.718 

0.038 

0.675 

2.680 

0.100 

1.597 

0.019 

3-8 

2 

0.946 

0.027 

0.508 

2.493 

0.160 

1.381 

0.031 

*  Configuration  1-10  was  not  flown  in  landing  approach  due  to 
uncontrollability  during  level  off  task 


The  results  presented  in  Table  XIII  demonstrate  some  of  the 
problems  that  were  encountered  during  the  flight  test  matching. 

One  of  the  problems  was  that  the  ramp  input  rise  time  and/or  the 
length  of  the  match  time  could  not  be  varied  as  easily  as  was  done 
in  the  analytical  matching.  During  the  analytical  matching,  these 
inputs  could  be  varied  continuously  until  the  best  match  was 
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achieved.  This  is  how  the  parameter  La  was  fixed  at  the  basic 
aircraft  value  and  how  the  cost  functions  were  minimized. 

Because  of  the  difficulties  involved  in  precisely  varying  the 
input  rise  time  and  the  length  of  available  match  time  some 
configurations  did  not  match  well  using  the  flight  test  data  and 
realistic  values  of  La  were  not  achieved  in  every  case.  A 
realistic  value  of  La  is  important  because  this  generally  effects 
the  accuracy  of  the  other  equivalent  parameters.  Table  XIII  shows 
that  four  of  the  configurations  (1-1,  1-3,  2-5,  and  2-7)  had 
unrealistic  values  of  1/Tq2  (and  therefore  La ) . 

In  the  cases  where  a  poor  match  was  achieved,  generally  the 
rise  time  of  the  ramp  input  was  ’’close  enough”  based  on  analytical 
predictions  but  not  enough  data  was  collected.  That  is,  the  length 
of  available  pitch  rate  response  data  was  insufficient  to  identify 
the  system  dynamics.  Depending  on  the  configuration,  anywhere 
from  3  to  10  seconds  of  data  were  generally  required  to  make  a 
good  match.  The  problem  was  that  to  get  a  good  response  from  the 
aircraft,  a  fairly  large  stick  force  input  (about  5  lbs)  was 
required.  This  much  input,  however,  caused  a  large  pitch  rate  up 
that  had  to  be  recovered  before  the  nose  was  too  high.  So,  there 
was  a  trade  off  between  getting  a  good  response  from  the  aircraft 
and  getting  enough  data  to  make  the  system  identification. 

Additional  problems  included  noise  and/or  turbulence  effects. 
In  some  cases,  the  rise  time  was  correct  and  there  was  enough  data 
to  make  the  match,  but  noise  in  the  data  or  turbulence  made  the 
data  unusable.  These  problems  were  minimized  however  because  the 
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flights  were  conducted  early  in  the  morning  in  smooth  air. 


Despite  the  problems  mentioned,  the  overall  results  of  the 
flight  test  data  matching  were  encouraging.  Table  XIV  compares 
the  predicted  handling  qualities  with  the  actual  pilot  ratings. 
These  results  show  that  8  of  the  13  conf igurations  were  correctly 
predicted  by  the  flight  test  matching  (62  percent  correct).  Of 
the  5  that  were  incorrect,  4  were  on  the  ’’border  1  ine”  in  time 
delay  as  was  discussed  previously  with  regard  to  the  analytical 
matching  results. 


Table  XIV 


Flight  Test  Matching 
and  Flight  Test  Pilot 

Results 

Ratings 

Conf iguration  Predicted 

Level 

Pilot 

Ratings 

Actual 

Level 

1-1 

2 

2/4 

1 

1-3 

2 

7/7 

3 

1-10 

3 

10 

3 

2-1 

2 

3/2 

1 

2-D 

2 

4.5/3 

2 

2-2 

2 

4/2 

1 

2-5 

3 

8/10 

3 

2-7 

2 

4.  5/5 

2 

3-1 

1 

2/3 

1 

3-3 

2 

3/3 

1 

3-5 

2 

6/5 

2 

3-6 

2 

5/6 

2 

3-3 

2 

3Z7/4/4 

2 
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The  overall  results  of  the  flight  test  matching  show  that 


good  results  can  be  obtained  when  large  enough  inputs  are  flown  in 
smooth  air  and  an  adequate  amount  of  data  is  collected  before  the 
pitch  up  is  terminated.  In  addition,  it  was  found  that  the  pilots 
could  fly  the  ramp  and  step  inputs  precisely  enough  by  hand  to  be 
used  in  the  LOES  program.  In  fact,  the  hand-ramp  inputs  were  used 
in  several  of  the  conf iguration  matches  tabulated  in  Table  XIII 
because  the  inputs  worked  as  well  or  better  than  the  auto-ramp 
inputs.  This  is  an  important  finding  because  it  means  that  the 
LOES  system  identification  program  can  be  used  on  virtually  an: 
instrumented  test  aircraft  since  an  auto-input  capability  is  not 
required.  All  that  is  required  is  a  DAS  that  can  measure  stick 
force  input  and  pitch  rate  output. 

Using  the  flight  test  data  an  additional  time  delay 
measurement  was  made  which  will  be  referred  to  as  "TD”  to  be 
distinguished  from  phase  delay  (Tp)  and  equivalent  time  delay 
( Te  or  Tq).  TD  was  found  by  measuring  the  delay  between  the  first 
onset  of  stick  force  and  the  resulting  response  in  pitch  rate. 
These  measurements  are  shown  graphically  in  Appendix  G. 

Figure  48  gives  a  comparison  between  Te,  Tp,  and  TD  which 
shows  them  to  be  numerically  similar.  As  was  discussed  earlier, 

Tp  tends  to  be  higher  than  Te  for  a  given  configuration.  Despite 
this  fact,  the  bandwidth  criteria  are  more  restrictive  on 
allowable  Tp  than  the  LOES  criteria  are  on  Te.  This  was  one  of 
the  reasons  for  the  recommendation  to  change  the  bandwidth 
criteria . 
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1-1  1-3  1-10  2-1  2-D  2-2  2-5  2-7  3-1  3-3  3-5  3-6  3*9 

CONFIGURATION 

Figure  48.  Comparison  of  Time  Delay  Parameters 

Figure  48  provides  a  good  illustration  of  the  importance  of 
time  delay  as  a  flying  qualities  parameter.  Configuration  1-10, 
for  example,  had  significantly  higher  values  of  time  delay  than 
the  other  configurations  and  was  also  the  worst  flying 
configuration  (1-10  was  not  evaluated  in  the  landing  task  because 
it  was  uncontrollable  at  altitude).  In  general,  the  handling 
quality  predictions  were  largely  determined  by  time  delay  (rather 
than  t3  p  or  u)sp ) .  Configurations  with  lower  time  delay  were 
generally  predicted  to  have  good  handling  qualities  and 
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configurations  with  higher  time  delay  were  predicted  to  have  poor 
handling  qualities.  The  flight  test  results,  in  general,  agreed 
with  these  predictions. 

The  overall  results  of  the  analytical  and  flight  test 
matching  have  shown  that  the  Least  Squares  program  is  well  suited 
to  perform  lower  order  equivalent  systems  matching  and  predict 
handling  qualities.  The  results  of  the  LOES  predictions  compared 
favorably  with  the  Bandwidth  results  and  with  the  actual  pilot 
ratings.  These  results  are  summarized  in  Table  XV. 


Table  XV 


Flight  Test  Handling  Quality  Prediction  Results 
(Percent  Correctly  Predicted) 

Actual  MIL-STD 
Boundaries 

Border  line 
Cases  Correct 

Bandwidth 
( Pre-Flight ) 

54 

69 

LS  LOES 
(Pre-Flight) 

69 

92 

LS  LOES 
(Flight  Data) 

62 

92 

VI.  Conclusions  and  Recommendations 


A  program  has  been  developed  for  matching  the  time  response 
of  high-order  aircraft  systems  with  lower-order  equivalent 
systems.  The  lower-order  equivalent  systems  are  extracted  in  the 
form  of  the  short  period  pitch  rate  transfer  function  to  take 
advantage  of  the  “classical’*  aircraft  data  base  in  MIL-STD-1797 . 
The  attractive  aspects  of  the  Least  Squares  LOES  program  are  its 
simplicity  and  its  adaptability  to  flight  test  parameter 
identification  requirements. 

Comparisons  between  the  Least  Squares  program  and  LONFIT  have 
shown  that  both  programs  give  numerically  similar  LOES  parameters 
for  a  given  high-order  system.  The  LOES  parameters  derived  from 
both  programs  were  found  to  approach  the  corresponding  basic 
aircraft  values  as  the  control  system  dynamics  were  made  “faster” 
compared  to  the  aircraft  dynamics.  That  is,  as  the  control  system 
dynamics  became  more  "transparent”,  the  equivalent  parameters 
predicted  that  the  aircraft  would  have  a  classical  low-order 
response.  This  is  consistent  with  LOES  theory. 

During  the  flight  test  program  the  Least  Squares  technique 
correctly  predicted  the  flying  qualities  levels  of  92  percent  of 
the  aircraft  configurations  tested  compared  with  only  69  percent 
correctly  predicted  by  Hoh’s  Bandwidth  method.  Additionally,  the 
Bandwidth  method’s  time  delay  parameter  ( Tp )  was  found  to  be 
higher  than  the  LOES  time  delay  parameter  (ts)  for  12  of  the  13 
configurations  tested.  Despite  this  fact,  the  MIL-STD-1797  Level 
2  boundary  for  Tp  i3  lower  than  the  Level  2  boundary  for  Te. 
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These  findings  resulted  in  a  recommendation  by  the  test  team  to 
change  the  Level  2  Bandwidth  boundary. 

It  was  found  that  the  short  period  numerator  term,  La,  could 
be  fixed  at  thd  basic  aircraft  value  by  varying  the  rise  time  of 
the  ramp  input  and/or  the  length  of  the  overall  match  time. 

Fixing  La  was  the  preferred  method  of  matching  since  this  gave 
more  realistic  values  for  all  LOES  parameters  and  also  gave  the 
most  correct  flying  qualities  predictions. 

Fixing  La  during  the  analytical  matching  was  easy  since  the 
rise  time  and  length  of  match  time  could  be  varied  as  needed. 

These  inputs  could  not  be  varied  as  easily  with  the  actual  flight 
test  data  however.  Generally  the  rise  time  of  the  input  could  be 
controlled  adequately  but  the  problem  was  in  getting  enough  data. 
There  was  a  trade  off  between  getting  a  good  response  out  of  the 
airplane  (which  required  a  fairly  large  stick  force  input)  and 
getting  enough  data  before  the  resulting  rapid  pitch  up  had  to  be 
recovered . 

An  important  finding  during  the  flight  test  matching  was  that 
the  pilots  could  hand  fly  the  required  stick  force  ramp  inputs. 

One  of  the  reasons  the  NT-33A  aircraft  was  used  for  the  test 
program  was  to  take  advantage  of  the  auto-input  capability.  It 
was  thought  that  the  pilots  could  not  fly  the  ramp  inputs  with  the 
precision  required  by  the  Least  Squares  program.  In  actual  flight 
test  applications  however,  it  was  found  that  a  smooth  ramp  input 
of  approximately  1  second  was  adequate.  The  pilots  were 
instructed  that  the  smoothness  of  the  input  was  more  important 
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than  the  actual  rise  time.  The  results  showed  that  in  many  cases 
the  hand-flown  inputs  worked  as  well  or  better  than  the  auto¬ 
inputs.  This  is  an  important  finding  because  it  means  the  Least 
Squares  program  can  be  used  with  any  DAS  instrumented  aircraft. 

Despite  the  problems  already  mentioned  with  the  actual  flight 
test  matching,  9  of  the  13  configurations  were  matched  well  (with 
La  close  to  the  basic  aircraft  value)  and  12  out  of  13  were 
correctly  predicted  (92  percent).  These  results  are  encouraging 
and  the  recommendation  is  to  conduct  more  flight  testing  with  the 
Least  Squares  LOES  program  to  determine  how  the  problems  can  be 
solved.  The  program  has  been  run  at  the  Air  Force  Flight  Test 
Center  (AFFTC)  Simulation  Facility  on  Matrixx  and  can  be  applied 
to  almost  any  flight  test  mission  using  only  about  5-10  minutes  of 
actual  flight  time.  The  test  aircraft  must  have  a  data 
acquisition  system  capable  of  recording  stick  force  input  and 
pitch  rate  output  at  10  samples  per  second. 

Parameter  identification  techniques  will  soon  be  a  part  of 
the  curriculum  at  the  USAF  Test  Pilot  School  and  this  may  be  a 
good  program  to  use.  Potential  starting  points  for  an  AFIT/TPS 
thesis  project  could  be  refining  the  overall  technique  to  work 
better  using  flight  test  data,  modifying  the  program  to 
simultaneously  match  pitch  rate  and  normal  acceleration,  modifying 
the  program  to  identify  1 atera 1 -directiona 1  equivalent  parameters, 
and  conducting  comparisons  with  MMLE  or  other  LOES  programs 
regarding  pre-flight  flying  qualities  predictions. 
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The  overall  test  results  of  the  Least  Squares  LOES  program 
have  shown  that  it  works  well  in  comparisons  with  LONFIT,  Hoh’s 
Bandwidth  method,  and  actual  pilot  ratings.  The  recommendation  is 
to  use  the  Least  Squares  LOES  program  in  applications  requiring 
flight  test  parameter  identification  and  handling  qualities 
predictions. 
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APPENDIX  A.  FREQUENCY  AND  TIME  RESPONSE  MATCHING  RESULTS 
This  appendix  contains  the  equivalent  parameters,  CAP,  and  the 
cost  functions  obtained  from  the  freqency  domain  matching  (LONFIT) 
and  the  time  domain  matching  (LEAST  SQUARES)  techniques.  The  data  is 
presented  for  the  LONFIT  results,  the  LEAST  SQUARES  results,  and 
finally  a  comparison  between  the  two  is  given  in  tabular  format. 
Graphical  comparisons  of  this  data  are  presented  in  the  discussion 
(chapters  III  and  IV). 


TABLE  A1 


LONFIT  FREQUENCY  RESPONSE  MATCHING  RESULTS 
FLIGHT  CONDITION  1* 


“fs 

La 

La 

K0 

^sp 

U)__ 

sp 

Te 

CAP 

costf 

cost^- 

(rps) 

- 

(s_1) 

( rps/ lb ) 

- 

( rps) 

(  sec ) 

l/(g*s) 

- 

- 

6.0 
6.  0 

Fix 

Free 

0.  428 
0.909 

-0.092 

-0.073 

0.245 
0.  147 

2.270 

2.517 

0.287 

0.269 

0.569 

0.329 

264.62 
179. 08 

1.017 
1. 652 

8.0 

8.0 

Fix 

Free 

0.428 

0.750 

-0. 106 
-0.090 

0.236 
0.  163 

2.404 
2.  572 

0.  252 
0.239 

0.638 

0.417 

178.29 

123.79 

0.851 

1.393 

10.0 

10.0 

Fix 

Free 

0.428 

0.681 

-0. 116 
-0. 101 

0.235 
0.  174 

2.481 

2.612 

0.226 

0.215 

0.  680 
0.473 

135. 16 
95.21 

0.745 

1.263 

12.0 

12.0 

Fix 

Free 

0.428 

0.644 

-0. 122 
-0. 108 

0.  235 
0.  181 

2.529 

2.640 

0.  206 
0.  196 

0.706 
0.  511 

111.56 
79.  50 

0.688 
1.  193 

18.5 

18.5 

Fix 

Free 

0.  428 
0.  595 

-0. 133 
-0. 120 

0.238 
0.  193 

2.  601 
2.686 

0.164 
0.  156 

0.747 
0.  572 

81.80 

59.95 

0.626 
1. 102 

31.0 

31.0 

Fix 

Free 

0.  428 
0.572 

-0. 139 
-0. 127 

0.240 

0.201 

2.  640 
2.712 

0.  127 
0.  120 

0.769 
0.  607 

69.36 
52.  12 

0.612 

1.071 

^Flight  condition  1  is  at  M  =  0.70,  Alt  =  35,000  feet,  VTo  =  681  fps, 
and  q  =  171  psf.  The  basic  aircraft  short  period  dynamics  for  the 
A-4D  at  this  flight  condition  are  La  =  0.428  sec-*,  £sp  =  0.225,  and 
wap  =  2.77  rad/s. 
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TABLE  A2 


LONFIT  FREQUENCY  RESPONSE  MATCHING  RESULTS 

FLIGHT  CONDITION  2* 

ufs 

La 

La 

Ke 

£sp 

“sp 

Te 

CAP 

COStf 

COStfc 

( rps ) 

- 

(  s-1 ) 

( rps/ lb ) 

- 

(rps) 

(  sec ) 

l/(g*s) 

- 

- 

6.0 

6.0 

Fix 

Free 

2.080 
8.  197 

-0.059 
-0. 027 

0.720 

0.369 

4.524 
5.  887 

0.220 
0.  168 

0.333 
0.  143 

57.908 
33. 290 

0.069 
0.  125 

8.0 

8.0 

Fix 

Free 

2.080 
5.  Ill 

-0.078 

-0.044 

0.654 
0.  354 

5.  192 
5.949 

0.  201 
0.  163 

0.  439 
0.235 

50. 142 
32.878 

0.073 
0.  133 

10.0 

10.0 

Fix 

Free 

2.080 

3.989 

-0.091 

-0.059 

0.613 

0.368 

5.631 

6.065 

0.  185 
0.  155 

0.  516 
0.312 

45.050 

31.650 

0.077 
0.  130 

12.0 

12.0 

Fix 

Free 

2.  080 
3.546 

-0. 102 
-0.069 

0.  587 
0.376 

5.  930 
6.204 

0.  172 
0.  146 

0.  573 
0.368 

41.742 

30.648 

0.079 
0.  129 

18.5 

18.5 

Fix 

Free 

2.080 

3.059 

-0. 120 
-0.088 

0.549 
0.  386 

6.433 

6.508 

0.  141 
0.  121 

0.674 

0.469 

36. 734 
29.051 

0.085 
0.  128 

31.0 

31.0 

Fix 

Free 

2.  080 
2.  853 

-0. 133 
-0. 101 

0.529 

0.391 

6.739 

6.726 

0.  112 
0.094 

0.740 

0.537 

34. 157 
28.263 

0.088 
0.  128 

*Flight  condition  2  is  at  M  =  0.85,  Alt  =  0  feet,  VjQ  =  950  fps,  and 
q  =  945  psf.  The  basic  aircraft  short  period  dynamics  for  the  A-4D 
at  this  flight  condition  are  La  =  2.08  sec-1,  ?SD  =  0.436,  and  uiSD  = 
7.05  rad/s. 


TABLE  A3 

LEAST  SQUARES  TIME  RESPONSE  MATCHING  RESULTS 
FLIGHT  CONDITION  1 


wfs 

La 

La 

^sp 

sp 

Te 

CAP 

COStf 

COStf. 

(rps) 

- 

(s-1 ) 

(rps/ lb) 

- 

( rps ) 

( sec ) 

l/(g*s) 

- 

- 

6.0 

Free 

0.649 

-0.053 

0.  120 

2.  196 

0.  261 

0.351 

380.76 

1. 835 

8.0 

Free 

0.568 

-0.063 

0.133 

2.  265 

0.  229 

0.  427 

317.82 

1.648 

10.0 

Free 

0.  523 

-0.071 

0.  143 

2.310 

0.  204 

0.  482 

282.33 

1.456 

12.0 

Free 

0.495 

-0.076 

0.  150 

2.342 

0.  185 

0.524 

259.24 

1.296 

18.5 

Free 

0.448 

-0.087 

0.  167 

2.405 

0.  146 

0.611 

218.04 

0.965 

31.0 

Free 

0.  413 

-0. 097 

0.  182 

2.  459 

0.  113 

0.691 

185. 80 

0.691 
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TABLE  A4 


LEAST 

SQUARES  TIME  RESPONSE  MATCHING  RESULTS 
FLIGHT  CONDITION  2 

wfs 

La 

La 

K0 

^sp 

wsp 

Te 

CAP 

COSt£ 

costt 

(  rps  ) 

- 

<s-l) 

( rps/ lb ) 

- 

( rps ) 

(  sec ) 

l/(g*s) 

- 

- 

6.0 

Free 

8.  416 

-0.009 

0.258 

3.  890 

0.074 

0.061 

433. 24 

0.026 

8.  0 

Free 

5.674 

-0.016 

0.243 

4.  133 

0.072 

0.  102 

384.44 

0.033 

10.0 

Free 

4.504 

-0. 021 

0.  245 

4.294 

0.066 

0.  139 

342.09 

0.038 

12.0 

Free 

3.874 

-0. 026 

0.253 

4.411 

0.060 

0.  170 

309.61 

0.041 

18.5 

Free 

2.  987 

-0.037 

0.283 

4.644 

0.044 

0.245 

245.06 

0.043 

31.0 

Free 

2.  451 

-0.050 

0.326 

4.854 

0.  029 

0.326 

190.27 

0.041 

TABLE  A5 


LONFIT  AND 

LEAST 

SQUARES 

COMPARISON 

FLIGHT  CONDITION 

1* 

wfs 

Type 

of 

La 

Ke 

?sp 

usp 

Te 

CAP 

COStjf 

costt 

(rps) 

(s_1) 

(rps/lb) 

- 

(rps) 

(  sec ) 

l/(g*s) 

- 

- 

6.0 

Lon 

0.909 

-0.073 

0.  147 

2.517 

0 . 269 

0.329 

179.08 

1. 652 

6.0 

LS 

0.649 

-0. 053 

0.  120 

2.  196 

0.261 

0.  351 

380. 76 

1.835 

8.0 

Lon 

0.750 

-0.090 

0. 163 

2.  572 

0.239 

0.417 

123.79 

1.393 

8.0 

LS 

0.568 

-0.063 

0.133 

2.265 

0.  229 

0.427 

317.82 

1.648 

10.  0 

Lon 

0.681 

-0. 101 

0.  174 

2.612 

0.  215 

0.  473 

95.21 

1.263 

10.0 

LS 

0.523 

-0. 071 

0.  143 

2.310 

0.204 

0.482 

282.33 

1.456 

12.0 

Lon 

0.644 

-0. 108 

0.  181 

2.640 

0.  196 

0.511 

79. 

1. 193 

12.0 

LS 

0.495 

-0.076 

0.  150 

2.  342 

0.  185 

0.  524 

259.24 

1.296 

18.5 

Lon 

0.  595 

-0. 120 

0.193 

2.686 

0.  156 

0.  572 

59.95 

1.  102 

18.  5 

LS 

0.  448 

-0. 087 

0.  167 

2.  405 

0.  146 

0.611 

218.04 

0.965 

31.0 

Lon 

0.572 

-0. 127 

0.201 

2.712 

0.  120 

0.607 

52.  12 

1.071 

31.0 

LS 

0.413 

-0.097 

0.  182 

2.  459 

0.  113 

0.  691 

185. 80 

0.691 

*La  free  for  all  configurations 
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TABLE  A6 


LONFIT  AND  LEAST  SQUARES  COMPARISON 
FLIGHT  CONDITION  2* 


wfs 

Type 

of 

La 

*0 

?sp 

wsp 

Te 

CAP 

COSt£ 

cost^ 

(rps) 

— 

(  s_1 ) 

( rps/ lb ) 

— 

(rps) 

( sec ) 

l/(g*S) 

— 

- 

6.0 

Lon 

8.  197 

-0.027 

0.359 

5.887 

0.  168 

0.  143 

33.290 

0.  125 

6.0 

LS 

8.416 

-0.009 

0.  258 

3.890 

0.074 

0.061 

433. 24 

0.026 

8.0 

Lon 

5.  Ill 

-0. 044 

0.354 

5.949 

0.  163 

0.235 

32. 878 

0.  133 

8.0 

LS 

5.674 

-0. 016 

0.243 

4.  133 

0.  072 

0.  102 

384.44 

0.033 

10.  0 

Lon 

3.  989 

-0. 059 

0.  368 

6.  065 

0.  155 

0.  312 

31. 650 

0.  130 

10.0 

LS 

4.504 

-0.021 

0.  245 

4.  294 

0.066 

0.139 

342.09 

0.038 

12.0 

Lon 

3.546 

-0.069 

0.376 

6.204 

0.  146 

0.368 

30.648 

0.  129 

12.0 

LS 

3.  874 

-0.026 

0.253 

4.411 

0.060 

0.  170 

309.61 

0.041 

18.5 

Lon 

3.059 

-0.088 

0.  386 

6.508 

0.  121 

0.469 

29.051 

0.  128 

18.5 

LS 

2.987 

-0.037 

0.283 

4.644 

0.044 

0.245 

245.06 

0.043 

31.0 

Lon 

2.853 

-0. 101 

0.391 

6.726 

0.094 

0.537 

28.263 

0.  128 

31.0 

LS 

2.451 

-0.050 

0.326 

4.854 

0.029 

0.326 

190.27 

0.041 

*La  free  for  all  configurations 
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APPENDIX  B.  LEAST  SQUARES  MATRIX*,  PROGRAMS 

The  computer  programs  used  for  the  least  squares  matching, 
the  extraction  of  continuous  state  space  matrices  from  the 
discrete  matrices,  and  the  calculation  of  the  LOES  parameters  and 
cost  functions  were  all  written  for  MATRIXx  (31).  The  MATRIXx 
programs  were  written  as  main  programs  with  subroutines.  The 
subroutines  are  executable  files  and  are  preceded  by  the  command 
"exec"  (meaning  execute).  The  first  program,  called  COST. FREQ, 
was  used  to  calculate  costf ,  cost^,  and  CAP  given  the  results  from 
the  LONFIT  (frequency  domain)  matching. 

To  run  COST.  FREQ  on  MATRIXx  an  input  variable  ’’in’*  must  exist 
in  the  MATRIXx  memory  where  in  =  <lalpha  gain  zeta  omega  tau>. 

The  vectors  "numhos”  and  ’’denhos”  must  also  be  defined  where 
numhos  is  the  numerator  of  the  HOS  transfer  function  and  denhos  is 
the  denominator  of  the  HOS  transfer  function.  The  program  uses 
the  input  vector,  which  contains  the  low-order  equivalent 
parameters  from  the  frequency  domain  matching,  and  forms  the  LOES 
transfer  function.  The  LOES  and  HOS  transfer  functions  are  then 
used  to  calculate  cost£  and  cost^.  The  LOES  parameters  are  also 
used  to  calculate  CAP,  the  control  anticipation  parameter  (defined 
in  chapter  I). 

COST. FREQ 

exec( ’ costf .freq’  ) 
exec(  ’ costt . freq ’  ) 
v=681; 

cap=(32. 174*in(4)**2)/(v*in(l ) ) ; 
out=<in  cap  costf  costt> 
return 
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'costf.freq'  -  (subroutine  used  to  calculate  cost£ ) 

num=in(2)*<l  in(l)>; 

den=<l  2* in( 3 ) * in( 4 )  in(4)**2>; 

<o , dbhos , phhos>  =  bode( numhos, denhos, .1,10,21); 
<o,dblos,phlos>=bode(num,den, .1,10,21); 
if  abs( phhos( 1 )-phlos( 1 ) )  >  180,... 

for  i=l : 21 , phi os( i )=phl os( i )-360; end 
for  isl : 21 , phlos( i ) =phlos( i )-57. 29578*o( i )*in( 5 ) ; 
cost=0; 

for  i  =  l s  21 , . . . 

cost=cost+(dbhos( i )-dblos( i  )  )**2  +  0. 01745* ( phhos( i )- 
phlos( i ) ) **2; 
costf = ( 20/21 ) *cost 
return 


’costt.freq’  -  (subroutine  used  to  calculate  costt) 

for  i=l : 100, u5( i )=5; end 
< t , yhos>= lsim( numhos, denhos, u5, . 1 ) ; 

< t , y los>  = lsim( num, den, u5, . 1 ) ; 

costt=0; 

for  i=l s 100, . . . 

costt=costt  +  ( 57. 29578*yhos(  i )  -  57. 29578*y lost i ) ) **2; end 

costt =costt/ 100 

return 


The  next  main  program  is  called  LEAST. SQUARES  and  was  used  to 
perform  the  time  domain  matching  for  comparison  with  L0NFIT.  To 
run  LEAST. SQUARES  on  MATKIXx  the  variable  ”shos”  must  be  defined 
and  exist  in  the  MATRIXx  memory.  "shos "  is  the  system  matrix  of 
the  H0S  defined  as: 


shos 


ahos 

bhos 

CH0S 

dhos 

where  A^os*  ®H0S*  ^HOS’  anc*  ^HOS  represent  the  state  space  form  of 
the  high-order  system.  LEAST. SQUARES  forms  the  input  vector  u 
depending  on  the  desired  rise  time  (various  rise  times  were  used 
but  the  version  shown  below  is  for  a  rise  time  of  1.1  sec).  The 
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input  u  and  corresponding  output  y  (found  from  the  HOS  time 


response  to  u )  are  used  to  perform  the  least  squares  matching. 

Note  that  within  the  subroutine  ’ls.trll’  is  another  subroutine 
’extract. ab ’  which  is  used  to  extract  the  continuous  A  and  B 
matrices  from  the  discrete  F  and  S  forms.  The  extraction  programs 
are  shown  last  (the  Sinha  and  Lastman  technique  and  the  variation 
of  this  technique  which  uses  the  power  series  expansion  of  e^T  ag 
explained  -n  chapter  IV). 

LEAST. SQUARES 


exec( 'Is. trll ’ ) 
exec( ’ roots ’ ) 
exec( ’ tau ’ ) 
exec( ’costf.time’  ) 
exec( ’costt.time’  ) 
v=681 ; 

cap*  (  32.  174* omega*:* 2 )  /  (  v*  lalpha ) ; 

out=<lalpha  gain  zeta  omega  tau  cap  costf  costt> 
return 


’ls.trll’  -  (performs  the  least  squares  matching  using  the 
discrete  input  u  and  corresponding  output  y; 
”trll”  refers  to  the  rise  time  of  1.1  sec) 


t  =  <. Is  . Is 10>  ’  ; 

U ( 1 ) =0 ; 

for  i=l: 11, u(i+l)=u(i )+0. 4545455; end 

for  i=13: 100, u(i)=5.0; end 

ns=7; 

<t, y>  =  lsim(shos, ns, u, . 1 ) ; 

for  i  =  3: 100, x( i-2, ;  )  =  <y( i-1 )  y(i-2)  u(i-l)  u(i-2)>; 

for  i=3: 100, capy ( i-2 , 1 )  =y( i , 1 ) ; 

theta=x  capy 

yout=x*theta; 

dnum2=<theta(3)  theta(4)> 

dden2=< 1  -theta(l)  -theta(2)> 

ns=2; 

<sd, ns> =sf orm( dnum2 , dden2) ; 
c  =  sd(3, :  ) ; 
n  =  10; 

exec( ’extract. ab; 2 ’ ) ; 
b»rea 1(b);- 
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sc=  <a , b; c> ; 

<  nun , den  >  =  tf orm( sc, ns) 

<o, dls, pls>=bode( num , den, .1,10,21); 
for  i=l : 100, u5( i )=5; end 
<t,yls>=lsim(num,den,u5, .1); 
return 


’roots’  -  (the  equivalent  parameters  lalpha,  zeta,  omega,  and 

gain  are  calculated  from  ”num"  and  "den”  -  the  output 
of  the  ls.trll  program) 

rd  =  roots( den  ) ; 
sigma=abs( real ( rd( 1  )  )  ) ; 
omega  =  abs( rd( 1 )  ) ; 
zeta=sigma/omega ; 
rn= roots ( num); 
lalpha=-rn; 
gain=num( 1 ) ; 
return 


’tau’  -  (the  equivalent  time  delay  tau  is  calculated  by  iteration 
in  the  frequency  domain) 

ns  =  7 ; 

<  o , d , p>  =  bode( shos, ns,  ,1,10,21); 

<o, dl , pi >=bode( num, den, .1,10,21); 
if  abs(p( 1 )-pl( 1 ) )  >  180,... 

for  i=l : 21 , pi ( i ) =pl ( i )-360; end 
pplot=<p  pi > ; 
plot(o,pplot,  ’ logx ’  ) 
tau=nt; 
dif f  =  1.0; 

whi le  dif f  >  0, . . . 
sspl=0; . . . 
ssp2=0; . . . 
plos=pl ; . . . 

for  1-1:21, pi os( i)=plos(i)-57.29578*o(i)*tau; end .... 
for  i  =  l : 21 , sspl  =  sspl  +  0. 01745* ( p ( i )-plos( i )) **2;  end  ,..  . 
tau=tau  +  0.001;... 
plos=pl ; . . . 

for  i  =  1 : 21 , p lost i )  =  p los( i ) -57 . 29578*o( i ) *tau ;  end  ,  .  .  . 

for  i = 1 : 21 , ssp2=ssp2  +  0 . 01745* ( p( i )-pl os( i ) ) * *2; end , . . . 

diff=sspl-ssp2; 

pplot2  =  <p  p 1 os > ; 

plot(o,pplot2,  ’logx’  ) 

tau 

return 
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’costf • time '  -  (calculates  costf) 
ns-7 ; 

<o, dbhos, phhos>  =bode( shos, ns ,.1,10,21); 

<o, dblos , phlos>=bode( num, den, .1,10,21); 
if  abs( phhos( 1 )-phlos( 1 ) )  >  180,... 

for  i=l:21,phlos(i)=phlos(i )-360; end 
for  i=l:21,phlos(i)=phlos(i )-57. 29578*0 ( i ) *tau ; end 
cost=0; 

for  i=l : 21 , . . . 

cost=cost+( dbhos( i )-dblos( i ) )**2  ■+■  0. 01745* ( phhos( i )- 
phlos( i ) ) **2; 
costf = ( 20/21 ) *cost 
return 


'costt. time ’  -  (calculates  costt) 

for  i =1 : 100 , u5( i , 1 ) =5 ; end 
ns  =  7; 

<  t , yhos  >  =  1 sim( shos , ns , u5 , . 1  )  ; 

< t , y los>= lsim( num , den , u5 , .1); 

costt =0 ; 

for  i=l : 100, . . . 

costt=costt  +  ( 57 . 29578*yhos( i )  -  57 . 29578*y 1 os ( i ) ) * *2 ; end 

costt =costt/ 100 ; 

return 


’ extract . ab; 1 ’  -  (performs  Sinha  and  Lastman  extraction  technique 

of  the  continuous  A  and  B  matrices  from  the 
discrete  F  and  G  forms) 

T=0. 1; 

<  F , G , H , I >= SPLIT ( SD , NS ) ; 

AT=0.5*(F-INV(F)); 

FOR  1=1:3,... 

FSTAR=EXP( AT) ; . . . 

AT= AT+INV ( F ) * ( F-FSTAR ) ; END 
A=AT/T 

Q=INV( A)*(EXP( AT) -EYE) ; 

B=INV(Q)*G 

RETURN 


’ extract . ab ; 2  *  -  (performs  Sinha  and  Lastman  extraction  for  the 

A  matrix  and  the  power  series  extraction  for 
the  B  matrix) 


T =0 . 1; 

<  F , G , H , I >  =  SPLIT( SD , NS  ) ; 
AT=0 . 5* ( F-INV ( F ) ) ; 

FOR  1=1:3,... 
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FSTAR  =  EXP( AT ) ;  .  .  . 

AT=AT+INV( F ) * ( F-FSTAR ) ; END 
A=AT/T; 

FACT=1 . 0; 

8=0NES( 2 )-0NES( 2 ) ; 

FOR  J=0: N, . . . 

L  *  J  j  •  •  • 

F ACT = FACT* ( L+l ) ; .  . . 

K=( ( A**L)*(T**(L+1 ) ) )/FACT;  .  .  . 

Q=Q+K; END 

FACT 

K 

A 

B=INV(Q)*G 

RETURN 


The  next  program  was  used  to  conduct  the  analytical  matching 
prior  to  the  flight  test  matching.  The  results  were  used  to 
predict  flying  qualities  levels  based  on  MIL-STD-1797  guidance  for 
equivalent  LOES  parameters.  The  version  presented  here  is  called 
LS.SOF  (Least  Squares  Second  Order  Filter).  The  program  builds 
the  basic  aircraft  transfer  function  for  configuration  3  (in  this 
case)  and  then  adds  the  high-order  dynamics  by  adding  a  second 
order  pre-filter  (pre-filter  dynamics  are  for  configuration  3-8  in 
this  case). 

After  the  HOS  transfer  function  is  completed  ("numhos"  and 
’’denhos"  exist  in  the  MATRIXx  memory)  the  program  accomplishes  the 
least  squares  matching,  extracts  the  LOES  parameters,  and  computes 
the  cost  functions  similar  to  the  LEAST . SQUARES  program  although 
there  are  some  differences.  The  most  important  difference  is  that 
this  version  is  interactive,  allowing  the  user  to  specify  rise 
time,  match  time,  and  initial  estimates  for  time  delay.  This 
interaction  speeds  up  the  overall  time  required  to  achieve  the 
best  match  since  several  combinations  of  rise  time  and  match  time 
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must  be  tried.  The  ’’best  match’’  was  defined  by  La  fixed  at  the 
aircraft  value  while  simultaneously  achieving  the  lowest  possible 
cost  functions. 


LS. SOF 

inquire  tr  ' ENTER  RISE  TIME  IN  SEC:' 

inquire  tm  'ENTER  MATCH  TIME  IN  SEC:' 

inquire  nt  'ENTER  INITIAL  TIME  DELAY  IN  SEC:' 

exec  ('  conf  ig.  3' ) 

exec ( ' nd. sof ' ) 

exec (' Is .match' ) 

exec ( ' roots . ' ) 

exec ( ' tau . ' ) 

exec (' costf .time' ) 

exec ('  costt .time' ) 

v=205; 

cap= (omega* *2) /4 . 5; 

out=<lalpha  gain  zeta  omega  tau  cap  costf  costt> 
return 


’config.3’  -  (builds  the  basic  aircraft  transfer  function  by 

reading  in  the  appropriate  stability  derivatives 
and  using  them  to  build  "numac"  and  "denac") 

XU=-0 . 041; 

XW=0 .11; 

XQ=0; 

XDE=0 .0032; 

ZU=-0 . 25; 

ZW=-0. 87625; 

ZQ=0 ; 

ZDE=1 . 1; 

MU=0  ; 

MW=-0 . 04098; 

MQ=-2. 32375; 

MDE=0 . 33685; 

AN=MDE ; 

BN  *  XDE*MU  +  ZDE*MW  -  MDE* (XU  +  ZW)  ; 

CN  =  XDE* (ZU*MW-MU*ZW)  +  ZDE* (XW*MU-XU*MW)  +  MDE* (ZW*XU-XW*ZU) ; 
NUMAC 1=<AN  BN  CN>; 

TTH1=MIN (ABS (ROOT (NUMAC1) ) ) 

TTH2=MAX (ABS (ROOT (NUMAC 1) ) ) 

N1=<1  TTH1>; 

N2*<1  TTH2>; 

AC2 “CONVOLVE (N1 , N2 ) ; 

NUMAC=MDE *1.1 * AC2 
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Uo=205; 

W=25; 

G-32 .174; 

CTHO=COS (0 . 0785398) ; 

STHO=SIN (0.0785398)  ; 

AD  =  -MQ  -ZW  -ZU; 

BD  =  XU* (MQ+ZW)  +  MQ*ZW  -  Uo*MW  -  XW*ZU  +  W*MU; 

CD  =  -XU* (ZW*MQ-Uo*MW) +ZU* (XW*MQ+W*MW) -MU* (Uo*XW+W*ZW-G*CTHO) +G*MW*STH< 
DD  =  G*CTHO* (ZU*MW-MU*ZW)  -  G*XU*MW*STHO; 

DENAC=<1  AD  BD  CD  DD>; 

RTS=ROOT (DENAC) ; 

OMEGAP=MIN (ABS (RTS)  )  ; 

OME  GAS  P  =MAX (ABS (RTS) ) ; 

SIGMAP=MIN (ABS (REAL (RTS) ) ) ; 

SIGMASP=MAX (ABS (REAL (RTS) ) ) ; 

ZETAP=S IGMAP /OMEGAP ; 

ZETASP=S IGMASP /OMEGASP / 

OUTAC=<ZETAP  OMEGAP  ZETASP  OMEGASP>; 

ZETASP=0 .50; 

OMEGASP=3 .20; 

RETURN 


'nd.sof 1  -  (builds  the  second  order  filter  transfer  function, 
the  feel  system  and  actuator  transfer  functions, 
and  the  HOS  transfer  function  "numhos”  and  "denhos") 

kf ilt=81 ; 
zl=0 . 7 ; 
wnl=9; 

numf ilt=kf ilt ; 

denf ilt»<l  2*(zl)*(wnl)  wnl**2>; 

numf s=84 . 5; 

denfs=<l  31.2  67 6>; 

numserv=75**2; 

denserv=<l  105  75**2>; 

a=conv (numf i It , numf s) ; 

b=conv (numserv, numac)  ; 

nvomhos=conv  ( a,  b)  ; 

s=<l  0>; 

numhos=conv  (sf  nvimhos)  ; 

e=conv (denf ilt , denf s)  ; 

dl=<l  2* (zetap) * (omegap)  omegap**2>; 

d2=<l  2* (zetasp) * (omegasp)  omegasp**2>; 

denac=conv (dl , d2 ) ; 

f=conv (denserv, denac)  ; 

denhos*conv (e, f ) ; 

return 
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’Is. match’  -  (conducts  the  least  squares  matching  using  the  time 
response  of  the  HOS  and  outputs  the  LOES  transfer 
function  "num”  and  ”den”) 

tm=tm*10; 
tr=tr*10; 
tinc=5/tr; 
u (1) =0; 

for  i=l : tr,u (i+1) =u (i) +tinc; end 
for  i= (tr+2) :tm,u (i) =5 . 0;end 
<t , y>=lsim (numhos, denhos, u, .1); 

for  i=3 : tm, x (i-2, : ) =<y (i-1 )  y(i-2)  u(i-l)  u(i-2)>; 

for  i=3 :tm, capy (i-2, 1) =y (i, 1) ; 

theta=x\capy 

yout =x*thet a ; 

dnum2=<theta  (3)  theta  (4)  > 

dden2=<l  -theta (1)  -theta (2) > 

ns=2 ; 

<sd, ns>=sf orm (dnum2 , dden2 ) ; 

c=sd (3, : ) ; 

n=10; 

exec ( ' extract . ab' ) ; 
b=real (b) ; 
sc=<a, b; c>; 

<num, den>=tform (sc, ns) 

for  i=l : tm, u5 (i) =5; end 

<t , yls>=lsim (num, den, u5, .1); 

<t, yhos>=lsim (numhos, denhos, u5, .1) ; 
ycomp=<yhos  yls>; 
plot  (t, ycomp) 
return 


Note  that  ’extract. ab ’  and  ’roots’  are  identical  to  ’extract. ab; 2 ’ 
and  ’roots’  which  were  already  presented  as  subroutines  of  the 
LEAST. SQUARES  program.  Therefore  these  are  not  repeated  here. 


’tau’  -  (the  equivalent  time  delay  tau  is  calculated  by  iteration 
in  the  frequency  domain) 

<o, d,p>=bode (numhos, denhos, .3,10,21) ; 

<o, dl, pl>=bode (num, den, .3,10,21)  ; 
if  abs  (p(l) -pi  (1) )  >  180,... 

for  i=l : 21, pi (i) =pl (i) -360; end 
pplot=<p  pl>; 
plot (o,pplot,  '  logx' ) 
tau=*nt ; 
diff-1.0; 
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while  diff  >  0,  .  .  . 
sspl=0; . . . 
ssp2=0; . .  . 
plos=pl; . . . 

for  i=l : 21, plos (i) =plos  (i) -57 . 29578*0 (i) *tau; end,  . . . 
for  i=l :21, sspl®sspl  +  0 . 01745* (p (i) -plos (i) ) **2/ end,  ..  . 
tau=tau  +  0.001;... 
plos=pl; . . . 

for  i=l : 21 , plos (i) =plos (i) -57 . 29578*0 (i) *tau; end, . . . 

for  i=l  :21,  ssp2=*ssp2  +  0 . 01745* (p (i) -plos (i) ) **2; end, .. . 

diff=sspl-ssp2; 

pplot2=<p  plos>; 

plot (o,pplot2, ' logx' ) 

tau 

return 


’costf.time’  -  (calculates  cost^) 

<o,dbhos,phhos>=bode (numhos, denhos, .3,10,21); 

<o, dblos,phlos>=bode (num, den, .3, 10,  21)  ; 
if  abs (phhos (1) -phlos (1) )  >  180,... 

for  i=l : 21, phlos (i) =phlos (i) -360; end 
for  i=l :21, phlos (i) =phlos (i) -57 .29578*0 (i) *tau;end 
cost=0; 

for  i=l :21,  . .  . 

cost— cost+ (dbhos (i) — dblos (i) ) **2  +  0 . 01745* (phhos (i) -phlos  (i) ) **2 

costf= (20/21) *cost 

dcomp=<dbhos  dblos>; 

pcomp=<phhos  phlos>; 

return 


’costt.time'  -  (calculates  costt) 

<t,yhos>=lsim (numhos, denhos, u5, . 1) ; 
<t,ylos>=lsim(num,den,u5, .1) ; 
costt=0; 
for  i=l :tm, . . . 

costt— costt  +  (57 . 29578*yhos (i)  -  57 . 29578*ylos (i) ) **2 ; end 

costt«costt/tm; 

return 


The  next  program  (LS.FLT)  was  used  to  conduct  the  least 
squares  matching  using  the  actual  flight  test  (DAS)  data.  In  this 
case  the  input  vector  *’u”  and  the  output  vector  ”y"  are  not 
calculated  from  the  HOS  transfer  function  but  must  already  exist 
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in  the  MATRIXx  memory.  Getting  ”u"  and  ”y”  from  the  aircraft’s 
data  tape  into  the  MATRIXx  memory  was  a  trick  in  itself!  To  do 
this  consult  the  local  expert  on  the  data  acquisition  system  and 
computer  system  you’re  using! 

The  subroutines  in  LS.FLT  are  very  similar  to  the  subroutines 
already  presented  and  are  therefore  presented  below  without 
explanation.  Note  that  for  the  actual  flight  test  matching,  the 
HOS  was  not  modeled  in  the  frequency  domain  and  therefore  cost^  is 
not  calculated. 


LS.FLT 

INQUIRE  TM  'INPUT  MATCH  TIME  IN  SEC:' 

exec ( ' Is . f ltmat ' ) 

exec (' roots . '  ) 

exec ( ' costt . fit' ) 

v=205; 

cap= (omega**2) /4 . 5; 

out=<lalpha  gain  zeta  omega  0.000  cap  costt> 
return 


’Is. f ltmat ’ 
tm=tm*10 

for  i=3 : tm, x (i-2, : ) =<y (i-1)  y(i-2)  u(i-l)  u(i-2)>; 

for  i=3 : tm, capy (i-2, 1) =y (i, 1) ; 

theta=x\capy 

yout=x*theta; 

dnum2=<theta (3)  theta (4) > 

dden2=<l  -theta ( 1 )  -theta ( 2 ) > 

ns=2 ; 

<sd,  ns>=*sform  (dnum2,  dden2)  ; 

c=sd ( 3 ,  : )  / 

n*=10; 

exec ( '  ext  r act . ab ' ) ; 
b*real (b) ; 
sc*<a, b; c>; 

<num, den>*tf orm (sc, ns) 
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um=u ( 1 : tm, 1) ; 

<t, ylos>=lsim (num, den, urn, .1); 

yhos=y (1 : tm, 1)  ; 

ycomp=<yhos  ylos>; 

plot (t, ycomp) 

return 


’ roots ’ 

rd=root (den) ; 
sigma=abs (real (rd (1) ) )  ; 
omega=abs (rd(l) ) ; 
zeta=3igma/ omega; 
rn=root (num) ; 
lalpha=-rn; 
gain=num (1) ; 
return 


’costt. f It ' 

costt=0; 

for  i=l : tm, . . . 

costt=costt  +  (57 . 29578*yhos (i)  -  57 . 29578*ylos (i) ) **2; end 

costt=costt/tm; 

return 


The  final  program  presented  is  called  DATA. SMOOTH  and  was 
used  to  smooth  the  flight  test  data.  This  was  done  by  taking  50 
samples/second  and  averaging  every  five  data  points  to  get  10 
samples  a  second.  The  least  squares  program  worked  very  well  with 
only  10  samples/second  but  required  smooth  data.  Therefore  this 

program  was  used  4_\roughout  the  flight  test  matching. 

DIM 

INQUIRE  M  'ENTER  COLUMN  LENGTH:' 

PR=PITCHX01' ; 

SF-LNGLSB01' ; 

T-<0 : . 02 : (M* . 02) - . 02>' ; 

PLOT (T,SF, 'CHART  5  95  55  100') 

PLOT (T, PR, 'CHART  5  95  0  45') 

INQUIRE  T1  'ENTER  T1 : ' 

INQUIRE  T2  ' ENTER  T2 : ' 

PR2«PR(T1/.02:T2/.02,1)  ; 

SF2»SF(T1/.02:T2/ .02, 1)  ; 

T2«<0: . 02 : T2-T1>' ; 

PLOT (T2,SF2, 'CHART  5  95  55  100') 

PLOT (T2,PR2, 'CHART  5  95  0  45') 
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DIM 

INQUIRE  M2  ' ENTER  COLUMN  LENGTH  (DIV  BY  5) : ' 

FOR  1=1 : (M2/5) ,  .  .  . 

J=I*5; . . . 

U (1/ 1) = (SF2 ( J-4) +SF2 ( J-3) +SF2 ( J-2) +SF2 ( J-l ) +SF2 ( J) ) /5 ; . . . 
Y (I, 1) = (PR2 (J-4) +PR2 (J-3) +PR2 (J-2) +PR2 ( J-l) +PR2 ( J) ) / 5 ; END 
T3=<0: .1: (M2/5-1) * . 1>' ; 

PLOT (T3,  U, ' CHART  5  95  55  100') 

PLOT (T3,Y, 'CHART  5  95  0  45') 

ESF=U(1, 1) ; 

EPR=Y (1,1); 

FOR  1=1:  (M2/5)  ,  .  .  . 

U (I, 1) =U (1,1) -ESF; . . . 

Y (I,  1) =-0.0174533* (Y(1, 1) -EPR) ; END 
PLOT (T3,U, 'CHART  5  95  55  100') 

PLOT(T3,Y, 'CHART  5  95  0  45') 

RETURN 
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APPENDIX  C.  SORTIE  SUMMARY  AND  PROJECT  PILOT  EXPERIENCE 


Twenty-five  flights  totaling  27.8  hours  were  flown  in 
support  of  this  project.  A  buildup  approach  was  used  for 
familiarizing  the  team  with  the  task  and  test  aircraft.  Eight 
T-38A  sorties  were  flown  to  practice  the  offset  landing  task. 
Five  NT-33A  HUD  flights  were  then  used  for  familiarization  with 
the  test  aircraft  and  at  the  same  time  fullfill  a  TPS  curriculum 
requirement.  Lastly,  twelve  NT-33A  data  sorties  were  flown  in 
direct  support  of  this  project. 
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TABLE  Cl 


SORTIE  SUMMARY 


DATE 

hlQ 

CREW 

MISSION 

12 

SEP 

89 

T-38 

LINDSEY/RUNYON 

OFFSET  LAND 

13 

SEP 

89 

T-38 

THOMAS /RUNYON 

OFFSET  LAND 

13 

SEP 

89 

T-38 

BAUM/RUNYON 

OFFSET  LAND 

15 

SEP 

89 

T-38 

LINDSEY/RUNYON 

OFFSET  LAND 

18 

SEP 

89 

T-38 

THOMAS /RUNYON 

OFFSET  LAND 

19 

SEP 

89 

NT-33 

THOMAS /EASTER 

HUD  EVAL 

20 

SEP 

89 

NT-33 

THOMAS /EASTER 

DATA  1 

20 

SEP 

89 

T-38 

BAUM/RUNYON 

OFFSET  LAND 

20 

SEP 

89 

NT-33 

BAUM/EASTER 

HUD  EVAL 

21 

SEP 

89 

NT-33 

THOMAS /EASTER 

DATA  2 

21 

SEP 

89 

NT-33 

LINDSEY/EASTER 

HUD  EVAL 

22 

SEP 

89 

NT-33 

BAUM/ EASTER 

DATA  3 

22 

SEP 

89 

NT-33 

MANNING /EASTER 

HUD  EVAL 

25 

SEP 

89 

NT-33 

LINDSEY/EASTER 

DATA  4 

25 

SEP 

89 

NT-33 

LIU/EASTER 

HUD  EVAL 

25 

SEP 

89 

T-38 

MANNING/RUNYON 

OFFSET  LAND 

26 

SEP 

89 

NT-33 

LINDSEY/EASTER 

DATA  5 

3 

OCT 

89 

NT-33 

BAUM/BALL 

DATA  6 

3 

OCT 

89 

NT-33 

THOMAS /BALL 

DATA  7 

4 

OCT 

89 

NT-33 

LINDSEY/BALL 

DATA  8 

4 

OCT 

89 

NT-33 

LINDSEY /BALL 

DATA  9 

5 

OCT 

89 

NT-33 

BAUM/BALL 

DATA  10 

10 

OCT 

89 

NT-33 

BAUM/ BALL 

DATA  11 

10 

OCT 

89 

NT-33 

THOMAS /BALL 

DATA  12 

16 

OCT 

89 

T-38 

LIU/RUNYON 

OFFSET  LAND 

TOTAL  SORTIE  COUNT 

SORTIES  HOURS 

T-38  TRAINING  8  9.3 

NT-33  TRAINING  5  7.8 

NT-33  DATA  12  1Q_uZ 

TOTAL  HOURS  27.8 


HOURS 

1.2 

1.2 

1.2 

1.1 

1.2 

1.6 

1.0 

1.2 

1.4 
0.9 

1.5 

0.8 

1.6 
1.0 
1.7 
1.0 
0.9 
0.8 
0.8 
0.9 
1.0 
0.8 
0.9 
0.9 
1.2 
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TABLE  C2 

PROJECT  PILOT  EXPERIENCE 


Pilot 

Aircraf  t 

Hours 

A 

C-  141 

2500 

B 

F/RF-4 

1000 

T-39 

50 

B-52 
T-  37 


2200 

150 


APPENDIX  D.  AIRCRAFT  DESCRIPTION  AND  TEST  INSTRUMENTATION 


The  NT-33A  variable  stability  airplane,  shown  in  Figure  Dl, 
is  a  T-33  jet  trainer  modified  with  a  Variable  Stability  System. 
The  VSS  can  be  divided  into  two  independent  parts.  The  first 
part,  the  variable  feel  system,  provides  a  variety  of  stick  and 
rudder  pedal  forces,  gradients,  and  displacements.  The  variable 
feel  is  provided  by  disconnecting  the  elevator,  aileron,  and  rud¬ 
der  controls  in  the  front  cockpit  from  their  respective  control 
surfaces  and  connecting  the  controls  to  separate  servomechanisms. 
The  second  part  of  the  VSS  is  the  response  feedback  flight  con¬ 
trol  system.  This  part  augments  the  normal  T-33  dynamics  to  rep¬ 
resent  those  of  the  vehicle  being  simulated. 


Figure  Dl .  NT-33A  Variable  Stability  Aircraft 

The  augmentation  is  accomplished  by  connecting  the  elevator, 
aileron,  and  rudder  control  surfaces  to  individual  servos .  These 
individual  servos  can  be  driven  by  a  number  of  different  inputs, 
such  as  the  aircraft's  artificial  feel  system  (pilot's  commends, 
position  or  force) ,  attitude  and  rate  gyros,  accelerometers, 
dynauaic  pressure  pickups,  angle  of  attack  vane  and  sideslip 
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n 


probe.  This  arrangement,  through  a  response-feedback  system,  al¬ 
lows  the  normal  t-33  derivatives  to  be  augmented  to  simulate  the 


handling  qualities  of  existing  or  hypothetical  aircraft.  A  block 
diagram  of  the  VSS  is  shown  in  Figure  D2 . 


Figure  D2 .  Variable  Stability  NT-33A  Block  Diagram 


The  original  T-33  nose  section  has  been  replaced  with  the  larger 
nose  of  an  F-94  to  provide  the  volume  required  for  the  electronic 
components  of  the  response- feedback  system  and  recording  equip¬ 
ment  .  The  physical  layout  of  the  control  system  is  shown  in 
Figure  D3.  Each  control  surface  has  an  electro-hydraulic  posi¬ 
tion  servo  which  is  actuated  by  inputs  from  the  VSS.  The  servos 
operate  in  parallel  with  the  normal  T-33  control  surface's  ac¬ 
tuating  mechanisms.  Each  of  the  surface  position  actuators  has  a 
hydraulic  limiting  circuit  which  limits  the  maximum  hinge  moment 
which  can  be  generated  by  the  flight  control  system.  The  rear 
cockpit  controls  have  a  direct  mechanical  connection  to  the 
aircraft  control  surfaces  at  all  times. 
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NT- 33  SAFETY  FEATURES 


The  primary  safety  feature  in  the  airplane  is  the  safety 
pilot;  however,  the  NT-33A  also  has  numerous  safety  provisions 
designed  into  the  variable  stability  equipment. 

Control  Interlock  System:  The  control  buttons  and  switches 
for  the  various  functions  of  the  VSS  are  wired  so  that  the  proper 
sequence  of  operation  of  these  controls  must  be  observed  to 
energize  the  various  parts  of  the  VSS.  No  action  will  result 
from  a  button  or  switch  activated  out  of  sequence.  For  example, 
interlocking  circuits  prevent  servo  engagement  prior  to  auto 
balance  engage. 

Automatic  Safety  Trips :  The  automatic  safety  trip  monitors 
the  servovalve  amplifier  error  signals  and  the  normal  and  lateral 
accelerometer  output  signals.  If  these  signals  exceed  preset 
values,  the  VSS  is  automatically  shut  off.  Safety  trip 
accelerometer  limits  have  been  set  as  follows: 

a.  n2  Pushover  (-  0.3g  on  the  g  meter) 

b.  nz  Pullout  (+  4 . 8g  on  g  meter) 

c.  ny  <+  0.25g) 

Audio-visual  VSS  Shut-off  Warning  System:  When  the  VSS  has 
been  disengaged  either  automatically  or  manually,  red  lights  will 
flash  in  both  cockpits  and  a  "beep,  beep"  will  be  heard  in  the 
interphone . 

Special  Aircraft  Limitation:  The  variable  stability  NT-33A 
is  limited  to  375  KIAS  with  a  never  exceed  speed  of  400  KIAS. 
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Special  Pilot  Emergency  Procedures:  In  addition  to  normal 
NT-33A  emergency  procedures,  the  following  pertain  specifically 
to  the  variable  stability  NT-33A: 

-Manual  VSS  trip:  The  project  pilot  can  manually  disengage 
the  VSS  from  the  control  surfaces  and  return  control  of  the 
aircraft  to  the  safety  pilot .  Disengage  switches  are  located  on 
both  the  centerstick  and  the  sidestick. 

-In  the  event  of  safety  pilot  incapacitation,  the  project 
pilot  cam  fly  the  aircraft  back  to  the  base  via  his  fly-by-wire 
controls  with  normal  T-33  characteristics.  This  is  accomplished 
by  actuating  the  red  guarded  safety  trip  bypass  switch  located  on 
the  left  side  of  the  VSS  engage  panel  and  sequentially  depressing 
the  adjacent  four  buttons  starting  from  the  left.  Subsequent 
buttons  are  pressed  after  the  light  below  the  previously  pressed 
button  is  lit. 

-If  a  feel  system  hardover  should  occur,  the  project  pilot 
can  activate  the  feel  system  hydraulic  bypass  switch  and  move  the 
control  stick  out  of  the  way  to  ensure  non-interference  during 
ejection . 
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TABLE  D1 


NT-33A  DIGITAL  TAPE  PARAMETERS 

DIGITAL 

CHANNEL 

NUMBER 

RECORDED  VARIABLE 

ENGINEERING  UNITS 

1 

Pressure  Altitude 

feet 

2 

Normal  Acceleration 

g ’s 

3 

Velocity  (Indicated) 

knots 

4 

Pitch  Rate 

degrees /second 

5 

Pitch  Angle 

degrees 

6 

Yaw  Rate 

degrees /second 

7 

Elevator  Stick  Deflection 

inches 

8 

Angle  of  Sideslip 

degrees 

9 

Event  Mark 

N/A 

10 

Radar  Altitude 

feet 

11 

Pitch  Error 

degrees 

12 

Roll  Rate 

degrees /second 

13 

Roll  Angle 

degrees 

14 

Longitudinal  Acceleration 

g’s 

15 

Roll  Error 

degrees 

16 

Elevator  Deflection 

degrees 

17 

Lateral  Acceleration 

g’s 

18 

Elevator  Stick  Force 

pounds 

19 

Vertical  Velocity 

f eet/second 

20 

Rudder  Deflection 

degrees 

21 

Total  Aileron  Deflection 

degrees 

22 

Change  in  Heading 

degrees 

23 

Lateral  Stick  Deflection 

inches 

24 

Angle  of  Attack 

degrees 

25 

Aileron  Stick  Force 

pounds 

26 

Rudder  Pedal  Deflection 

inches 

27 

Rudder  Pedal  Force 

pounds 

28 

Time  Rate  of  AOA  Change 

degrees /second 

APPENDIX  E.  AIRCRAFT  CONFIGURATION  IDENTIFICATION 


The  plots  in  this  appendix  demonstrate  the  fidelity  between 
the  computer  generated  data  and  the  actual  flight  test  data. 

That  is,  these  plots  show  how  close  the  theoretical  data  (i.e.  the 
estimate  of  how  the  aircraft  will  fly)  agrees  with  how  the 
aircraft  actually  flies.  This  is  important  since  the  prediction 
of  flying  qualities  and  the  verification  of  the  Least  Squares  LOES 
program  were  dependent  on  the  analytical  data. 

The  CALSPAN  Corporation  accomplished  configuration 
verifications  for  the  three  baseline  configurations  (1-1,  2-1,  and 
3-1).  These  verifications  are  shown  in  Figures  E1-E3.  The  first 
plot  on  these  pages  shows  the  theoretical  input  (8e  in  degrees) 
compared  to  the  actual  input  (Fs  in  pounds).  The  next  two  plots 
show  the  corresponding  theoretical  and  actual  responses  in  pitch 
rate  and  angle  of  attack. 

The  remaining  configurations  were  verified  as  part  of  the 
flight  test  program  and  are  shown  in  Figures  E4-E16.  The  first 
plot  in  these  figures  shows  the  actual  stick  force  input.  The 
second  plot  shows  the  corresponding  actual  pitch  rate  output 
compared  to  the  theoretical  pitch  rate  output.  The  theoretical 
output  in  these  figures  was  computed  using  the  actual  input.  This 
was  possible  because  MATRIXx  (31)  can  accept  arbitrary  inputs. 

Note  that  the  agreement  between  the  theoretical  and  actual 
aircraft  responses  was  very  good  for  both  the  CALSPAN  and  MATRIXx 
verifications.  This  indicates  that  the  analytical  models  were 
adequate. 
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ANGLE  OF  ATTACK  PITCH  RATE  8«  (DEG)  or  Fmm  (LBS) 


F LT  445?  KNOTTS  PRRRRG  9/15/89  RFTPS  CRL 
REC  9  DER  =  4 9 0 ,  DERD=450,  DEQ=500  (  1-1) 


0.000  DEG'SEC  0.000  RfiD/SEC 


I  !  i  i  i  i  i  i  i  '  I  i  i  i  i  i  i  ;  *  i  1  i  . I  i  i  :  i  i  i  i  i  i  i  i  i  ■;  j  i  i  i  i 

3  12  3  4 

TIKE  (SEC) 


‘ire  El. 


System  Verification  by  CALSPAN 


Configuration  1-1 
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PITCH  RATE  (RPS) 


TIME  (SEC) 

Figure  E4.  System  Verification  -  Configuration  1-1 
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PITCH  RATE  (RPS) 


Configuration  1-10  was  not  flown  in  the  landing 
evaluation  because  it  was  highly  uncontrollable 
at  altitude.  Therefore,  there  was  no  data  collected 


for  system  verification 


Figure  E6.  System  Verification  -  Configuration  1  10 
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PITCH  RATE  (RPS) 
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Figure  E7.  System  Verification  -  Configuration  2-1 
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0) 
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Figure  E8.  System  Verification  -  Configuration  2-D 
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Figure  E9.  System  Verification  -  Configuration  2-2 
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PITCH  RATE  (RPS)  Fes  (LBS) 


2 


Figure  E10.  System  Verification  -  Configuration  2-5 
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Figure  Ell.  System  Verification  -  Configuration  2-7 
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PITCH  RATE  (RPS)  Fes  (LBS) 


TIME  (SEC) 


TIME  (SEC) 

Figure  E12.  System  Verification  -  Configuration  3-1 


187 


PITCH  RATE  (RPS) 


0  12  3  4  5 

TIME  (SEC) 


Figure  E13.  System  Verification  -  Configuration  3  3 
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0  1 


0  1 


Figure  E14.  Syste 
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Theoretical  Response 
Actual  Response 
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Verification  -  Configurati 


PITCH  RATE  (RPS) 


PITCH  RATE  (RPS) 


APPENDIX  F.  AIRCRAFT  PARAMETERS  AND  STABILITY  DERIVATIVES 


This  appendix  contains  the  NT-33A  longitudinal  and  lateral 
parameters,  feel  system  and  actuator  parameters,  and  stability 
derivatives.  Table  FI  shows  the  longitudinal  and  lateral 
dynamics.  Note  that  wap  and  £ap  were  varied  to  produce  the 
different  flying  qualities.  The  only  other  variable  between 
configurations  was  the  pre-filter  dynamics  (shown  in  Chapter  V, 
Table  VIII).  Table  FI  also  shows  the  feel  system  and  actuator 
dynamics  which  were  held  constant  throughout  the  flight  test 
program.  Table  F2  shows  the  stability  derivatives  for  the  three 
baseline  conf igurations  (1-1,  2-1,  and  3-1). 
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TABLE  FI 


NT-33A  PARAMETERS 

PARAMETER 

VALUE 

unsp  (rad/sec) 

variable 

^sp 

variable 

n/a  (g/rad) 

4.50 

1/Tq2  (1/sec) 

0.70 

u)np  (rad/sec) 

0.  17 

0.  15 

1/T01  (1/sec) 

0.08 

u)ncj  (rad/sec) 

1.30 

«d 

0.20 

<p/B 

1.50 

Tr  (sec) 

0.30 

Fes/in  (lbs/in) 

6.50 

Fas/in  (lbs/in) 

3.00 

A 

Feel  systems:  -  (in/lb) 

s2  +  2( . 6 ) ( 26 )s  +  262 


Elevator:  A  =  84.50 
Aileron:  A  *  169.00 
Rudder:  A  =  11.47 


752 

Actuators:  -  (deg/in) 

s2  +  2( .7)(75)s  ♦  762 
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Table  F2 


NT-33A  Flight  Test  Stability  Derivatives 


Paraneter 

1-1 

2-1 

3-1 

(■> 

8P 

1.00 

2.00 

3.20 

0.75 

0.75 

0.50 

Xu 

-0.041 

-0.041 

-0.041 

xw 

0.11 

0.11 

0.11 

x. 

0 

0 

0 

xa 

0.0032 

0.0032 

0.0032 

e 

Zu 

-0.25 

-0.25 

-0.25 

Zw 

-0.74939 

-0.77925 

-0.87625 

0.0 

0.0 

0.0 

l6 

1.1 

1.1 

1.1 

e 

"a 

0.0 

0.0 

0.0 

", 

-0.002134 

-0.01128 

-0.04098 

M, 

-0.75061 

-2.22075 

-2.32375 

M* 

0.33685 

0.33685 

0.33685 

APPENDIX  G.  LEAST  SQUARES  LOES  MATCHING  RESULTS 


This  appendix  contains  the  Least  Squares  matching  results 
from  the  flight  test  program.  The  analytical  (pre-flight)  results 
are  shown  in  Figures  G1-G13.  These  results  were  obtained  using 
computer  generated  HOS  data  based  on  the  known  NT-33A  dynamics. 

The  HOS  and  LOES  time  responses  were  then  matched  using  the  Least 
Squares  program.  The  first  plot  on  each  page  shows  a  comparison 
between  these  time  responses.  The  next  two  plots  show  a 
comparison  between  the  HOS  and  LOES  frequency  responses.  The 
corresponding  LOES  parameters  and  cost  functions  are  tabulated  at 
the  bottom  of  each  page.  Also,  the  predicted  handling  quality 
level  (based  on  MIL-STD-1797  guidance  for  equivalent  parameters) 
is  listed  for  each  configuration. 

The  flight  test  results  (Figures  G14-G25)  differ  from  the 
analytical  results  in  that  the  HOS  pitch  rate  response  was  taken 
from  actual  flight  test  data  instead  of  computer  generated  data. 
The  Least  Squares  program  was  then  used  to  match  the  HOS  and  LOES 
time  responses  (shown  on  the  first  plot  for  each  configuration). 
The  next  two  plots  show  how  the  ’’actual”  time  delay  (referred  to 
as  TD  to  distinguish  it  from  equivalent  time  delay  Te)  was 
measured.  The  frequency  responses  (and  therefore  tq)  were  not 
measured  during  the  actual  flight  test  program.  A  detailed 
comparison  between  the  analytical  and  flight  test  results  is  given 
in  Chapter  V. 
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TIME  (SEC) 


FREQUENCY  (RPS) 


CONFIG  LEVEL  l'T92  K®  t-p  w0p  tj  CAP  coatf  coat^ 

( 1/aec ) ( rpa/ lb )  (rpa)  (aac)  <l/*aec) 

1-3  2  0.706  0.007  0.396  0.769  0.185  0.131  156.25  0.062 


Figure  G2. 


Least  Squares  LOES  Matching  Results  (Analytical  Data) 

Configuration  1-3 


PHASE  (DEG) 


TIME  (SEC) 


TIME  (SEC) 


FREQUENCY  (RPS) 


CONFIG  LEVEL  l/Te2  K*  t-p  “«p  T*  CAP  eomtg  co*fct 
(l/>*c)(rpi/lb)  (rpa)  <aac>  (l/«a«cl 

2-1  1  0.690  0.015  0.574  1.751  0.044  0.662  74.33  0.006 

Figura  G4.  Laaat  Squaraa  LOE3  Matching  Raaulta  (Analytical  Data) 

Configuration  2-1 
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1 


1/Ti2  K,  Zap 
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0.709  0.013  0.543 


wap  CAP  coitf  coatt 

(rpa)  lame)  (1/gaac) 

1.609  0.070  0.575  108.24  0.009 


Figure  G5. 


Least  Square*  LOES  Matching  Result*  (Analytical  Data) 

Configuration  2-D 
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PHASE  (DEG) 


TIME  (SEC) 
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0 

-100 


1  10 
FREQUENCY  (RPS) 

1/T*2  k*  r.p  ▼#  CAP  cost f  coatt 

( 1/aac )( rpa/ lb )  (rpa)  (aac)  (1/gaac) 

0.715  0.011  0.471  1.518  0.107  0.512  186.35  0.019 

Figura  G6.  Laaat  Squaraa  LOBS  Matching  Raaulta  (Analytical  Data) 
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Tigura  C7,  Laaat  Squaraa  LOES  Hatching  Raaulta  (Analytical  Data) 

Configuration  2-5 
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PHASE  (DEG) 


Mi  i-i  1  k  A_A.i.-i..I  n,  U 

012345678 
TIME  (SEC) 


FREQUENCY  (RPS) 
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FREQUENCY  (RPS) 


CONFIG  LEVEL  1/T#Z  rap 
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2-7 
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Figur*  G8.  Laaat  Squaraa  LOES  Matching  Raaulta  (Analytical  Data) 

Configuration  2-7 
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CONFIG  LEVEL  1/762  *6  tap  “*p  T0  CAP  coat  •  coat- 
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Figure  G9.  Laaat  Squaraa  LOSS  Matching  Raaulta  (Analytical  Data) 

Configuration  3-1 
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Figura  G10.  Laaat  Squaraa  LOES  Matching  Raaulta  (Analytical  Data) 

Configuration  3-3 
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Figura  Gil.  Laaat  Squaraa  LOES  Matching  Raaulta  (Analytical  Data) 

Configuration  3-5 
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Figure  G13.  Laaat  Squaraa  LOES  Hatching  Raaulta  (Analytical  Data) 
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Flours  G15.  Least  Squares  LOES  Matching  Results  (Flight  Test  Data) 
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figura  G18.  Laaat  Squaraa  LOBS  Matching  Raaulta  (flight  Taat  Data) 
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Figura  021 .  Laaat  Squaraa  LOES  Matching  Raaulta  (Flight  Taat  Data) 
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Figure  G25 .  Laaat  Squares  LOES  Matching  Results  (Flight  Test  Data) 
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APPENDIX  H. 


BANDWIDTH  THEORY  AND  PREDICTION  RESULTS 


The  Bandwidth  method  has  been  proposed  as  a  handling 
qualities  requirement  by  MIL-STD-1797  (3).  The  Bandwidth 
method  is  simple  to  use  because  it  assumes  a  ’’gain  only”  pilot 
model  and  involves  only  the  use  of  open  loop  pitch  to  stick 
force  (9/Fs)  Bode  plots. 

Bandwidth  can  be  loosely  defined  as  the  maximum  frequency 
at  which  closed  loop  compensatory  tracking  can  take  place 
without  threatening  the  stability  of  the  aircraft;  i.e.  the 
maximum  open  loop  crossover  frequency.  Hence,  a  large  value  of 
bandwidth  is  generally  desireable  to  achieve  superior  tracking 
performance. 

The  reason  for  including  bandwidth  in  this  study  was  to 
compare  the  results  of  the  bandwidth  predicted  handling  quality 
levels  with  the  least  squares  LOES  results.  The  bandwidth 
method  sets  up  boundaries  for  Level  1,  2,  and  3  handling 
qualities,  and  thus  is  well  suited  for  a  comparison  with  the 
LOES  method.  The  bandwidth  criteria  are  based  on  maximum 
crossover  frequency  and  system  phase  delay,  while  the  LOES 
method  uses  requirements  for  equivalent  £sp,  wsp,  and  tq.  The 
approach  taken  in  this  study  was  to  see  how  well  the  bandwidth 
boundaries  and  the  phase  delay  parameter  correlated  with  the 
LOES  boundaries  and  equivalent  time  delay. 

The  following  discussion  of  the  bandwidth  theory  is  taken 
from  MIL-STD-1797  (3)  and  Hodgkinson  (32).  Crossover 
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frequency,  directly  determined  by  pilot  gain,  is  a  rough 
measure  of  the  rapidity  of  the  closed  loop  response. 
Physically,  the  pilot  will  increase  his  gain  to  track  more 
rapidly  moving  targt  with  acceptable  error.  However,  the 
pilot  cannot  indefinitely  increase  crossover  frequency  by 
increasing  gain,  because  he  will  eventually  lose  closed  loop 
stability  (when  the  phase  margin  of  the  open  loop  system 
becomes  negative).  The  pilot  would  like  to  choose  a  value  of 
crossover  frequency  which  allows  him  to  double  his  gain  and 
still  provide  adequate  phase  margin.  A  reasonable  crossover 
frequency  would  then  be  one  which  provides  at  least  6  dB  of 
gain  margin  and  45  degrees  of  phase  margin. 

The  above  crossover  frequency  is  the  bandwidth  frequency 
(ojbw)  and  is  shown  in  Figure  HI.  To  find  wBW  for  a  system, 
first  find  the  frequency  at  which  the  phase  shift  is  -135 
degrees;  this  is  “BWphase*  Then  find  the  amplitude  at  which 
the  phase  shift  is  -180  degrees  and  add  6  dB.  The  frequency 
corresponding  to  this  amplitude  is  wByga^n.  The  smaller  of 
these  two  frequencies  is  uigy. 

Handling  qualities  and  pilot  ratings  are  not  dependent  on 
bandwidth  alone;  the  shape  of  the  phase  curve  at  frequencies 
above  wgy  becomes  important  as  well.  If  the  phase  curve  drops 
off  rapidly  at  frequencies  above  wgy,  the  aircraft  will 
generally  receive  poor  pilot  ratings,  since  an  abrupt  loss  in 
stability  margin  is  produced  when  the  pilot  attempts  to 
increase  the  crossover  frequency.  One  measure  of  rapid  phase 
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equivalent  system  time  delay,  unlike  wBW,  cannot  easily  be 
measured  by  hand.  Phase  delay  (Tp),  a  parameter  that  measures 
phase  rolloff  by  hand,  is  defined  in  Figure  HI.  Usually,  Tp  is 
numerically  similar  to  Te. 


Frajjency  (r&v'sec) 

UBW  =  min C “BWphase 1  wBWgain^ 

Tp  =  -  [<p(  )+180]  /  C57 . 3x2^30] 


Figure  HI.  Definition  of  Bandwidth  and  Phase  Delay  Parameters 


The  bandwidth  criteria  suggests  that  systems  with  high 


attainable  crossover  frequencies  and  without  rapid  phase 
rolloffs  should  have  good  handling  qualities.  Figure  H2  shows 
flying  qualities  boundaries  based  on  bandwidth  frequency  and 
Tp.  The  boundaries  are  referred  to  as  the  bandwidth  criteria. 


b)  Category  C  Flight  Phases 
Figure  H2.  Bandwidth  Requirements  ( MIL-STD-1797,  Ref  HI) 


The  bandwidth  method  can  then  be  summarized  as  follows: 

1.  Determine  w18q  from  the  Bode  plot  of  0/ Fg. 

2.  Find  wBWphase  ( f req  where  phase  margin  is  45  degrees). 

3.  Find  magnitude  0/Fg  at  and  add  6  dB. 

4.  Find  (freq  where  above  magnitude  occurs). 

5.  wgy  =  win  [wgyp^ase’  wBWgain-*' 

6.  Tp  -  -  [<p(  2wj.qo  )  +  180  ]  /  [  57 . 3x2u>280  3  - 

7.  Find  Tp  and  wBW  on  the  Tp  vs  wBy  plot  to  determine 

the  predicted  level  of  handling  qualities. 
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Figure  H3.  Bandwidth  Results  for  Configuration  1-1 
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Figure  H4.  Bandwidth  Results  for  Configuration  1-3 
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Figure  H5.  Bandwidth  Results  for  Configuration 


Figure  H6 .  Bandwidth  Results  for  Configuration  2-1 
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Predicted  Flying  Qualities  -  Level  2 


Figure  H7 .  Bandwidth  Results  for  Configuration  2-D 
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Figure  H8.  Bandwidth  Results  for  Conf iguration  2-2 
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Figure  H9.  Bandwidth  Results  for  Configuration  2-5 


231 


FREQUENCY  (RPS) 


wBW  s  min[wBWphaae,  wBWgain3  =  2.00  rps 

t-  *  - Co( 2w^gQ ) +180) / [57 . 3x2w^gQ] 

=  - [-235+ 180] / [2( 3. 25 ) ( 57 .3)1 
=  0.148  sec 

Predicted  Flying  Qualities  -  Level  2 


Figure  H10.  Bandwidth  Results  for  Configuration  2-7 
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Figure  H12.  Bandwidth  Results  for  Configuration  3-3 
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Figure  H14.  Bandwidth  Results  for  Configuration  3-6 
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Figure  H15.  Bandwidth  Results  for  Configuration  3-8 
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The  most  widely  accepted  method  for  specifying  flying  qualities 
of  highly-augmented  aircraft  is  the  low-order  equivalent  systems 
technique.  This  technique  matches  the  frequency  response  of  a  high- 
order  system  (HOS)  with  a  lower-order  equivalent  system  (LOES).  The 
LOES  is  generally  in  the  form  of  a  “classical"  (unaugmented) 
aircraft  so  that  comparisons  with  the  classical  data  base  can  be 
made.  Most  LOES  research  has  been  accomplished  using  frequency  response 
matching,  however  there  are  also  methods  available  for  matching  the 
time  responses  of  the  HOS  and  LOES.  Time  response  matching  is  an 
attractive  option  for  flight  test  applications  since  all  that  is  required 
is  output  versus  input  data. 

This  thesis  develops  and  tests  a  Least  Squares  LOES  time  response 
matching  program.  The  program  was  tested  analytically  at  the  Air 
Force  Institute  of  Technology,  Wright-Patterson  AFB,  OH.  During  the 
analytical  testing,  a  HOS  was  modeled  using  an  A-4D  aircraft  with 
servo-actuator  and  feel  system  dynamics.  Lower  order  equivalent 
parameters  were  calculated  using  the  Least  Squares  program  and  the 
results  were  compared  to  those  obtained  using  a  frequency  response 
matching  program  (LONFIT).  Various  flight  conditions  and  feel  system 
dynamics  were  used  to  change  the  HOS  dynamics  and  the  effects  of  these 
changes  on  the  two  matching  techniques  were  investigated. 

The  Least  Squares  LOES  program  was  also  tested  during  a  flight 
test  program  the  the  USAF  Test  Pilot  School,  Edwards  AFB,  CA.  During 
this  test,  a  variable  stability  NT-33A  aircraft  operated  by  the 
CALSPAN  Corporation  was  used  to  model  the  higher-order  aircraft 
dynamics.  Thirteen  configurations  were  tested  using  system  identification 
and  offset  landing  tasks.  The  configurations  were  chosen  to  span  Level  1, 
2,  and  3  flying  qualities  by  varying  the  short  period  natural  frequency, 
the  short  period  damping  ratio,  and  pre-filter  dynamics.  The  LOES 
Least  Squares  program  was  used  to  predict  the  flying  qualities  of  the 
configurations  based  on  MIL-STD-1797  guidance  for  equivalent  time  delay, 
short  period  natural  frequency,  and  short  period  damping  ratio.  Parameter 
identification  was  also  accomplished  using  the  actual  flight  test  output 
versus  input  data.  The  results  of  the  flying  qualities  predictions  were 
compared  to  predictions  given  by  Hoh's  Bandwidth  method  for  specifying 
flying  qualities. 


